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ABSTRACT 


We have given below an operational method for solving the equation of radiative transfer in a “‘gray” 
plane-parallel atmosphere. We have assumed an expansion for the Ergiebigkeit, B(r), of the form: 


= o +r +> AKuds) , 
0 


where K,,(7) is the exponential integral of the nth order, defined by 


K,(z) -f — dx. 
1 x 


The coefficients, A;, to any required degree of approximation, can be evaluated from a system of linear 
equations. We have obtained solutions according to various independent assumptions and find the opera- 
tional method to be at once simple and capable of a high degree of accuracy. For instance, with four 
terms in the expansion above, we have reached an accuracy of two parts in a million. A comparison with 
Hopf’s exact solution shows that our method gives consistent accuracy both near the boundary and in 
the deeper layers of the atmosphere, a feature not reached by any extant approximate method, with the 
possible exception of the Ambarzumian-Chandrasekhar method of solution by functional equation. The 
latter method has, however, so far been applied only to the emergent intensity at the boundary of the 


star. 
We believe that our operational method is capable of extension to more complicated radiative transfer 


problems, e.g., to the case of nongray atmospheres. 

The equations describing the distribution of thermal neutrons in a pile are identical 
with those that treat the intensity of radiation in a star. As a result of the rekindled in- 
terest in the problem, several new methods of attack have become available. The pro- 
cedure given herewith has the double advantage of simplicity and high accuracy. 

The basic equation of radiative transfer, for material stratified in planes, is 


(1) 


where 7 is the optical depth, u the cosine of the angle between the ray and the outward 
normal, / the specific intensity, and B the emission function. / is a function of 7 and y; 


1 


4 
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B depends only on +r. For radiative equilibrium, the absorbed and emitted radiations 
must be equal, or 


1 1 
an f an f = 2) 


From equations (1) and (2) and certain boundary conditions we are to determine both 
I and B. The chief boundary condition is that the intensity of the ingoing beam vanish 
for r = 0, at the outer surface of the slab. Introducing the symbolic operator 

a 


OT 
we obtain the symbolic linear equation, 


(1-uD)I=B, 


whose symbolic solution, by the ordinary methods of differential equations, is 


1 
(5) 


The first term on the right represents the general solution and the second the particular 
solution. 

The constant of integration, A, differs for the outgoing and ingoing beams. Denote 
the former by /* and the latter by /~. At some very deep level, 7 = 11, let 7 = 7}. Then 


1 
A 5? |... , (6) 


Similarly, for the ingoing beam, / = 0 for + = 0, or 


1 


Now, unless /* or B increases exponentially with 7, a condition that is unacceptable, we 
must have A+ — 0 as 7;—> ©, for a slab of infinite thickness. Thus we get for our two 


equations: 


Substitute these equations in equation (2). Then: 


fred + f ras Bas 2 f Bay 


In the second integral on the left-hand side let us replace » by —u. Then we obtain 


a symbolic integrodifferential equation for B. The quantity in brackets is a function of 
u alone. The definite integral on the right is, therefore, a function only of 7, which we in- 


> 
(4) 
- 
| 
= 
51 [t= B 
| 
(8) be 
: 
| 
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dicate as f(r). The operator integral on the left is a function only of D. Hence equation 
(11) is really in the familiar symbolic form: 


F(D)B= f(r). (12) 


We obtain our general solution by setting f(r) = 0 and trying 
B= e™. (13) 


Since 
F(D) = (m), (14) 


1 
F(m) = f in| tanh m=0. (15) 


The equation, 
tanh m=m, (16) 


has only two roots that will satisfy condition (15). There is an infinite number of conju- 
gate imaginary roots, but these would correspond to values of @ greater than 7 and hence 


can be disregarded. 
The two significant roots are coincident at m = 0. Hence the linear solution: 


B=a+6r, 


to which we must add the particular solution, 


1 1 
We have changed the variable of integration from yp to y, in order to avoid confusion. 


The bracketed factor is a function of y, which we may call ¥(y). Hence F (D) operates 
only on the exponential. Thus, setting m = —1/y, in equation (15), we have 


¥ (y) 


In the present paper we shall not make any special use of the form of the denominator. 


After all, ¥(y) depends in a peculiar way upon B, whose value we still do not know. The 
important conclusion is that we may write T in the form 


1 
J, (y) dy, 


and let the undetermined function, ¢, carry the burden. Set 


1 
B=atbrt e-$(y) dy 
0 
back in equation (10). Then 


=0. 
0 


(22) 


The left-hand side of this equation must vanish for all values of r. 
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Before we attempt solution of these equations, let us set equation (21) back in our 
symbolic equations (8) for 7. Then 


1 
=a+br+ bp +f (y) dy, 


(24) 


In equation (24) we have changed yw to —yu, so that 0 < w < 1 for both equations. Now 


calculate the integral: 
L = (25) 


L(n) = f+ (= (26) 


The factor (— 1)" enters because we changed the sign of « for the negative beam. The 
result is 


a +f ¢ (y) dy | 


L(0) =2J =2B, 


= Constant , 


By definition, 


L(2) =2K, 
where J and K are the conventional parameters and 7¥ is the flux. 


Break the first integral on the right-hand side into partial fractions, dividing the 
terms in yw” by the respective denominators. The result is 


tf 
l=1 


+ (= 1) = 1) | (y) dydu +0(n), 
1 


where 


Q(n) =(- f| —2+ yr (y) dydu 


1 1 i 1 
0 


(30) 
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With n = 0, this expression agrees with Q(0), previously defined in equation (22). Be- 
cause of the exponential term, Q(m) — 0 as r—> ©. Further, we can write 


1 
Q(n) = 


where D is the operator of equation (3). Thus 
Q(1) = {Q(0) dr = Constant. (32) 


But Q(1) can vanish at infinity only if the constant itself is zero. Thus, proceeding step 


by step, we find, in general, that 
O(n) =Q. 


For L(1), the other integral vanishes identically, and we find that 
=3 (34) 
4 ’ 


which defines 5 in terms of jf}. We are to define a and $(y) in terms of 8, so as to satisfy 
for all 7 and ». The exponential factor in equation (30) is troublesome, but we can avoid 
the difficulty by adopting a set of weighted mean values through the layers. Multiply 
equation (30) by r*dz and integrate from 0 to ©. The resulting equations agree with that 
obtained from equation (30) when we set 7 = 0 and replace n by n + k + 1. 

Thus, if we can determine ¢(y) to satisfy the equation 


-f (a — bu) =0 


for all m, we shall have the correct solution. 
Let us write equation (21) in the form 


(r+ 


and redefine a and $(y) with respect to a multiplicative constant, so that 
1 
= dy. 
q (7) at f (y) dy 


As Hopf! has shown, g(r) is a slowly varying function of r. 
qg(~) =0.71044609=a, 


These results will provide a check on our later solutions. 

In accord with the convention of equation (37), we shall set 6 = 1 in equation (35) 
and carry out the integration with respect to uw. In the second bracket we divide yu" by 
(u — y). The fractional remainder just cancels the similar term in the first bracket, so 
that the apparent singularity at y = uw disappears completely. Thus: 


1+y a 1 


\ 


' Mathematical Problems of Radiative Equilibrium (‘Cambridge Tracts in Mathematics and Mathe- 
matical Physics,” No. 31 [Cambridge: At the University Press, 1934}). 


ca 
‘ 
5 
(31) 
7 
— 
— 
a 
i 
4 
(37) 
é 
(38) 
q (0) 7350269. (39) 
Vs 
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This is as far as we can carry the analysis without some explicit assumption about the 
function ¢(y). In the absence of other information, let us assume that ¢(y), whatever 
its explicit form may be, can be expanded in the series: 


o(y) = (41) 
0 


Negative values of j are clearly excluded because they would cause equation (40) to be 
infinite, when n = 0. When we set this equation in equation (40), the various terms of 
the coefficient of A; become: 

2 


1 
dy 
0 


(43) 
1 ( 1)*-1 


(n+ j even), 


1 


k=1 


1 


1 


1 n 1 n+j+1 1 j+1 1 
Write C,,; for the coefficient of A ;, so that 


a 1 


n+j+1 ( 1)*-1 n+j+l 1 j+1 n 1 
n+j+2 
for n + 7 odd, and 
n+j+1 n 1 n+j+1 1 J+1 1 


(-1) 


(48) 


for n + 7 even. Equations (46), one equation for each value of , are linear in the 4,’s. 
The C,,’s are numerical quantities, simply calculated. Some of the values appear in 


Table 1. 


6 
0 
(42) 
j 

= 
] t 
(44) 
ae 
a 
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ee 

AS 
(46) 
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(47) 
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As a first test of the efficacy of the procedure, the solution of four simultaneous equa- 
tions with four unknowns gave the values of a, Ao, A1, and As. This assumption amounts 


to breaking off the series after three terms. The last column of Table 1 was not used. 
At this point, we introduce a slight change of notation. With ¢(y) defined as in equa- 


tion (41), we obtain 


= 
0 


The quantity K,(7) is the exponential integral of the mth order, defined by 
e 
= ae. 


TABLE 1 
VALUES OF Cy; 


2 3 


— 1. 500000000 —0. 704568546 | —0. 458333333 —0.339407795 
—0. 204568546 — .125000000 — .089407794 — .069444444 
+0. 125000000 + .077258872 | + .055555555 + .043280149 
+0. 243925539 +0. 166666667 +0. 126613482 +0. 102083333 


The Mathematical Tables project, under the supervision of G. Placzek,” has prepared ex- 
tensive tables of the functions K,(7), for m = 0 to m = 20, and for r = 0 to r = 10. 
Satisfactory asymptotic expansions exist for other ranges. 

In the present study we break the series after a finite number of terms. Because of the 
fact that the behavior of K,(7) resembles that of Kn+:(7), the last term of the series ef- 
fectively represents the functional character of the neglected terms. 

Table 2 summarizes the results of several solutions, obtained according to various in- 
dependent assumptions. The solutions are expressed in the following form: 


0 


Although we made no attempt to fit the constants of Solution I, our series has repro- 
duced the theoretical values g(~ ) and g(0), of equations (38) and (39), with an accuracy 
better than one part in 10°. If, in Solution II, we adopt the exact value for g(@), we get 
q(0) accurate to one part in 10°. The remaining solutions, with one more term, are accu- 
rate to about two parts in a million. 


2 The Functions En(x) = - e-7“u-"du (Doc. No. MT-1) (Chalk River, Ontario: National Research 
Council of Canada, Division of Atomic Energy, [1946]). 


2 
(49) 
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There is an interesting and somewhat obscure point in the solution of the equations. 
The determinants resulting from the formal solution of the simultaneous equations (46) 
have a tendency to vanish. Formally, the tendency indicates an incipient indeterminacy. 
A linear dependence must exist either between the equations or between various parame- 
ters 1 ;. One consequence of this peculiar condition is the fact that the A,’s are not neces- 
sarily unique. Various linear combinations may give the same functional results to a high 
degree of accuracy. 

A shift of only one unit in the fifth decimal place of a produces a change of three units 
in the second decimal place of Ae, in Solution II. Solutions III and V were obtained by 
identical procedures. Both satisfy the initial conditions to nine decimal places. The dif- 
ference is that Solution V was carried out to twenty places of decimals, to remove the 
indeterminacy as far as possible. In Solution VI, the constant A 3 was arbitrarily adjusted 
to give an exact fit for g(0). The shift amounted to only one part in a million. 

Computations made from Solution IV with calculated values from Eddington,* Chan- 
drasekhar,* Placzek,® and LeCaine® appear in Table 3. Mark’s’ exact values, computed 
by machine methods, are given for comparison. 


TABLE 2 
CONSTANTS 


Solution a=q(@) Ao Ai Az A; | q(0) 


7104513 —0.2666602 | 0.3468140 | —0.1196287 | 
.71044609* | — .27086231 | .3764936 | — .1514577 |............| .5773447 
.7104460°* | — .28070408 | .54030420 | — .63919530 | 0.36208123 | .57734932 
71044930 — .27186451 .41860309 — 31239946 . 134339605 | . 56773815 
.71044609* | — . 28069496 .54016112 | — .63878155 | .36177922 | .57734932 


0.71€44609* | —0.28069496 | 0.54016112 | —0.63878155 | 0.36178017 | 0.57735027* 
| | | 


* Exact value assumed. 


From equations (8), (36), and (49) we easily derive the following formula for the 
emergent intensity: 
A, , Ae ( Ay 
(0, =(a+ - 41-53-53) 
(52)* 


ut+l 


A; 
+( 42+ — (Agu — Ayu?+ Am  — Agu*) log 


In the computation we shall adopt the coefficients as given by Solution IV. 


3 The Internal Constitution of the Stars (Cambridge: At the University Press, 1930), p. 333, eq. (231.2); 
see also S. Chandrasekhar, Ap. J., 99, 189, Table 1, 1944. 

4 Ap. J., 100, 82, Table 1, 1944. 

5 The Neutron Density near a Plane Surface (Doc. No. MT-16) (Chalk River, Ontario: National Re- 
search Council of Canada, Division of Atomic Energy [1947]). 

6 Phys. Rev., 72, 566, 1947. 7 Phys. Rev., 72, 563, 1947. 

® We have obtained two separate formulas for 7* and J~, where Chandrasekhar attempts to use one 
formula for the complete range of u (Ap. J., 100, 81, 1944). For example, he gets, in the second approxi- 
mation: 


. 0.116675 
2% 25 — 


Note that Chandrasekhar’s expression for 7(0, 4) has an infinite discontinuity atu = — 1/1.9723 and 
therefore does not satisfy the boundary condition (vanishing of the inward intensity) at this point. The 


= 
| 
| od 
| 
| 
NO 
4 
5 
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Table 4 gives in successive columns the values of /+(0, 4) as computed from equation 
(52), the exact values? (Hopf), Chandrasekhar’s fourth approximation,!° and the values 
as deduced from the H(u) function in Ambarzumian’s functional equation. Chandrasek- 
har" has computed H() by iteration from Ambarzumian’s functional equation 


dat 


H =1+ 34H 
+ (nu) =. 
TABLE 3 


q(r) 


Computed Correct Eddington 
(IV) Mark (2d Approx.) 


Chandrasekhar | 


(4th Approx.) | 


Placzek 


0.5773382 | 0.57735031 | 0.583 | 0.577 0.5774 
.6105623 |. | .5974 
.6278099 |, 626 .6139 
.6400963 645 .638 .6274 
6495643 648 6386 
.6634362 |. 662 .6557 
. 6686998 .670 .6621 
‘6731871 | 67: 673 672 6676 
.6803844 680. | 680 6761 
6858845 683 086 6823 
.6901802 . 686 .6870 
6935930 689 695 
.6963401 696: 692 698 6933 
6985750 | .6985 694 6954 
. 7093432 . 703 .7058 
. 7096034 : . 7064 
. 7097999 .7098 0.711 70985 

0.7104494 0. 7104461 Correct value 0.7069 | 0.7104457 

assumed 


0. 
0. 
0.70.. 
0.80. 
0.90. . 
\ 
&. 
2.0. 
22... 


and we have computed /(0, «) from the relation: 


1(0,n) (u). (54) 


third approximation for /(0, u) has two such infinities. Further, these points of infinite discontinuity mul- 
tiply without limit in the process of infinite approximation. In view of the existence of these infinities, the 
limiting function for 7(0, 4) obtained by Chandrasekhar’s method is not the exact Hopf function but dif- 
fers from it at an enumerably infinite set of points (of measure zero). In other words, the two functions 
are equal almost everywhere. It is precisely because of this property that the differences between the two 
functions, though real and not to be ignored from a purely mathematical viewpoint, are nevertheless physi- 
cally insignificant. 


*S. Chandrasekhar, Ap. J., 99, 189, Table 2, 1944. 
1 [bid., 100, 82, Table 2, 1944. 
 Jbid., 106, 143, Table 1, 1947. 
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Mark” has obtained an asymptotic expansion near the boundary, from which he has 
constructed Table 5 for g(r). 

The high accuracy of the approximations and the general simplicity of the results ob- 
tained by the present method are gratifying. The foregoing procedure, developed for 


TABLE 4 


(0, 


Ambarzumian’s 
Functional 
Equation 


Chandrasekhar 


Computed IV Exact (Hopf) (4th Approx.) 


0. 4330036 | 0.43301 0.43301 
0. 5400318 0.54011 : 0.54014 
0.6279790 0.62802 : 0.62799 
0.7112144 | 0.71123 : 0.71122 
0.7920933 0.79210 0.79211 
0.8715594 0.87156 0.87156 
0.9500923 0.95009 : 0.95007 
1.0279502 1.02796 1.02797 
1. 1053614 1.10536 1.10534 
1. 1823876 1.18238 1.18238 
1. 2591280 1.25912 1.25911 


Mark's Asymp- 
totic Expansion Correct 


0.5773 0.5773 
5882 | 
5054 
6287 
0.6522 


strict radiative equilibrium, can be extended to the case of nongray atmospheres. We 
have to replace equation (15) by a more general equation, wherein m is the root of the 
expression: 


1 
= 
8 +—tanh-'m=0. 


We shall return to these more general problems in a later paper, if further study indicates 
that the method is useful. 


12 Op. cit., p. 562. 


‘ 
1 | | 
| 
r | : 
Ae 
| Solution IV 

6106 
a 


TRANSFER OF RADIATION 11 


In a recent paper, Kourganoff'* has used the expansion (51) of our paper and deter- 
mined the constants from the conservation of the net flux. He has fitted his solution for 
sixteen values of the optical depth 7 by the method of least squares. Kourganoff has indi- 
cated in a very recent letter that he has reached an accuracy of about two parts in a 
million in a six-term expansion. 


We would like to thank Mr. David Murcray, Mr. Willard Shrode, and Mr. John G. 
Wolbach for having helped us with the extensive numerical work in this paper. 


'3CR, November 15, 1948, p. 1020; Contr. T.A.P., Ser. A, No. 34. 
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SUPERTHERMIC PHENOMENA IN STELLAR ATMOSPHERES 
V. ON EMISSION LINES AT HIGH KINETIC TEMPERATURE 


RicHARD N. THomas* 
Institute for Advanced Study 
Received March 17, 1949 


ABSTRACT 


The treatment of the radiation field of a high-kinetic-temperature chromosphere is continued. The 
conditions for emission lines to be formed and for an effective photosphere to be found in the chromo- 
sphere are investigated from a purely formal standpoint. For a quantitative investigation of the behavior 
of the Balmer series, particularly with regard to the transient appearance of emission lines,a detailed 
investigation of inelastic collision cross-sections seems required. 


If one admits the quite general existence of an outer stellar atmosphere that takes the 
form of an extended high-kinetic-temperature chromosphere-corona, it appears possible 
to resolve several of the difficulties connected with the anomalous presence of emission 
lines; for the existence of such a high kinetic temperature implies a source of energy other 
than the photospheric radiation field. Hence we find the possibility that, in the chromo- 
sphere, energy from this other source eventually becomes radiation in a limited spectral 
region. In short, the atmosphere departs from the condition of radiative equilibrium. 


The principal problem may be stated as follows: What is the contribution of the 
chromosphere to the general radiation field that the customary methods applied to the 
classical atmosphere predict? The question may be broken into two parts: (1) Does the 
chromosphere add to, or subtract from, the intensity calculated at the photospheric re- 
gion? (2) When does the effective photosphere move into the chromosphere? 

The terms “‘photosphere,” etc., lose much of their meaning in the above context; so 
hereafter we refer to those regions of the atmosphere where 07,/dr < 0 as the lower at- 
mosphere or classical photosphere, and the regions where 07,/dr > 0 as the chromo- 
sphere. With the term “effective photosphere” we shall attempt to retain the customary 
connotation of the term ‘“‘photosphere,” viz., the region of optical depth 1. Thus, when 
the effective photosphere in any spectra] region lies in the chromosphere, we conclude 
that the classical-atmosphere description breaks down completely. For the majority of 
the stars, however, it appears that the energy distribution in at least the visible region of 
the continuum refers to the classical photosphere. The phenomena then investigated are 
the spectral lines and perhaps the far ultraviolet region of the spectrum. 

In this preliminary discussion we consider purely formal answers to questions 1 and 2. 
That is, the chromosphere as defined above is completely characterized only when the 
population of each spectroscopic state is explicitly given; for the situation is, in general, 
far removed from thermodynamic equilibrium. We assume this description to be given 
by the 8, (introduced in Paper IT)! and express the formal solutions to questions 1 and 2 
in terms of these b,. We do not, however, attempt to calculate the b, but assume that 
their values are given. Solutions! for the solar chromospheric 6, will be used for numerical 


examples only. 
* Frank B. Jewett Fellow for 1948-49. 
1 Paper II, 4p. J., 108, 142, 1948; Paper III, 4p. J., 109, 480, 1949. 
12 
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SUPERTHERMIC PHENOMENA 13 


I, THE DIFFERENTIAL CHROMOSPHERIC CONTRIBUTION AND QUESTION 1 


Elsewhere it has been shown? to be extremely likely that the 6, exceed unity, in this 
case of high kinetic temperature. We have the defining equation for 6,: 


h? X 


where JN, is the population of level » and the other symbols are conventional, and x, is 
measured relative to the ground state of the ion whose density is V;. With the exception 
of metastable levels, we expect 5,’ > b, for n > n’.3 

Menzel has earlier remarked‘ on the condition for appearance of emission or absorp- 
tion lines when 6, # 1. He states the condition as emission or absorption according as 
b, exceeds or is less than ,’. Fortunately for the hope that the high 7, will help clear up 
the emission-line problems, Menzel’s criterion is not generally valid. We may reach 
more complete criteria by considering the expression for the difference between emission 
and absorption in a line. By writing down the Einstein transition conditions and using 
equation (1), we obtain, for the net emission, 


bn 
AE = Const. 1-—S 


where we write: S = ratio of intensity at v,’, to that of a black body at T,, and T, 
and 7, are the kinetic and radiation temperatures, respectively. Thus we obtain the 


conditions: 


Te 


< 
AES 0 according as },, > b,, —1 +S" 


Thus the modification of Menzel’s conditions is quite large when 7, differs much from T, 
even if we were to consider AE alone as determining the emission or absorption. Equa- 
tion (3), however, only provides an answer to question 1 above, relating to the positive 
or negative value of AE. We note that AE < 0 is both a necessary and a sufficient con- 
dition for an absorption line to become still deeper in the chromosphere. Thus, even 
though T, > 7,, a high enough value of 6,’ will provide blacker lines than the classical 
theory would predict. The result has a clear connotation for those arguments which ex- 
plain the absorption by contrasting the optical depth within and slightly adjacent to the 
line. For large enough 6,’, the optical depth at the line center would correspond to a 
temperature far greater than that of the continuum, and we should then expect emission 
by the usual arguments. This point is a useful one to illustrate the effect which departures 
from thermodynamic equilibrium can have on equilibrium arguments. 


Il. THE INTEGRATED CHROMOSPHERIC EMISSION AND QUESTION 2 


We see that AE > 0 is a necessary, but not sufficient, condition for the formation of 
an emission line. For an emission line, the integrated net emission must fill up the original 
absorption line. Following Milne,’ we may write down the equation of transfer in terms 


2 4p. J., 109, 500, 1949. 

3 We adopt the convention n’”’ > n > n’. 

4 Pop. Astr., 47, 17, 1939. 

6 Handbuch der Astrophysik (Berlin: Julius Springer, 1930), chap. ii. 
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of the detailed processes.* We introduce the b,, however, from equation (1), and obtain 


2hv® by’ ho/kT, 
=-T/, e 1 (4) 
where 

dr= 1 ds. (5) 


n 


Taking ds as measured positively outward, 7 = 0 at the bottom of the chromosphere, 
and 7» as the optical thickness of the top of the atmosphere along the direction s, there 


results: 
= I,(0) e-~+—; e wf ev -1] dr. (6) 


Any significant value for 7» will indicate that the chromosphere contributes to the 
structure of the absorption line. The perturbation by the chromosphere may be stated 
as the following: Suppose we are given some height in the classical photosphere above 
which lie Nw atoms of hydrogen. For rather low photospheric temperatures 7,, we will 
have essentially (mz, is the density of atoms in the ground state): 


for the number of atoms in state n’ above the photosphere, and 8, is the isothermal 
logarithmic €ensity gradient at 7,. On the other hand, the existence of a chromosphere 
at high kinetic temperature beginning at the same level will give, essentially (my, is the 
density of ionized hydrogen), 


since effectively all the hydrogen will be ionized. Using the condition 
ny 


we have 
N,! 


Nn’ Ti? 


Thus for high enough 8,’ a very small 7) computed on a classical basis may become a ro 
greater than unity on the high-temperature-chromosphere basis. For example, substitut- 
ing the solar chromospheric values 7, = 3.5-104 and N, = 2-10", we obtain 


Bre 0.5-10-*b,. (10a) 
No, 
When we solve for the solar chromosphere,' neglecting collisional excitation from level 2, 
we obtain by as some 10‘. Thus the optical depth at the line center has increased by a fac- 
tor of 50. For the solar case it appears that the centers of the early Balmer lines are 


formed rather high in the chromosphere. 


® We note, by contrast with Milne’s further development of this equation, that we are very much 
concerned with cyclic processes and therefore cannot solve for the ratio b,/b,- independently of the values 
of the other d,. In the following, we imply always the possibility of solution for the b, by a development 
similar to that carried out for the hydrogen atmosphere of Papers II-III. 
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Returning to equation (6), we may further illustrate the situation by considering an 
isothermal atmosphere with constant b,. We obtain 


For large ro, the condition for a chromospheric contribution sufficient to make J,(ro) 2 
is 


1) > 1) (0). (12) 


Condition (12) is just that given by equation (3). Thus, under this restricted situation 
equation (3) becomes necessary for emission and sufficient for nondepletion of the classi- 
cal photospheric field by the chromosphere. For emission, we require I,(ro) > I,(T,). 
The condition here is simply (12) with S(O) = 1, if to > 1. 

For the nonisothermal case, the same conclusions hold if the integrand of equation (6) 
is, apart from the e’ factor, a sufficiently slowly varying function of 7, or in any event if 
it is an increasing function of r. 

Finally, we note that it has been shown in an earlier paper! that the condition 

arty: | >1 
ds : 


leads for the ratio of chromospheric field to classical photospheric field: 


Ratio = "= an (0), (14) 


when we neglect the induced emission. The condition S,’, 2 1 insures that, somewhere 
in the chromosphere, absorption lines have disappeared. 

The same sort of criterion may be applied to the continuum, where we now replace b, 
by unity when we consider the free-bound continuum. Condition (14) follows directly 
from equation (6) if we lay down the condition equivalent to (13), namely, that in an 
optical thickness defined by the value (d/nN./ds)™ for ds, the integrand of equation (6) 
does not materially change, apart from the e’ factor. 

Again by way of illustration we may quote the results! for the solar chromosphere. 
The essential region of interest is in the Lyman lines; for the results are not so sensitive 
to the neglect of collisional excitation from level n = 2. We obtain a very large ro when 
we choose as the base of the chromosphere that region for which Wildt obtains’ V, = 
2-10" and T, = 3.5-10*. We obtain a value S;, >> 1 in equation (14). Thus we con- 
clude that the Lyman region in the sun has an energy output far exceeding that of a black 
body at 6000°. The exact numerical value of S;, depends on the be value, if b: is suffi- 
ciently large. In any event, however, the Lyman contribution will not be less than 20d, 
for Si2, at the chromospheric level in question. 


III. THE MARGINAL APPEARANCE OF EMISSION LINES 


The considerations on the Balmer lines are of somewhat greater interest than those on 
the Lyman lines because of their more accessible location in the spectrum. In particular, 
when one considers emission lines, those occurrences* where the early members of the 


TAp. J., 105, 36, 1947. §Struve, A.J., 53, 108, 1948. 
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series are in emission and the later ones in absorption become of interest. The necessary 
conditions for such a situation are, from equation (3) [neglecting S and 1]: 


—he/kT WkTe —hW/kT (15) 
Balmer limit . 
Thus the essential question—for given 7, and T,—is the 62/6, ratio. 


From the work on the solar chromosphere, where collisional excitation from level 2 
has been neglected, one finds for the above relation on H,: 


bs Te = 1 
bs ew — 1 +o], 


where 
O<n, p<l 


and e is of the form: 


€1 


== 0 1 
€ <a <i, 


and 5S» is S evaluated at H, and averaged over all directions. Thus S2 is somewhat less 
than the observed residual intensity of H.. We see that in this case we expect no emis- 
sion in the Balmer series 2s a consequence of the presence of a high-kinetic-temperature 
chromosphere. 

If, however, we add to the equations of condition on b, a term involving collisional 
excitation from level 2, we must divide the expression (16) by a term of the form 


a7) 


where C' is some constant depending on the cross-section for collisional excitation from 
level 2. For sufficiently great V, and T,, we then see the chance of bringing the earlier 
chromospheric lines into emission. 

Such a situation as pictured may well provide the possibility of treatment of the ap- 
pearance of emission lines in RR Lyrae; for, at the stage of the velocity-curve where the 
emission lines occur, whether one speaks of shock waves or waves of finite amplitude, one 
must admit the presence of a 7, an order of magnitude or more greater than the usually 
accepted T,. Also, V, must rise. One then has the requisite conditions for the appearance 
of emission lines. A quantitative treatment of the situation should provide the series 
member where the emission changes to absorption, if one is given JT, and N,. Or, not 
knowing these parameters, one might investigate the possibility of inferring them from 
the observations on series members where emission changes to absorption. As a prelimi- 
nary step, however, one requires the evaluation of collisional cross-sections in the Balmer 
lines. 

In summary, we find, upon inquiry into the purely formal expression for the condition 
that emission lines appear in a star possessing a high-kinetic-temperature atmosphere, 
that a knowledge of the inelastic collision cross-section is essential. So far as we can go in 
purely formal development, the solutions indicate that a typical case of the appearance 
of emission lines—RR Lyrae—is compatible with the formal relations. It would then 
appear that the calculation of thermal collision cross-sections for H is the next step in 
this study. 


I am greatly indebted to Dr. Lyman Spitzer, Jr., for reading and discussing the paper 
during its preparation. 
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SELF-CONSISTENT FIELD WITH EXCHANGE FOR Ca xu1* 
W. H. Davist 


Brown University, Providence, R.1. 
Received March 14, 1949 


ABSTRACT 


Since Ca XItI appears to be responsible for one of the coronal lines, it seemed desirable to find the "3 
Fock wave functions for the 3P ground state of this ion. These wave functions are tabulated. Two rough i 
checks have been made with experiment. One check is the comparison of the Hartree energy parameters 
with observed energy difference, and the other is the comparison of the calculated and observed energy 
differences corresponding to the 3P;—3P, transition responsible for the coronal line. 


INTRODUCTION 


In 1945 Edlén! identified certain spectral lines in the sun’s corona with transitions in 
highly ionized atoms. One of these lines was attributed by him to the transition 3P;—3P2 : 
in the configuration (1s)?(2s)*(2P)* of Ca x1m (a calcium atom ionized twelve times). ; 
This line is forbidden by the ordinary selection rules and, according to Edlén, occurs i 
only in sources of very low density. It was decided to make a theoretical investigation, : 
using the Fock method, of the 3P ground state of Ca xm, and to compare theoretical 
and experimental energy differences. 


METHOD OF SOLUTION 


The Fock method was chosen for the investigation of Ca x11, since it provides the 
most accurate wave functions for ions or atoms with as many electrons as Ca xu. How- 
ever, before using the Fock method proper, it is necessary to find the wave functions 
by the Hartree method, since these provide a good first approximation to the Fock i 
solutions. The fundamental difference between the Hartree and Fock methods is that 
the former uses a resultant wave function which is a determinant of the 1-electron wave if 
functions, The resultant Fock wave function, therefore, satisfies the Pauli exclusion v1 
principle, but the resultant Hartree wave function does not. i 

The Fock and Hartree methods have been applied to many atoms. The details of the ' 
method can therefore be found in many places.? Hence we shall content ourselves with f 
a few general statements. We assume that the 1-electron wave functions are products E 
of two functions. The first gives the dependence on the angular position found in the hy- ‘ ; 
drogen atoms, and the second gives the as yet unknown dependence P,,:(r) on the dis- : 
tance from the nucleus. The differential equations which the P,,,(r) functions must satisfy 
are then determined from the variation principle of quantum mechanics. For the special 
case of the 3P ground state of Ca xu we get the following Fock equations: 


1 


* Part of a dissertation presented to the Graduate School of Brown University in partial fulfilment 
of the requirements for the degree of Doctor of Philosophy. 


t Now at the University of Buffalo, Buffalo, New York. 
1M.N., 105, 323, 1945. 
2D. R. Hartree, Proc. Cambridge Phil. Soc., 24, 89, 1926. 
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18 W. H. DAVIS 


where 
= [Wo (1s, 1s) + 2Wo(2s, 2s) +4Wo(2p, 2p), 


Ti, (17) = [Wo (1s, 2s) — Pos (7) +[27Wi (1s, 2p)] Pop (1) , 

= [2Wo (1s, 1s) +Wo (2s, 2s) +4Wo (2p, 2p)], 

T2,(r) = [Wo (1s, 2s) —e1s, 25] Pis (7) + [2 (2s, 2p)] Pop (1), 

2Vap(1r) [2Wo (1s, 1s) + 2Wo (2s, 2s) +3Wo (2p, 2p) — 3 (2p, 2p)], 


(1) rW, (Is, 2p) Pis + rW, (2s, 2p) Po» (1), 


Wx (nl, n'l') = 2f Lf P.,(p,) Pia) 


The Hartree equaticns for the ground state of Ca x11 differ from the above Fock equa- 
tions only by the omission of the 7,:(7) and the W2(2p, 2p) terms. 

For the actual integration of the Fock and Hartree equations, analytical methods are 
impracticable. Experience of others has suggested that the Milne method of numerical 
integration is the most useful when an electrical calculator is available. 

The general technique of solving the Fock and Hartree equations is one of successive 
approximations. We make successive estimates of the V,.(r) and én, n’’ until the wave 
functions satisfy the boundary conditions as well as the differential equations. The 
T,.(r) terms in the Fock equations make normalization by a constant multiplier impos- 
sible. This necessitates the adjustment of one additional constant and thus renders the 
solutions of the Fock equations even more lengthy than the solutions of the Hartree 
equation. 


RESULTS 


The Hartree wave functions have been obtained for the (1s)?(2s)*(2p)4 and the 
(1s)?(2s)(2p)* states of Ca x11. The Fock wave functions have been obtained for the 3P 
ground state of Ca xm. The Fock wave functions and the corresponding ¢ energy param- 
eters are listed in Table 1. (The e’s are not the true energies but only approximations to 
the true energies.) The solutions are accurate to about four decimal places or better for 
small r and about three decimal places for large r. The € parameters are accurate to 
two or three decimal places. The ratios a¥,/a%, below have the following significance. 


The a”, are the quantities lim P,./r*! for the Hartree case, and the a%, are the corre- 


sponding quantities for the Fock case. The ratios of these quantities may be useful in 
solving the Fock equations for other highly ionized atoms: 


F F F 
Gis 


1.0250. 


We know that an excited state should have less energy than the ground state. This 
condition is satisfied by the € parameters. The greater energy would also correspond to 
a greater concentration of charge. This condition is predicted by the wave function also, 
since the ground-state wave functions are shifted toward the origin as compared to the 
excited state of wave functions. 
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The difference between the Fock solution and the Hartree ground-state solution is 
rather small, but not negligible, for this tightly bound ion. As usual, the Fock solution 
predicts a more compact atom and, as a result, also gives larger e parameters. This is a 
step in the right direction, since it is well known that the Hartree solutions in general 
give an atom which is not compact enough. 


TABLE 1 
FOCK WAVE FUNCTIONS FOR 3P STATE OF Ca XIII 


Prs(r) Pop(r) 


S 


+0.2708 
+0.0858 
—0. 1020 
—0. 2875 
—0.4672 
—0.6384 
—0.7992 
—0.9483 
—1, 
—1. 
—1. 
—1. 


SSSSSSSSSSS 


—1. 
—1.7 
—1. 
—1. 
—1. 
—1.73 
—1. 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


COMPARISON WITH EXPERIMENT 


We have two rather rough checks with experiment. The first check is a comparison 
of the energy of transition from the excited (1s)?(2s)(2p)® state to the ground state, 
(1s)*(2s)*(2p)*, as deduced from the Hartree ¢ parameters, with that deduced from the 
corresponding spectral lines (this is not the coronal line) (Table 2). The agreement is 
only within 20 per cent. This large discrepancy is not surprising, since, according to the 
Hartree theory, the € parameters give only approximate values for the energies. To get 
a more accurate test, we should have to evaluate the energy integrals.* 


3D. R. Hartree and M. M. Black, Proc. R. Soc., London, A, 139, 311, 1933. 


' Pur) | Pir) Pir 
0 0 | 0.046 | —1.6061 | 1.3950 
005......| 0.7920 | 40.2399 | 0.0064 || 0.034 —1.5209 | 1.3151 = 
010......| 1.4336 | +0.4328 | 0.0243 0.025 —1.4481 | 1.2344 
O15......| 1.9465 | +0.5840 | 0 .....| 0.018 —1.3630 | 1.1542 
020......| 2.3496 | +0.6984 | 0 0.013 —1.2766 | 1.0754 
025......| 2.6594 | +0.7806 | 0 0.010 —1.1904 | 0.9987 
030......| 2.8001 | +0.834 | 0 0.007 —1.1055 | 0.9247 
ae 035......| 3.0840 | 40.8635 | 0 0.005 —1.0230 | 0.8538 
040......| 3.1618 | 40.8713 | 0 
3.2228 | 40.8599 | 0 —0.8675 | 0.7226 
050......| 3.2448 | +0.8328 0 01002 | -0:7275 | 0.6063 
055....../ 3.2348 | +0.7919 | 0 0.001 | —0.6044 | 0.505 
060......, 3.1986 | +0.7394 | 0 0.001 | 0.498 | 0.418 
065......| 3.1412 | +0.6770 | 0 
075......, 2.9796 +0.5295 | 0 —0. 268 0.229 
msO......) 2.8821 | +0.4471 | 0 
2.4329 1 0.109 | 0.097 
1.9845 | 1 70.068 | 0.061 
0.7095 | 70.008. | 0.009 
7588 || 0.003 | 0.005 
7669 0.002 | 0.003 
7527 | 0.001 | 0.002 
5463 | | ep =54.24 
A127 | 
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We can make a second check by comparing the energy difference corresponding to 
the 3P;—3Pz transition in the ground state of Ca x11 as calculated by the perturbation 
theory! with that determined from the observed coronal line. In the theoretical calcula- 
tions we consider the spin-orbit interaction (Russel Saunders coupling) as the per- 
turbing term. When we perform this perturbation, making use of the properties of the 
operators involved, we get the following relation: 


h? 1 OV» 
AEs, =2 16m2c? I Pap (7) E Or 


where V2,(r) is the potential in which the 2p electron is assumed to move and AE3p,_ 3p, is 
the energy difference desired. When we make use of the Fock P,, functions to evaluate 
the above equation, we get 22,240 cm~'. The experimental value computed from the 
coronal line is 24,464 cm. The discrepancy of about 10 per cent is not surprising, con- 
sidering the assumptions that we have used. Probably the greatest part of this error lies 
in the neglect of the magnetic interaction of different electrons. However, this compari- 
son does give a rough check on the wave functions. 


TABLE 2 


EXPERIMENT* From ¢ PARAMETERS 


Experimental 
Identification Value | Identification 
(Cm!) 


618,246 
503.782 


*B. Edlén, M.N., 105, 323, 1945. 
CONCLUSION 


The Fock wave functions have been found for the 3P ground state of Ca x1.° Two 
rather rough checks have been made with experiment, which show that the wave func- 
tions are at least reasonably close to the true ones. These Fock wave functions could be 
used to evaluate other quantities of astrophysical interest. For example, the author in- 
tends to calculate the transition probability for the ‘‘forbidden” coronal line attributed 
to Ca x11. In addition, the Fock wave functions could be useful in making good initial 
approximations to the solutions for other highly stripped atoms responsible for coronal 
lines. It is intended at Brown University to carry through the Fock solution for Ca x11. 


The author wishes to express his sincere thanks to Professor A. O. Williams, Jr., who 
suggested the subject of this thesis and who supervised the general development, as well 
as to Professor R. T. Beyer, who made the initial estimates for the potential field of 
Ca 


4E. U. Condon and G. H. Shortley, The Theory of Atomic Spectra (Cambridge: At the University 
Press, 1935). 


5 On request the author will be pleased to send tables of potential functions or any other data required. 
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THE MILKY WAY AT THE JUNCTION OF GEMINI, 
MONOCEROS, AND ORION 


Bart J. BoK, MARGARET OLMSTED, AND Betty D. BOUTELLE 


Harvard College Observatory 
Received April 15, 1949 


ABSTRACT 


The present study is based largely on photographically determined color indices and magnitudes of 
stars in the Henry Draper Extension tor the section between galactic longitudes 160° and 173°, latitudes 
—3° and +9°. 

The structural features of the region are surveyed, and an attempt is made to determine what is the 
greatest well-determined distance for an object in this section of the Milky Way. The conclusion is 
reached that no known clusters or nebulae, O stars, c stars, or B stars are found at distances greater than 
3000 parsecs and that the only known objects that are possibly at greater distances from the sun are 
three cepheid variables, one of which may be as distant as 5000 parsecs. It seems unlikely that, in the 
anticenter direction, the Milky Way system stretches more than 4000-5000 parsecs beyond our sun. 

The space reddening is highly variable over the section. For Field 1, at / = 171°, b = 0°, we find 
relative freedom from overlying obscuration, with a value of A,, S 0.5 mag/kpc to a limit of 3000 
parsecs indicated. 

For Field 2, at / = 166°, 6 = 0°, we find, however, A,g = 1.2-1.5 mag/kpc to a limit of 1500 parsecs, 
in agreement with an independent spectrum-color survey (Ap. J., 109, 139, 1949) by McCuskey. Fields 
3 and 4, at galactic latitudes +3° and —2°, adjacent to Field 2, yield space reddening comparable to that 
found for Field 2. Fields 5, 6, and 7, with / = 164°-170° and 6 = +7”, are all relatively free from ob- 
scuration, with A,,g = 0.5 mag/kpc to a limit of 2000 parsecs representing a fair average result. 

Incomplete color data for three fields with galactic latitudes —3° to —7° give evidence for the presence 
of strong near-by overlying absorption, presumably caused by extensions of the Taurus-Orion complex 
of dark nebulae. One field, No. 11, is reserved for subsequent and more detailed study; it is located near 
1 = 171°, b = +3°, and falls in the region of clusters and bright nebulosity south of Wolf’s great dark 
nebula. The space reddening appears to be quite irregular. 

In the preliminary analysis of the star densities, considerable stress is laid upon a comparison with the f 
results obtained previously for a region at / = 185°, b = 0°. Negative density gradients, with an average 
density 0.5 © ata distance of 1000 parsecs, are found for the B8-AO stars; the Al—A7 and FO-F§8 stars 
yield average densities equal to 0.5 © at distances of 500 parsecs or less. For the late-type giants the 
average space densities are fairly constant over the first 500 parsecs from the sun. : 


cor 


In 1945 Bok and Mrs. Rendall-Arons'! published the results of an investigation on the i 
stellar distribution for a section of the Milky Way in Monoceros. The section under in- ; 
vestigation was between galactic longitudes 176° and 186°. The present investigation 
deals with the section between galactic longitudes 160° and 173°, directly north of the ? 
one previously investigated. The basic materials for this new study are spectral types ; 
from the Henry Draper Extension® and blue-red color indices determined from plates 
taken with various cameras at the Oak Ridge station. 

The Milky Way at the junction of the three constellations, Gemini, Monoceros, and 
é Orion, has the distinctive smooth appearance of the anticenter regions. The section un- 
der investigation lies between a = 6-7 hours; 6 = +5° to +20°; it is shown nicely on 
Chart 35 of the Ross-Calvert Adlas of the Milky Way. The southern boundary of the sec- 
tion is roughly the line between Betelgeuse and Procyon, whereas y Geminorum is the 
northern apex. To the west of the section, beginning at galactic latitude —3°, are the 
complex obscured regions associated with the Orion and Taurus nebulosities. Since we are 
—in accordance with the plan outlined in the earlier paper—primarily interested in re- 
gions free from obvious obscuration, we have omitted regions south of galactic latitude 


1 Ap. J., 101, 280, 1945; Harvard Reprints, No. 273. 
2 Harvard Ann., Vol. 100, No. 5, 1931. 
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= —3°, whereas to the north of the galactic circle we have gone to b = +9°. We had 
originally intended to study colors and spectra for eleven centers, but a preliminary in- 
vestigation led to the elimination of four centers affected by rather irregular local ob- 
scuration. The seven remaining centers border the conspicuous obscuration to the north 
of S Monocerotis, which has been studied in a classical paper by Wolf.’ The obscured re- 
gion has also been studied by Andrews,‘ who finds that the dark nebula at a distance of 
about 800 parsecs produces a total absorption of 1.5—2.0 mag. Before we proceed to de- 
scribe the results of our own analysis, we shall describe briefly what, on the basis of the 
available published material, may be learned about the galactic structure for this section 
of the Milky Way. 

I. STRUCTURAL FEATURES 


The distribution of the dark nebulosity over the region has been described very well 
in Wolf’s paper. The three principal centers of bright nebulosity are near 12 Mon (a 
cluster, NGC 2244 and bright nebulosity, NGC 2237), near S Mon (= 15 Mon; NGC 
2264), near R Mon, and near BD+10°1159 (NGC 2245 and 2247). Struve and Elvey® 
and Struve, Elvey, and Linke® list, altogether, nine patches of nebulosity in this section 
of the Milky Way. The emission nebulosity near S Mon (Sp. O7s) has been found by 
Struve and Elvey and by Greenstein and Henyey’ to have strong [O 11] 3727 and rela- 
tively weak Ha; several near-by patches have the same characteristics. NGC 2237 near 
12 Mon shows strong Ha, HB, Hy, Hé, N: + Ne, and [O 1] 3727. Finally, NGC 2245, 
with an associated star of spectral class B8, shows an absorption spectrum. 

The T Tauri variable, R Mon, and the associated variable nebula discovered by Hub- 
ble have recently been studied by Greenstein,® who has also given references to earlier 
studies of the star and nebula. Greenstein places the star and nebula at a distance of 
400-500 parsecs, at about the distance of the near-by star S Mon and its associated bright 
and dark nebulosity. 

- The region near 12 Mon is especially interesting. Five of the eleven Henry Draper 
O stars between a = 5°40™-7500™ and 6 = +3° to +20° belong to NGC 2244, and three 
of the remaining six O stars are within 2° of the cluster. The distance and total absorp- 
tion of NGC 2244 can be reliably estimated from Zug’s® spectra and colors and from the 
Stebbins, Huffer, and Whitford colors.!° The distance to NGC 2244 is of the order of 
1300-1400 parsecs, and the total absorption between our sun and the cluster is 1.3 mag. 
This leads to a value for the photographic absorption per kiloparsec of 


Avg = 1.0 mag/kpc . 


The color excesses of three O stars directly north of NGC 2244 indicate an average photo- 
graphic absorption of the order of 


Ay, = 0.5 mag/kpc . 


These absorptions are confirmed by the colors of B stars. Since NGC 2244 lies at the 
northern border of a region of obscuration, it is not surprising that the value of A), for 
centers a few degrees to the north of the cluster is smaller than that for the cluster itself. 


3 Festchrift fiir Hugo v. Seeliger (Berlin: Julius Springer, 1924), p. 313. 

4 Harvard dissertation, unpublished; for summary see Pub. A.A.S., 7, 211, 1933. 
5 Ap. J., 89, 522, 1939. 

6 Ap. J., 90, 301, 1939. 7 Ap. J., 87, 80, 1938. 


8 Centennial Symposia (‘‘Harvard Observatory Monographs,” No. 7 [Cambridge, Mass.: The Observa- 
tory, 1948]), p. 25; and Ap. J., 107, 377, 1948. 


® Lick Obs. Bull., No. 454, 1935. 
10 4p. J., 91, 20, 1940; Mt. W. Contr., No. 621. 
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A small value of A,, ($0.5 mag/kpc) is also found for the cluster NGC 2264. Trum- 
pler" places this cluster at a distance of 450 parsecs. A slightly larger value might be 
preferable, but the cluster is not likely to be at a distance in excess of 600 parsecs. To this 
distance the reddening is slight and the absorption small though patchy. 

Trumpler" lists twelve galactic clusters in the region under investigation. With the 
exception of three, he places these clusters at distances under 3000 parsecs. The most re- 
mote cluster, NGC 2141, has an estimated distance of 4860 parsecs. We doubt that this 
distance can be accepted without reservation. The cluster lies within the limits of one of 
the four centers which were eliminated from our investigation because of excessive over- 
lying obscuration. Color-index measurements by Miss Olmsted indicate that this over- 
lying absorption amounts to 2 mag. at 1000 parsecs and that Trumpler’s distance must 
be reduced considerably. Trumpler describes the cluster as a ‘“‘well resolved thin but rich 
cluster of exclusively faint stars (15-18 mag).” 

The other two distant clusters in Trumpler’s list are NGC 2259 and NGC 2269, both 
with estimated diameters of 3’-4’ and distances of 3900 parsecs. The first of these is in a 
clear patch bordered by heavy obscuration north and west of S Mon, whereas the second 
is near obscuration directly east of 12 Mon. Again the distance estimates by Trumpler 
are most likely to be too large. Cuffey’® has shown this to be the case for several small 
clusters for which Trumpler and others had estimated the distance from angular diam- 
eter and appearance alone. A case in point is NGC 2194, for which Cuffey'’ finds a dis- 
tance of 1200 parsecs from the color-magnitude array, whereas Trumpler’s estimate was 
2520 parsecs. Cuffey finds, incidentally, for NGC 2194 a space reddening equivalent to a 
total photographic absorption of 2.7 mag. at the distance of the cluster. 

It is interesting to note that no cluster in the section is known with certainty to be at 
a distance from our sun in excess of 3000 parsecs. Since the section of the Milky Way un- 
der investigation (J = 160°-173°) is one of the least obscured in the general direction of 
the anticenter, considerable interest attaches to the greatest distance to which galactic 
objects can definitely be traced in this section. Among the known objects that are pre- 
sumably most remote and for which fair distance estimates can be made are the O stars, 
the B stars, the c stars, and the cepheid variables. 

We have already commented on the distribution of the eleven known O stars in this 
section of the Milky Way. Nine of these O stars are in or near 12 Mon, and, with the 
possible exception of HD 45166, these are certainly not farther than 1800 parsecs from 
the sun. S Mon is another O star, well within 1000 parsecs, and the single remaining O 
star (m = 7.1) isalso relatively near by (650 parsecs). 

The B stars are also scarce among the fainter stars in the section / = 160°-173°. The 
Henry Draper Extension for an area of 340 square degrees yields only nine new BO stars 
and a single Oe5 star. Bok and Miss Wright" have already drawn attention to the pau- 
city of faint early B stars for a field in Monoceros (/ = 180°). In Table 1 we have listed 
side by side the distribution of the BO—-B2 and B3-BS stars for an area of 67 square de- 
grees in Monoceros (/ = 180°, 6 = 0°), for an area of 65 square degrees at / = 165°, 
b = 0°, and for an area of 60 square degrees at / = 165°, b = +5°. It was noted by Bok 
and Miss Wright that the numbers for the Monoceros field represented a rapid thinning- 
out of the B stars with increasing distance from the sun. This is even more marked at 
b= +S°. 

The total number of known ¢ stars is eight, whereas one might have expected as many 
as twenty on the basis of averages for the entire galactic belt. The three most remote c 
stars are at / = 166°, b = +7°, all with apparent magnitudes between 8.0 and 8.5 and 


"| Lick Obs. Bull., No. 420, 1930. 

2 Harvard Ann., Vol. 106, No. 2, 1938. 

18 Ap. J., 97, 93, 1943. 

"4p. J., 101, 309, 1945; Harvard Reprints, No. 274. 
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spectral types B9 and AOp. Allowing for a possible absorption of 0.5 mag/kpce (see be- 
low), we place them at a distance of about 2600 parsecs from the sun; there are no known 
c stars beyond the 3000-parsec limit. 

Data for faint cepheid variables are very incomplete for this part of the sky. Dr. 
Payne-Gaposchkin has kindly permitted us the use of a manuscript catalogue of known 
galactic cepheids. For the interval / = 160°-180°, b = —10° to +-10°, she lists twenty-one 
long-period cepheids and three cluster-type variables. Her derived distances, corrected 
for an absorption of 0.7 mag/kpc, place all but eight of these at distances less than 3000 
parsecs from the sun. Five of these eight stars are south of the galactic circle in Orion, 
where the average photographic absorption is in excess of the average value assumed by 
Dr. Payne-Gaposchkin. The remaining three stars are VW, VZ, and WW Mon, discov- 
ered by Hoffmeister, with estimated distances of 3500, 5500, and 3700 parsecs. These 
three stars, all near? = 169°, 6 = +1°, require further study and are the only known ob- 
jects to lie beyond 3000 parsecs from the sun. 


TABLE 1 
DISTRIBUTION OF FAINT B STARS 


B3-B5 
Gem-Ori Gem-Ori | | Gem-Ori Gem-Ori 
l= 165° 
b= 0° 


9.0-10.0)... 
10.0-11.0 


It will be of considerable general interest to pursue further the problem of the extent 
of our galactic system in the direction of the anticenter. Spectral absolute-magnitude sur- 
veys to faint limits provide the most obvious means for extending the search; infrared 
surveys, like those now under way at the Warner and Swasey Observatory, should be 
especially valuable. None of the known novae in our section of the Milky Way are be- 
yond the 3000-parsec limit, but a systematic search for faint novae in some of the more 
transparent low-latitude fields might be profitable. But, obviously, the greatest interest 
would attach to a search for faint cepheids, preferably to be combined with a study of 
the color excesses of the variables of known period. Star-count surveys to the faintest at- 
tainable magnitude limits will be helpful in the analysis for the distance to which our 
Milky Way system stretches beyond our sun. 

To complete the general picture of what we may learn about galactic structure solely 
on the basis of published material, we should turn now to an over-all description of the 
absorption characteristics of the section on the basis of the color excesses by Stebbins, 
Huffer, and Whitford.!° It will, however, be more profitable to discuss this material in 
conjunction with our own color data. 


II. SPECTRAL TYPES AND COLOR INDICES 


The spectral types for the stars for which the color indices have been measured are the 
published ones from the Henry Draper Extension.? The color indices for the several fields 
were determined mostly from plates taken with the 8-inch Ross (IR) camera of the Oak 
Ridge station. The blue series were made on Cramer Hi-Speed plates without a filter; the 
red series on Eastman 103a-E plates with a ciné-red filter. The principal observer with 
the IR camera was Guillermo Haro, now of the Tacubaya Observatory, Mexico. 
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As noted below, a group of special check plates was taken with the 24-33-inch Jewett- 
Schmidt telescope (observer: Bok) and another with the 16-inch Metcalf telescope (ob- 
server: Henry A. Sawyer). 

The procedure of measurement and reduction follows closely that described in the pa- 
per by Bok and Mrs. Rendall-Arons.' The principal difference is that practically all the 
photometry is based on direct polar comparisons and that hardly any use was made of 
comparisons with the C regions. This proved to be a fortunate circumstance, since it was 
shown in a study by Bok and Miss Sawyer" that the present standard red magnitudes 
of the C regions are in error by perhaps as much as 0.3 mag. The standards used in the 
present study are again those of Trans. J.A.U., 1, 71, 1922, for the photographic magni- 
tudes and those of the Harvard Ann., Vol. 89, No. 5, 1935, for the photored magnitudes. 
The recent checking of the basic red magnitudes by Nassau and Burger" supports the 
standards used in the present research. 

Corrections for differential extinction were generally quite small; they were applied 
as in the earlier investigation. 

Special care was taken to avoid the introduction of errors arising from color equation. 
As a standard precaution, all normal polar-sequence stars are plotted as black dots in 
our graphs, whereas the “‘r” stars are shown by red dots. This precaution was applied 
equally to photographic and photored polar reduction curves. 

As before, the principal interest attaches to the color excesses of the stars of spectral 
types B and A. In agreement with the earlier investigation, we have assumed normal 
color indices as follows: 


Inside the circle on each plate free from distance correction, we have measured, in 
addition to all B and A stars, a number of stars of later spectral types. It was hoped that 
these stars would provide an additional check on possible zero-point errors in our color 
indices; such errors could be produced by photometric nonuniformity of the sky during 
exposures and similar effects. While the observed average color indices of these stars are 
listed in the tables for each field, they are, unfortunately, of little value until reliable 
normal color indices for unreddened late-type stars become available. For a rough check 
on possible systematic errors affecting the color indices of all stars in a field, we have 
used the following normal blue-red color indices: 


KO0-K2 
K5-M 


These are more positive than the means for the Harvard Standard Regions, as given by 
Mrs. Gaposchkin" by about the amount of the correction suggested by Bok and Miss 
Sawyer." In deciding on these average values, we have, furthermore, attempted to take 
into account the difference in color between late-type giants and dwarfs. Because of the 
obvious uncertainties in the normal colors, the observed colors for the F—M stars have 
been used only in a general rough check of the zero-point of the colors. The difference in 
observed color index between the B and A and the F stars is always, however, a depend- 
able indicator for space reddening. 


16 Harvard Bull., No. 918, 1946. 
16 4p. J., 103, 25, 1946. 
7 Harvard Ann., 89, 113, 1935. 
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FIELD 1 
Center (1900): a = 6°24™; 6 = +7°40’; 1 = 171°; b = 0° 


Field 1 is in a smooth and regular part of the Milky Way. Directly south of it lies the 
complex cluster NGC 2244 and the associated emission nebulosity NGC 2237 (12 Mon 
belongs to this group), for which we have already found an indicated absorption of the 
order of 1 mag/kpc to a distance of 1300 parsecs. We mentioned further that the absorp- 
tion amounts to only 0.5 mag/kpc in the clear part north of NGC 2244. 

Since Field 1 occupies a key position (low latitude, no obvious overlying obscuration) 
in this section of the Milky Way, we have thought it worth while to check the color in- 
dices for the HDE stars from independent series of polar comparison plates taken with 
the 8-inch Ross (IR), the 16-inch Metcalf (MC), and the 24-33-inch Jewett-Schmidt (J) 
telescopes at the Oak Ridge station. 

The number of comparisons, the probable error of a magnitude determination from 
a single plate, and the corresponding error for the mean of the plates of a given series are 
shown in Table 2. In the last line of this table we have listed the systematic differences 
from the mean for the plates of each series. Several interesting conclusions may be de- 
rived from Table 2. It is to be noted that the Jewett-Schmidt telescope gives the poorest 
results, as far as accidental errors are concerned. The reason is probably that very short 
exposure times, 1 minute and less, had to be used to obtain images of the field stars that 
were not too dense for measurement. Polar comparisons of 3-10 minutes and longer yield 
considerably smaller accidental errors than those of Table 2. 

The smallness of the accidental errors for the MC telescope should not be taken at its 


Fic. 1.—This chart represents the section of the Milky Way at the junction of Gemini, Monoceros, 
and Orion. The epoch of the chart is 1900; the centennial precession for the center of the chart is: +5.5 
min. in right ascension; —3’ in declination. The positions of all HD stars to photographic magnitude 7.5 
are shown as black disks. The O and B stars with known color excesses by Stebbins, Huffer, and Whitford 
(Mt. W. Contr., No. 621, 1940) are indicated by black disks surrounded by a circle for the stars of magni- 
tude 7.5 and brighter and by open circles for the fainter stars. The cepheids VW Mon, mag. 13.2-14.2, Per. 
1.53 day; VZ Mon, mag. 14.3-15.2, Per. 5.08 days; WW Mon, mag. 12.6-14.1, Per. 4.66 days and the nebu- 
lar variable R Mon are shown. The stars y Gem, & Gem, a Ori, u Ori, E Ori, v Ori, «Mon, S (= 15) Mon and 
12 Mon are marked on the chart. The large numbers 1-7 mark the centers of the areas used in the study 
of spectra and color indices. LF8 is the center for which spectra and color indices have been determined at the 
Warner and Swasey Observatory (Ap. J., 109, 139, 1949). Galactic clusters are indicated by broken cireles 
with diameters equal to the diameter of each cluster according to Trumpler (Lick Obs. Bull., No. 420, 
1930); they are identified by small letters a to p, as listed below. Diffuse nebulosities, according to Ceder- 
blad (Lund Medd., Ser. II, No. 119, 1946), are presented as small squares for diameters less than 6’ and 
squares of the dimensions given by Cederblad for the larger ones; they are identified by small letters 
g to z, as listed below, with C indicating continuous nebular spectra and E emission spectra. The more 
conspicuous patches of dark nebulosity are shown by full-drawn lines of irregular shape. 


Galactic Clusters Diffuse Nebulosities 
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face value. The few MC plates were taken within a very short time of one another, and 
the listed error is a measure principally for the internal agreement of the measures— 
which is excellent. Only the very best nights were used for the MC observations, whereas 
for the IR and J we had to be satisfied with a mixture of good and excellent nights. We 
conclude that our Schmidt-type telescope had better not be used in the future for polar 
comparisons of very short exposure times. 

The measured area had a radius of about 2°. No differential systematic errors were 
found depending on a star’s position in the field, and almost complete freedom from “‘dis- 
tance correction” seems indicated. 

The data shown in Table 2 depend on thirty-four stars measured in common on all 
IR, MC, and J plates. Many more stars were measured on the IR plates, two hundred 
and forty-eight of them in all. The resulting mean color indices are given in Table 3. 


TABLE 2 


No. of plates.......... 
Av. p.e. per measure. . .| 
Av. p.e. per star 

Syst. diff. from mean... 


TABLE 3 
MEAN OBSERVED COLOR INDICES IN FIELD 1 


SPECTRAL TyPE 


B8-A0 A1-A7 


| +0.02(1) | 

.09(6) .38(2) | +1.17(4) |........... 

+0. 32(1) | — .01(5) .23(5) | +0.20(3) | +1.40(4) | +1.81(3) 
| — .14(1) .15(10) | — .06(9) .20(2) | +0.52(2) | +1.27(3) | +1.58(2) 
—0.06(1) | t -06(16) | t .35(9) | +0.60(8) +1.30(4) | +1.85(4) 
08(19) | + .06(26) | + .35(16) | +0.88(8) | +1.23(8) | +1.61(5) 
to 09(7) | +0. 13(21) | to. 49(5) | +0.68(3) +1.65(3) 


win 


An inspection of the mean color indices for stars of spectral types F-M would seem to 
indicate that the observed color indices in Field 1 are too small by at least 0.10 mag. 
This conclusion is, however, not confirmed by Table 2, where the IR color indices are, if 
anything, slightly more positive than the mean for the MC and J. We have, therefore, 
not applied any systematic correction to the average color indices of Table 3. We have al- 
ready commented on the lack of reliable normal colors for the F—-M stars and about the 
possible effects on the means of the unknown ratio of dwarfs to giants. 

We conclude from Table 3 that the Al—A7 stars give no indication for the presence of 
space reddening, with the possible exception of the most distant group. This means, es- 
sentially, freedom from reddening effects to a distance of 1000 parsecs. The fainter B8- 
AO stars show definite indications for space reddening, with a total absorption of 0.4 mag. 
at 1200 parsecs and 0.7 mag. at 1500 parsecs. 
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Fic. 2.—The galactic cluster NGC 2244 and associated diffuse nebulosities NGC 2237 and NGC 2238. A3x 
enlargement of a 75-minute exposure on 103a-O backed plate taken with the 24-33-inch Jewett-Schmidt telescope 
of the Oak Ridge Station, 
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The results from the early B stars confirm, with two exceptions, the smallness of the 
interstellar absorption. The positive color index of a single B5 star (mp, = 9.77) need 
hardly concern us, since it is offset by much smaller excesses for two fainter earlier B 
stars. Our own color excesses of early B stars suggest a total photographic absorption of 
0.9 mag. at 3000 parsecs. The other exceptionally high absorption is found for two B 
stars observed by Stebbins, Huffer, and Whitford !° which would indicate a total photo- 
graphic absorption of 0.8 mag. at 900 parsecs. Since these two stars are not far from NGC 
2244, the somewhat greater absorption is not surprising. The four remaining excesses by 
Stebbins, Huffer, and Whitford lead to the following absorptions: 


We conclude, therefore, on the basis of all available evidence, that the total photo- 
graphic absorption derived from space reddening for the direction of Field 1 is repre- 
sented by 

A,, = 0.3 mag /kpe 


to a limiting distance of 3000 parsecs. Our own as well as other colors show that the ab- 
sorption for this direction is not so high as 


Aye = 0.5 mag /kpc ; 


but, at the same time, the presence of some reddening seems definitely established. 


FIELD 2 
Center (1900): a = 6"13™; 6 = +12°3;] = 166°; b = 0° 


This field lies directly at the galactic circle. The field has some irregular absorption 
overlying it, but the irregularities are slight. Inside the limits covered by our color sur- 
vey we find the galactic cluster NGC 2194, which, according to Cuffey" is highly red- 
dened; Cuffey finds a reddening corresponding to a total absorption of 2.7 mag. at 1200 
parsecs. 

Our basic plate material for the color-index survey is considerable, comparable in ex- 
tent to that for Center 1. It consists of polar comparison plates made with the 8-inch 
Ross (IR), 16-inch Metcalf (MC), and the 24-33-inch Jewett-Schmidt (J) telescopes at 
Oak Ridge. The number of comparisons, the probable error of a magnitude determina- 
tion from a single plate, and the corresponding error for the mean of the plates of a given 
series are shown in Table 4. In the last line of this table we have listed the systematic dif- 
ferences from the mean for the plates of each series. There is no marked variation with 
magnitude in the systematic differences, and we conclude that the color indices on the 
IR, MC, and J system yield essentially similar results; we note again that none of the 
systematic differences is as high as 0.10 mag. 

We also note here the smallness of the size of the derived accidental errors of the mag- 
nitudes derived from the MC plates and the relatively large probable errors for the Jew- 
ett plates. This confirms our conclusions that polar comparisons depending on very short 
exposures with the Schmidt camera produce rather large accidental errors. The IR mag- 
nitudes have probable errors intermediate between those from the MC and J plates. 

The results of Table 4 depend on only 34 stars measured in common on all IR, MC, 
and J plates. Altogether, 185 stars were measured on the four blue and the five red IR 
plates. The resulting mean color indices are given in Table 5. 

A comparison between the average norma! color indices for the various spectral in- 
tervals and the observed means shows that the averages for the F, G, and K stars are 
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about as we expect them to be; this result confirms the absence of systematic errors in 
excess of 0.10 mag. in the color indices derived from the IR plates. 

The color excesses derived from the IR polar comparisons show that considerable space 
reddening is present in Field 2. On the basis of the observed color excesses for the BO-BS, 
B8-A0, and A1-A7 stars, we find the total photographic absorption to vary as follows 


with distance: 
1000 1500 (parsecs) 
1.5 2 (mag.) 


This is somewhat smaller than the total absorption found by Cuffey for NGC 2194, 
but somewhat larger than the absorption indicated by the five B stars measured in this 
section by Stebbins, Huffer, and Whitford.'® 


TABLE 4 


No. of plates 

Av. p.e. per measure.........| 
Syst. diff. from mean 


TABLE 5 
MEAN OBSERVED COLOR INDICES IN FIELD 2 


SPECTRAL 


FO-F8 GO-G5 
0.35(1) | +0.27(2) | +1.08(1) 
+ .34(1) | 


+ 
+ 
+1.16(2) 
+ .03(1) -04(3) | + . 100: +2.21(1) 
+ .41(3) + +1.82(3) 
+ 
+ 
+ 


| 


+0.75(1) 


“SIN 


+0. 36(3) +1.83(2) 
50(15) 


0. 35(8) -84(5) | +1.63(2) 


It is of interest to compare our results with those obtained recently by Dr. Sidney W. 
McCuskey for Field LF8 in the Warner and Swasey Observatory survey.'* The total ab- 
sorptions derived from the Warner and Swasey color excesses for LF8 (a = 6"14™; 
6 = +14°0), which were communicated to us in advance of publication by Dr. Mc- 
Cuskey, are: 

Distance;..6.:046.2. 500 1000 1500 2000 2500 (parsecs) 
¢ 0.3 1.0 2.0 2.4 2.4 (mag.) 


The agreement between these two wholly separate determinations (different spectral 
classifications; different telescopes used for determination of magnitudes and color in- 
dices; wholly independent reduction of results) is quite satisfactory. 


18 Cf. McCuskey and Seyfert, Ap. J., 106, 1, 1947. 
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We conclude that, according to the various lines of evidence, the total photographic 
absorption in this field can be represented by 


Apg = 1.2-1.5 mag/kpe 


to a limiting distance of 1500 parsecs, with McCuskey’s results indicating a total absorp- 
tion of 2.4 mag. at 2500 parsecs. 

It is necessary to comment on the dissimilarity of space reddening for two fields so 
relatively close to each other as Fields 1 and 2. They are only a little more than 5° apart 
in galactic longitude and both are directly on the galactic equator; but, in spite of this, 
Field 1 shows only a total absorption of less than 1 mag. at 2500 parsecs, whereas for 
Field 2 the corresponding absorption is more nearly of the order of 3 mag. Since this rath- 
er surprising result is given by color excesses derived with a variety of techniques and 
telescopes, there can be little doubt about its reality. Apparently, Field 1 happens to be 
in a direction not affected by the Taurus and Orion complexes of dark nebulosity, where- 
as Field 2 is under the influence of extensions of these large clouds of galactic obscura- 
tion. 

FIELD 3 


Center (1900): a = 6 = +1499; ] = 165°; b = +3° 


FIELD 4 
Center (1900): a = 5857"; 6 = +1591; 1 = 162°; b = —2° 


These two fields are a little distance north and south of the galactic circle, at very 
nearly the same galactic longitude as Field 2. 

Field 3 lies in a rich and fairly uniform part of the Milky Way, bordered by regions 
with irregular obscuration. It is marked by an almost complete absence of early B stars 


and a scarcity of B8—A0 stars, but the number of A1—A7 stars is fairly large. 

The color indices of 196 stars in Field 3 were determined from five blue and five red 
polar comparison plates made with the 8-inch Ross (IR) camera. With one exception, the 
systematic differences, individual plate minus mean of five, were small; one blue plate 
yields magnitudes that are systematically too large by 0.33 mag. as compared to the 
mean of the five blue plates. We have used in our analysis the uncorrected color indices, 
and our derived color excesses could be too high by 0.10—-0.15 mag. The observed mean 
color indices are shown in Table 6. The color indices for the late-type stars are close to 
normal,and we have no reason forthe application of a systematic correction to the observed 
average color indices of the early-type stars. The derived color excesses indicate a total 
absorption slightly in excess of 1 mag. at 1000 parsecs. This is confirmed by two near-by 
reddened stars (650 parsecs) for which Stebbins, Huffer, and Whitford'® have measured 
excesses, but it is to be noted that, at a somewhat greater distance from tke center of 
Field 3, the same authors have an unreddened star at very nearly the same distance. Our 
single early B star yields an absorption of 0.8 mag. at 1500 parsecs. We conclude that 
the color excesses for Field 3 indicate a slightly smaller absorption than for Field 2, with 


A yg = 1.0 mag/kpc 


being a fair average. 

Field 4 at 6 = —2° lies ina rich field, again bordered by obscuration. A long-exposure 
Schmidt plate shows that there is some irregularity in the surface distribution of the 
faint stars, with the largest surface densities for the eastern half of the region. 

The color indices of one hundred and eighty-eight stars in Field 4 were determined 
from four blue and four red polar comparisons with the 8-inch Ross (IR) camera. The 
plates showed good agreement. The mean observed color indices are shown in Table 7. 
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The color indices of the late-type stars again do not suggest the need for systematic cor- 
rections. The data for Field 4 seem to be quite consistent. The three faint early B’s in 
the field are all quite reddened, with an absorption of at least 2 mag. indicated at 1500 
parsecs. This result is confirmed by the value of 1.4 mag. at 1000 parsecs derived from 
the B8—AO and A1-A7 stars. The Stebbins, Huffer, and Whitford'® color excesses for this 
section are generally small, but the stars are bright and near by. We shall accept the 
value of A,, derived from our own color excesses: 


Ape = 1.4 mag/kpe 
to a limiting distance of 1500 parsecs. 


TABLE 6 
MEAN OBSERVED COLOR INDICES IN FIELD 3 


SPECTRAL TyPE 


Al-A7 


.45(1) 
+ .21(1) 77(1) 
45(3) | 


‘66(1) | +1.31(1) 
+ : .53(6) | +0.70(3) | 
37014) | 1.0106) | 
| 40.4202) | + 137 + .55(19) +0.89(11) 
| 40.4003) | +0.69(1) +0.67(2) | 


O7(1) 


| 


| 


TABLE 7 
MEAN OBSERVED COLOR INDICES IN FIELD 4 


SPECTRAL TYPE 


G0-G5 


40.811) | 4+0.95(2) | 

-00(1) 7201) | +0.96(2) 
+ -52(3) | +0.63(1) 

.13(6) | + .21(4) 738) +1.01(3) 

+0.46(2) | + .42(13) | 5719)  +0.98(7) 

| +0. 73(1) | -26(4) + .41(23)) + .47(3) | +0.87(9) | 
ee 39(3) +0.43(13) 46(2) | +1.17(3) 


_ 


ATH 


— 


| 


FIELD 5 
Center (1900): a = 633"; 6 = +16°8;1 = 164°; b 


FIELD 6 
Center (1900): a = 6°40"; 6 = +15°0; 1 = 167°; b 
FIELD 7 
Center (1900): a = 646™; 6 = +12°1;] = 170°; b = +7° 


These three fields are in the zone of galactic latitude +5° to +10° and are free from 
strong obvious overlying obscuration. Field 5 is probably affected by some irregular 


2 
: | po-Bs |  B8-A0 ms «=| | Ko-K2 | KS-M 
; 
| Bo-Bs | Bs-Ao | ai-A7 | FOFS | mum | Ko-K2 | Ks-M 
10 +1.55(3) | +1.90(2) 
10 +1.59(4) | +2.08(1) 
a +1.47(3) | +2.12(6) 
= +6 
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overlying absorption, which, however, seems limited to a considerable number of very 
small areas. There is some suggestion of star deficiency in the northernmost parts of the 
field. Fields 6 and 7 are both quite uniform, and the stellar distribution on long-exposure 
Schmidt plates gives no indication of overlying irregular obscuration. 


TABLE 8 
MEAN OBSERVED desain OR INDICES 


SpecTRAL 


Field 5 


+0. 
+0. 


07(1) 


.04(3) 
.02(6) 
+ .29(9) 
+0. 


+0. 23(2) 


+ .14(3) 
— .02(3) 
+ .06(13) 
+ .18(15) 
+ .33(19) 
+0.45(7) 


+0.01(5) 
+ .12(4) 
+ .16(2) 
— .04(12) 
+ .03(18) 
+ 06(19) 
+0. 29(2) 


32(1) 


+0. 63(1) 


+0. 26(4) 
+0.25(4) 
+0. 32(1) 
+0.51(12) 
+0.75(9) 
+1.25(1) 


+1.48(1) 
+0.92(2) 
+1.04(2) 
+1.21(2) 
+0.88(4) 
+0.85(3) 
+0. 64(13) 
+0.55(3) 


Field 6 


+0.69(1) | 
+40.45(5) 
+40. 30(1) 
+0. 39(4) 
+0. 39(18) 
+0.63(11) 
+40.60(13) 


| 
| 


+0.81(3) 
+ .61(2) 
+ .91(6) 
+ .50(4) 


+1.08(1) 


+0.44(11) |.... 


0.62(2) 
(: 


tt+tt+++t 


0. 80(1) 


+1.56(1) 
+1.79(1) 
+2.20(2) 


+2.01(5) 


The color indices of Table 8 are all determined from blue and red polar comparisons 
taken with the 8-inch Ross (IR) camera. In the case of one or more large residuals, special 
check plates were taken on which remeasures were made of a selected sequence of about 
thirty stars. The final derived color indices are shown in Table 8. About two hundred 
stars were measured in each field. 

The color excesses appear to be relatively small in all three fields. There are no indi- 
cations that the extended Taurus-Orion obscurations reach to as far as galactic latitude 
+7°; they are apparently confined to latitudes south of +3°. The average color indices 


flee B8-A0 | A1-A7 | G0-GS | K0-K2 Ks-M 

8.25..........{ +0.08(1) | | | 4 
: $.75........... — .2Q) 
| | +0.90(1) |............ 
9.25..........| | | +1.42(1) | 
10(7) | + 78(4) 41.59(3) | +1.69(1) 
10(9) + .36(9) +1.27(5) | +1.69(4) 
10.75. . O8(7) | | + .28(11) +1.28(5) +1.42(3) 
| + .08(11) .98(6) | +1.43(3) | +1.42(3) 
+ .25(8) .57(2) | +1.07(10) | +1.62(4) 
+ .25(14) | 64(5) +1.50(2) | +1.44(8) 


34 BART J. BOK, MARGARET OLMSTED, AND B. D. BOUTELLE 


of the F-M stars vary irregularly; but, keeping in mind the small numbers of stars in 
each average and the uncertainties in the spectral classification and in the admixture of 
giants and dwarfs, we conclude that in only one case is there a justification for the ap- 
plication of an over-all systematic correction; for Field 6 we have applied a systematic 
correction of +0.10 mag. to all excesses derived from Table 8. 

In Field 5 the A1—A7 stars show no effect of reddening to m = 11.0. The B8—A0O stars 
begin to show evidence for reddening at m = 10.0, with a maximum indicated value of 


A yg = 0.8 mag/kpc 


toa distance of 1000 parsecs. There is only a single early B star in the field; the indicated 
total absorption at 2100 parsecs is 0.6 mag. This would support a considerably smaller 


value of Ap, than the one given above. 
The observed excesses in Field 6 are quite small. The brightest B8-A0 and Al—A7 


stars indicate some near-by reddening, but this is not confirmed by the more distant 
stars. There are no data from clusters or early B stars. A value of 


= 0.5 mag/kpc 


to a limiting distance of 1000 parsecs is a good compromise. There are no data for greater 
distances. 

In Field 7 we again have definite evidence for slight reddening. The B8—A0 and A1l-A7 
stars indicate the presence of an absorption of about 0.8 mag. at 1000 parsecs. There are, 
surprisingly, six early B stars in this field with estimated distances of 1200-1900 parsecs. 
The color indices of none of these stars a previously been measured, but from our re- 


sults we derive a value of 


to a limiting distance of 2000 parsecs. 
There appears to be no reason for separate analysis of these three fields. In what fol- 


lows, we shall assume that 
= 0.5 mag /kpc 
to a limiting distance of 2000 parsecs describes average conditions for the section at 
1 = 163°-173°; b = +5°-+9°. 
We may add here that Hubble’s counts of faint galaxies'® indicate a total photographic 
absorption of 2.4 mag. at 6 = +5° and of 1.5 mag. at b = +10°. 


FIELDS 8-11 


The original program for the section under investigation included four more fields. 
Every one of these fields was found to be affected by irregular absorption, and no at- 
tempt was made to complete the analysis for these fields. A brief summary of the pre- 
liminary results follows. 

Field 8 (a = 5"52™; § = +8°;1 = 167°; 6 = —7°).—Long-exposure plates made with 
the Jewett telescope show this field to be affected by considerable overlying absorption. 
The “wishbone” dark nebula, shown nicely on Chart 33 of the Ross-Calvert Allas of the 
Milky Way, in part overlies this field. The dark nebulosity is apparently an extension of 
the Taurus complex and is presumably relatively near the sun. 

The color excesses for this field indicate the presence of considerable reddening. At 
300 parsecs the total photographic absorption derived from the color excesses already 
amounts to 0.7 mag.; the total absorption at 1000 parsecs is of the order of 2 mag. 

Field 9 (a = 6'O2m: 6 = +10°5;1 = 166°; b = —3°).—This field shows a fairly uni- 
form surface distribution of the stars, but it is clearly located on the edge of a region of 


19 Ap. J., 79, 8, 1934; Mt. W. Contr., No. 485. 
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considerable obscuration. The measured color indices give evidence for considerable red- 
dening. At 400 parsecs from the sun the indicated photographic absorption is already 
equal to 1 mag., whereas the total photographic absorption at 1000 parsecs is about 
2 mag. With such a high absorption, the field is not a suitable one for detailed study of 
the density distribution. 

We have already noted above that the cluster NGC 2141, estimated by Trumpler to 
be at a distance of 4860 parsecs, falls within the limits of this field. If our estimated 
photographic absorption is correct, this distance should be reduced to about half this 
value. 

Field 10 (a = 6'10™; § = +2°5;1 = 174; b = —5°).— The measured color indices 
for this field indicate considerable reddening near by, with 1 mag. of photographic ab- 
sorption at 550 parsecs. The available evidence does not indicate much of an increase of 
absorption with distance beyond 550 parsecs. 

Long-exposure photographs give evidence for irregular near-by obscuration overlying 
most of the field. There are some stars with photoelectric colors by Stebbins, Huffer, and 
Whitford"® bordering on this field. These give evidence of little reddening, but the long- 
exposure photographs show them to be in regions free from obvious overlying obscura- 
tion. The selection of the center for color studies was not a happy one, and we shall 
therefore not include this field in our density analysis. 

Field 11 (a = 6435"; § = +9°;1 = 171°; b = +3).—This field is probably the most 
interesting one of the whole group, but it happens to border on the great dark nebula 
and is a region of highly irregular obscuration. S Mon and the associated cluster and 
nebulosities are in this field, and the three remote cepheid variables VW, VZ, and WW 
Mon are near it. 

The results from the available data on color excesses are inconclusive. Three of the 
stars with excesses measured by Stebbins, Huffer, and Whitford!® are unreddened 
(limit 1000 parsecs), whereas a fourth star would indicate 2 mag. photographic absorp- 
tion at 720 parsecs. Our own color indices indicate an average absorption of the order of 
1 mag/kpc, but there are indications of considerable fluctuation over the field. This field 
should be studied as a part of a more detailed analysis of the distribution of spectra, 
colors, and related star counts for the regions bordering the dark nebula of Wolf. We 
plan to undertake this study in the near future. 


III. DENSITY FUNCTIONS 


From the data on the excess reddening of the stars in the Henry Draper Extension, we 
have succeeded in deriving approximate values for the photographic absorption for 
seven fields in the section under investigation. It is now of interest to analyze the density 
distributions for the various spectral subdivisions in this same section of the Milky Way. 

In 1945 Bok and Rendall-Arons! presented the results of a similar analysis for a field 
at / = 185°, not far from the area under investigation. A complete analysis of the present 
material will be postponed until the star counts for the section by Father James Con- 
nolly, S.J., become available, but it is of interest now to make a comparative study of 
the stellar distribution for three areas in the present section and for “Region D” of Bok 
and Rendall-Arons. For purposes of comparison, we reproduce here, as Table 9, Table 13 
of the paper by Bok and Rendall-Arons. 

Our measured color indices refer in each case to relatively small fields of fairly uniform 
surface distribution, which appear to be similar to their surroundings in structure and 
in freedom from overlying localized obscuration. The derived absorptions may therefore 
be considered as applying to areas larger than those actually covered by the color survey. 
Furthermore, our color analysis has indicated the existence of fairly similar conditions 
in several neighboring fields -overed by the color survey; we found in the preceding sec- 
tions that the results are very similar for Fields 2, 3, and 4 and also for Fields 5, 6, and 7, 
but we found that Field 1 had to be treated separately. 
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The enlarged regions for which the counts with respect to spectral subdivisions were 
actually carried out and the assumed absorption characteristics are tabulated in Table 
10. For comparison purposes we have also given the data of Bok and Rendall-Arons for 


TABLE 9 


SUMMARY OF STAR DENSITIES (NUMBER OF STARS PER 1000 
CUBIC PARSECS) IN MONOCEROS REGION D 
(Corrected for A,,=0.6 Mag/Kpc) 


SPECTRAL TyPE 


| gGO-gM 


TABLE 10 
THE ENLARGED REGIONS USED FOR THE DENSITY ANALYSES 


REGION 


Cc Mon ‘‘D”’ 


Assumed value of A,ye.........| 0.3 mag/kpc 1.4 mag/kpe | 0.5 mag/kpc | 0.6 mag/kpc 
Limit to which Aj. applies... | 3000 pes 1500 pes 
Basic color data from fields. . 1 2-3-4 | 5-6- 
Limits of a andé..... 19 Gh37m—6b57™ 
+i—+10"" | | 
and and 
6"20™-6b32™ =| 620™—6537™ 6525™—7h05™ 
+4°-+10° +10°-+16° | +13°-+18° —7°-+3° 
and 
6>32™Gh49™ 
+4°-+8° 
Average gal. long. and lat......| 1=171°; b=0° = 165°; b=0° | 1=186°;b=+7° 1=182°; b=0° 


| j 
Total area (sq. deg.).........| 3 100 


‘Region D.” It will be noted that Regions A, B, and D are all centered upon the galactic 
equator, whereas C lies 7° to the north of it. 

Counts were made in the Henry Draper Catalogue and Extension between the limits 
of a and 6 indicated in Table 10. To supplement the data for the brighter stars and to 
provide the much-needed link between the average distribution for the immediate 
vicinity of the sun and the stellar distribution at greater distances, we have again used 


~ = - - fot. 
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B8-A0 AI-A7, | FOF? | FSF) 
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the counts from a paper by Seydl” for the stars brighter than m = 7.0 in the zone of 
galactic latitude 0° to + 10°. All counted numbers were reduced to a unit half-magnitude 
interval and to a total area of the sky of 100 square degrees. The derived values of 
log A’(m) were plotted and smooth curves drawn to represent the combined results from 
Seydl’s counts and our own HD and HDE counts. The results are shown in Tables 11-16; 
the magnitudes used are photographic. 


TABLE 11 
DERIVED VALUES OF LOG aXe) FOR THE B8-AO STARS 


MAGNITUDE 


| 


| 
ss| | 6s | 10 | 75 


0.50) 0.78) 1.06 1.29 1.50. 1.72, 
.20, 0.50, 0.78) 1.14, 1.43, 1.66) 1. 
22,0. 51. 0.77) 0.98, 1.20, 1.40) 

0. 54 0.80) 1.02) 1. 30 56) 


TABLE 12 
DERIVED VALUES OF LOG ithe FOR THE Al—A7 STARS 


9.11) 9.50 
9.11 9.50 
9.10, 9.49 
9.11 9.49 


SRB! 


MAGNITUDE 


| 10.0 | 10.5 


32 
30 


30) 9.47, 9-63 9.82 


0.32) 0.59, +0.89 
.57} — .82 
+ .78 

58) +0.82 


10 
08 


9.48) 0.65 9.86 


9.49 9.70 


1.70) 
1.70 
1.70 


1.95, 2.16) 2.36) 2.53 
1.98) 2.27 2.52) 2.56 
2.00 


1.87 


2.25} 2.43} 2.55 


| 9. 9.87 hee 0. | 2.44) 2.63 


TABLE 13 
DERIVED VALUES OF LOG A’(m) FOR THE FO-F8 STARS 


MAGNITUDE 


REGION 
65 | 7.0 | 75 | 80 | 8s | 9.0 


os 
0.68) 0.92) 1. 18) 1.48, 1.78 2.08 2.30 2.42 
.54 1.17) 1.48) 1.74) 1.94) 2.08) 2.18 2.20 
18) .48| 1.04) 1.30; 1.59 1.91) 2.14) 2.31) 2.43 
0.42} 0:6] 0.91) 1-16 1-66] 1.90) 2.2 
| | | | 


11.5 


0.18 


B8-AO Stars.—Table 11 shows that for a distance of 1000 parsecs from the sun the 
densities in D are considerably in excess of those in A, B, and C. We note, for example, 
that, in spite of a larger value of Ay, in D as compared to A, the value of log A’(m) at 
m = 11.0 is greater for D than for A. Whereas in D we previously found average densi- 
ties 0.6 © at 1000 parsecs, we have in A an average density of only 0.4 © at 1000 parsecs. 
The computed space densities in B agree well with those for A. The effect of the differ- 


20 Pub. Nat. Obs. Prague, No. 6, 1929. 


| 
Bs | 4.0 | 4.5 | 5.0 8.0 | 8.5 | 9.0 | 9.5 | 10.0 | 10.5 | 11.0 F 7 
9.86 0 1.91, 2.09, 2 2.38) 2.51 
9.86, 0 1.98} 2.11) 2 2.19) 2.00 
9.86 0 1.75] 1.92) 2.15) 2.19 
D. 0 2.03) 2.22) 2 2.58} 2.75 
| 
Gion | | | | | | } | | 
| 4.0] 4.5 | 50] 5.5 | 60] 6.5 7.0 8.0 | 8.5 | 9.0 9.5 
| 9 1.08} 1.37 
ay 
| : 
9.19} 9.44] 9.69) 9 
9:19] 9:45 9.69 9 
D............, 9.19} 9.44] 9.69 9 
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ence in absorption for A and B is shown nicely in the small value of log A’(m) for 
m = 11.0. The effect of galactic latitude is noticeable in C. The average density at 1000 
parsecs equals 0.3 © or less, with a continued rapid falling-off indicated. The presence 
of steep negative density gradients for the B8-AO0 stars at / = 160°-175° is established. 


TABLE 14 
DERIVED VALUES OF LOG A‘(m) FOR THE GO-GS5 STARS 


MAGNITUDE 


| | ] | | | | | | 
i | | | | | | | | 
5.0 | 5.5 | 6.0 | 65 | 7.0 | 7. S| 9.0 98 | 10.0) 10.5 11.0 | 11.5 


| 8.70, 8.99 9.28 9.54 9.85 0.16) 0.49, 0.83 1.15 1.46 1.72 
8.83 9.08 9.33 9.61 9.88 .52) .83 1.14 1.46 1.70. 
05) .32 


TABLE 15 
DERIVED VALUES OF LOG A’(m) FOR THE KO-K2 STARS 


MAGNITUDE 
REGION 


9.22 9.56 | 0.06 | 0.32! 0.55 0.80 1.03 | 1.27 1.511 
| 9.25 | 9.54] 9.83] .11] .36] .64] .88| 1.14] 1.36 | 1.55] 1. 
| 9.31/9.55| 9.82) .09| .94 1.22 | 1.48 | 1.73 | 1. 

1. 


(9.34 9.56 9.78 0.05 0.35 0.63 0.89) 1.18 1.46 | 1.75. 


TABLE 16 
DERIVED VALUES OF LOG A’(m) FOR THE K5—M STARS 


MAGNITUDE 


REGION 


42 | 9.69 | 9.97 | 0.25 | 0.53 0.85 | 1.12 | 
36 | 9.59) 9. 56 | 0. 
3 
40 


| 
| 
| 
| 
| 


Al-A7 Stars.—The values of log A’(m) in all four fields show a striking similarity. 
A density analysis with the aid of an (m, log w) table, corrected for absorption effects 
(which have only relatively little influence because of the small distances involved) gives 
densities of the order of 0.5 © at a distance of 500 parsecs. 

FO-F8 Stars.—The numbers of counted stars were so small that we have considered 
it advisable to combine the results for the FO-F2 and F5-F8 stars. There is still consid- 
erable local variation in the star numbers, with a maximum range of a factor of 2 in the 
counted numbers. 


1.96 2.17, 2.34) 2.51 
: 1.91) 2.44) 2.28}: 
1.91] 2.18) 2.41)... 
i 8.83) 9.32) 9.54 9.78, 0.04, 0.32 0.60 0.91) 2 

10.5 | 11.0 
75 | 1.98 | 2.23 
76 | 1.92 | 2.07 

i { 
tog 5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 | 10.0 | 10.5 | 11.0 = 

| 9.09 | 9 1.47 1.70 | 1.95 
BE D...........| 8.93 | 9.16 | 9. | 1.52 | 1.76 | 1.93 eo 

| 
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A density analysis shows that the density gradients in A, B, and C do not have quite 
so large negative values as the gradient previously derived for D. The average space den- 
sity of the FO-F8 stars for our section at 400 parsecs from our sun is, however, still less 
than 0.5 ©, thus confirming the existence of steep negative gradients found from the 
earlier investigations. 

G-M Stars.—Because of the uncertainties in the percentages of dwarfs at successive 
apparent magnitudes, we shali not give here a detailed analysis for space densities. Sec- 
tions 4 and 5 of the paper by Bok and Rendall-Arons show the kinds of results that may 
be obtained, if one makes certain assumptions about the percentages of dwarfs; the 
analysis refers there to ‘Region D.” The most striking difference between D and the new 
regions, A, B, and C, is that the new regions show a marked excess of GO-G5 stars over D. 

For the KO-K2 and K5-M stars, the numbers in A and D are remarkably similar. 
Region B is poorest in faint K—M stars, with the deficiency becoming especially great 
for the K5-M stars with m = 9-11. 

Preliminary density analyses confirm the results obtained previously for D, viz., that 
the average space densities are fairly constant over the first 500 parsecs from the sun. 
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ABSTRACT 


This paper, like previous ones of the series, is primarily concerned with the study of the fluctuation 
of the stellar luminosity function. The field under investigation has an area of 11 square degrees with a 
center at R.A. = 20500™, Dec. = +38°0 (1950). It is relatively free from space absorption, with the 
exception of about 1 square degree at its southern boundary. 

The photographic and photored magnitudes and the spectral types of 1230 stars were obtained. Stars 
of spectral type FO and later, with a few exceptions, were classified either as giants or as dwarfs down to 
the limiting photographic magnitude of 12.2. For nearly all stars brighter than 11.0 mag., four luminosity 
classes were used. 

The linear absorption coefficient in the clear region (areas A, B, and C) was derived from 683 stars 
and found to be 0.53 mag. per kiloparsec. 

The division of the stars into giants and dwarfs in each spectral group indicates that there is a rapid 
decrease in the ratio of dwarfs to giants between the FO-FS group and the G8-K3 group. In the former 
we have about 83 per cent dwarfs and in the latter about 7 per cent. 

Space densities, corrected for absorption, were computed for each spectral group in areas A, B, and C 
(Table 10 and Fig. 5). The densities were used to compute luminosity functions at distances of 100, 200, 
400, and 600 parsecs. These show very little variation with distance. Likewise, no significant differences 
exist when our results are compared with the standard van Rhijn luminosity function, with the exception 
of the range of absolute magnitude between —1.5 and +1.0. Here the maximum excess in our case in 
log is 0.3. 

The color excesses of the stars in the obscured region (area D) indicate that an interstellar cloud sets in 
at about 900 parsecs with a maximum absorption of about 0.9 mag. over and above the general absorption 
of the other three areas. A comparison between the general star counts to the seventeenth photo- 
graphic magnitude in area D and in the other areas yields nearly the same value for the total absorption. 

An analysis of the space density of the stars in area D indicates that the density for the early-type 
stars is between two and three times greater than the corresponding density in the other areas. 


I. INTRODUCTION 


The present investigation is a part of a general program of this observatory dealing 
in a broad sense with the study of the fluctuations of the stellar luminosity function in 
fields near the galactic plane. The field under discussion (LF3b) has for its center the 
point whose 1948 co-ordinates are a = 20"00™, 6 = 38°0 (J = 41°8, b = +322). It is 
circular in shape with a diameter of 3°7 and an area of 10.9 square degrees. The ‘Great 
Rift” is 4° east of its center. A portion of a parabola-shaped obscuration is included near 
its southern boundary; otherwise the area appears to be free of pronounced dark nebu- 
lae. Two other regions (LF3 and LF3a), which are part of the same general program, 
partially overlap LF3b (Fig. 1). 

This field was chosen not only on account of our desire to extend and strengthen our 
data on this interesting part of the Milky Way but also for its relatively high transpar- 
ency. A group of remote cepheid variables discovered by Baade! are in this region.’ An 
extensive survey of faint M stars in the same region is now in progress at this observa- 
tory. A portion of the survey has been completed.’ 


1A p. J., 79, 475, 1934. 

2 Oort and Oosterhoff, B.A.N., 9, 325, 1942. 

3 Nassau and van Albada, Ap. J., 109, 391, 1949. 
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Fic. 1.—A reproduction of the Ross Aflas. The circles represent: LF3 (upper left), LF3a (lower), and 
LF 3b (upper right). P Cygni is the central star in LF3. 


Fic. 2.—Division LF3b into four areas. In this reproduction north is to the left 
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II, OBSERVATIONAL DATA 
MAGNITUDES 


The photographic-magnitude sequence was determined by direct comparison with 
the North Polar Sequence on Eastman IIa-O plates. Twelve plates yielded a mean error 
for the mean magnitude of each star of +0.07. (Mean errors are used throughout this 
paper.) The photored magnitudes of the same stars were obtained by direct comparison 
with a standard photored‘ sequence at the pole. Eastman 103a-E plates and a Wratten 
No. 22 filter were used. Ten independent determinations were made which yielded for 
the mean error of the mean magnitude of each star the value +0.06. 

The magnitudes of the field stars were measured on three, and sometimes four, pairs 
of plates. In each pair, one plate was taken with the telescope east of the piers and the 
other with the telescope west. The average mean error of the adopted photographic mag- 
nitude of each field star was + 0.10, and the corresponding photored value was + 0.08. 
No serious systematic differences in the magnitudes in different parts of the field were 
observed. Nearly all the photometry was done with a graduated scale of stellar images. 
The photoelectric densitometer was used in determining the magnitudes of the brighter 
stars on plates of shorter exposure time. 

Although we were able to classify stars fainter than photographic magnitude 12.15, 
the limit of completeness of our survey is taken as 12.0 mag. The total number of stars 
available for analysis is 1230. 


SPECTRA 


The criteria for classification have been described in previous papers.® * Three or 
more plates of good quality were examined independently by both of us. These plates 
(Eastman IIa-O emulsion, exposure 30 minutes) were taken with the 4° objective-prism 
arrangement on the Burrell telescope. The width of the spectra was about 0.25 mm. 
Shorter-exposure plates were taken for study of the bright stars. In many cases we have 


examined as many as five plates, including some taken with the 2° prism, in order to 
clear up discrepancies in the classification or to overcome difficulties with overlapping 
spectra. The adopted spectral class for each star is therefore the average of several de- 


terminations. 
All stars of spectral class FO and later were classified either as giants’ (III) or dwarfs 


(V), with some exceptions when the stars were too faint or their spectra badly overlap- 
ping. For stars of spectral classes FO and later, with photographic magnitudes brighter 
than 11.0, four luminosity classes were used. 


COMPARISON OF DATA WITH OTHER INVESTIGATIONS 
Our magnitudes and spectral classes were compared with the corresponding data ob- 
tained from the two overlapping fields, LF3* and LF3a;? no significant differences were 
observed. In the case of LF3a, for example, we had in common 215 stars. The average 
systematic difference (LF3a minus LF3b) in mp, for these stars was +0.05 + 0.01. 
The average color difference of the 115 stars between magnitudes 9 and 11 for which 
colors were available in both fields was less than 0.01 mag. 


4 Nassau and Burger, Ap. J., 103, 25, 1946. 
5 Nassau and Seyfert, Ap. J., 103, 117, 1946. 
6 Nassau and van Albada, Ap. J., 106, 20, 1947. 


7 The luminosity notation of Morgan, Keenan, and Kellman, Am Allas of Stellar Spectra (Chicago: 
University of Chicago Press, 1943), is used throughout this paper. 


8 Nassau and Harris, unpublished. 
® Annear, unpublished. 
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Only 81 of our stars, mostly bright, are in Schalén’s'® field in Cygnus. A comparison 
between his photographic magnitudes and ours showed a difference (Uppsala minus 
Case) of 0.08 + 0.02. Comparing the spectra, we found that stars which we had classi- 
fied as AO-AS were classified later by Schalén, in some cases, as late as F5. This might 
prove of considerable interest, since many of the stars in this region are metallic-line 
A stars. In the later classes, beginning with F8, there is good agreement bet ween the two 
systems of classification. The agreement in the luminosity classification in the two sys- 
tems is also satisfactory. 

Comparing the 93 stars between magnitude 9 and 11 common to the Henry Draper 
Extension“ and our own list, we found a systematic difference (HDE minus Case) of 
+0.C4 + 0.02 mag. As far as we were able to ascertain, there is no scale difference be- 
tween the two magnitude systems in this region. We found no systematic difference 
between the two systems of classification, with the exception of the spectra for early-type 
stars. There is a tendency for our BS stars to be classified as B8, B9, or even AO in the 
HDE. From G8 to K5 we classified earlier by one subdivision. The agreement in the 
range AO-G5 and among the M-type stars is good. 

Our field includes 50 stars in common with the Henry Draper Catalogue. With this 
very limited material it is difficult to draw significant conclusions, but it appears that the 
tendency observed in the HDE in the case of B5 stars is present in the HD as well. 


III. INTRINSIC COLORS 
EARLY-TYPE STARS 


Our present knowledge of the intrinsic colors of early-type stars appears to be limited. 
For their determination we drew upon three sources. 

Seares and Joyner” have obtained the intrinsic photovisual colors (C) of B3-M stars 
near the pole. The number of stars at the extremes of the range of spectral types is too 
small to give by themselves reliable mean intrinsic colors, but from AO to KO the number 
of stars in each spectral class is sufficient. Their work has established the intrinsic color 
index of an A5 star as zero, and we have adopted this result. Table 1, column 2, gives 
their values for the early-type stars, which must be strengthened, if possible, by data 
from additional sources. 

From King’s extra-focal photometry of 100 bright stars'* we obtained a second estimate 
of the colors. His values, reduced to zero at A5, are shown in Table 1, column 3. 

The photoelectric colors (C;) of stars brighter than visual magnitude 6.0 and more 
than 15° from the galactic equator given by Stebbins and Huffer'* were reduced to the 
international scale (C) by means of the ratio” 


—=1.90. 1) 
ra 1.90 ( 


1 


The zero point was adjusted so that the intrinsic color of a B8 star is equal to —0.29 
mag., the value obtained from the two previous determinations. It was felt that this 
procedure partially overcame the effect of the hydrogen absorption at AO. The results 
thus obtained are given in Table 1, column 4. Our adopted photovisual colors (col. 5) 
were obtained by a judicious averaging of the values given by these three determinations. 

In order to reduce these color indices to our photored system, we utilized the data ob- 
tained for 23 stars between magnitude 6.5 and 10.9 (mostly bright stars of spectral classes 


10 Uppsala Medd., No. 55, 1931. 

1 Harvard Ann., Vol. 100, No. 1, 1925. 

2 4 p. J., 98, 261, 1943. 

13 Harvard Ann., Vol. 85, No. 3, 1923. 14 Pub. Washburn Obs., Vol. 15, Part 5, 1934. 
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between AO and AS at the pole) for which photovisual and photored magnitudes were 
available.‘ The ratio of the photored color index (RJ) to the photovisual color (C) is 


as (2) 


The intrinsic red indices derived by the use of this ratio are given in Table 1, column 6, 
The same values are repeated for the sake of completeness in Table 3, second column. 


LATER SPECTRAL TYPES 


Intrinsic colors for the rest of the main-sequence stars as well as for the giants were 
obtained as follows. The red magnitudes of 202 stars within about 2° of the North Pole 
have been determined (unpublished). The photographic magnitudes of these stars and 
their color indices on the international scale have been given by Seares, Ross, and 
Joyner.'® For the spectra we used the classifications given by Nassau and Seyfert.® The 


TABLE 1 
INTRINSIC COLOR INDICES 


—0.47 2 —0.29} —0.29! —0.30} —0.40 
— AO....) —0.15).......| —0.19; —0.15} —0.20 


—0.39} —0.31} —0.36, —0.38} —0. 50 | 


mean ratio of RJ to C for all these stars was found to be 1.33. The ratio was also calcu- 
lated for all spectral groups having 10 or more stars each, and no systematic variation 
with spectral class was found. Since the coefficient of interstellar absorption in the direc- 
tion of the North Pole has been determined,* ” it was possible to confirm the validity 
of the ratio 1.33 in the case of intrinsic colors. With this ratio and the color indices ob- 
tained by Nassau and Seyfert,' the intrinsic red color indices of all stars B8 and later 
were computed. The resulting intrinsic red indices (RJ;) are given in Table 3, beginning 
with spectral class A2. The colors of GO III and G2 III stars were taken from the 
values of the corresponding dwarfs by adding 0.13 mag. 


IV. PHOTOGRAPHIC ABSORPTION FROM RED EXCESS 


The regions LF3, LF3a, and LF3b contain 31 stars (mostly BO) for which Stebbins, 
Huffer, and Whitford'* have determined photoelectric colors. In the same paper Steb- 
bins and his associates give normal colors for stars between BO and AO. Thus an estimate 
of the ratio of photored excess (E,) to the corresponding excess on the photoelectric 
scale (E;) could be made. It was found to be 3.12. We are indebted to Mr. Harris® and 
Mr. Annear® for the colors of the majority of the stars used in the determination. 

The ratio E,/E, was also obtained independently from the data of the six-color pho- 
tometry of Stebbins and Whitford" as follows. From their paper we have E; = 0.130, 


16 Pub. Carnegie Inst. Washington, No. 532, 1941. 
16 4p. J., 91, 20, 1940. 
17 1p. J., 98, 20, 1943. 
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Apg — Ai = 0.813, and A, — A; = 0.411. (The effective wave length of our red mag- 
nitudes is 6200 A.) Therefore, 
E, = Ay — A, = 0.402, 


E, _ 3.10. (3) 


To obtain the ratio of the total photographic absorption (Apg) to the red excess, we 
have adopted the mean value Ayg/£; = 8.2 given by van Rhijn,'* which was the result 
of an exhaustive discussion of the observational material of many investigators. Com- 
bining this with equation (3), we have 


(4) 


Stebbins and Whitford!” showed that the range of A,,/£, is from 7.1 to 10.9. Therefore, 
the range of A,,/E, is from 2.3 to 3.5. In many of the regions of relatively high obscura- 
tion in which investigations have been carried out at this observatory, values of this 
ratio higher than that adopted yield unreasonable values of the absorption. 


V. ABSOLUTE MAGNITUDES 


EARLY-TYPE STARS 


The problem of assigning absolute magnitudes to early-type stars presents many well- 
known difficulties. From relatively recent determinations we have chosen three which 
are semi-independent in observationa] material and method of treatment. Nordstrém!* 
has collected and discussed the observational material up to 1936. His assigned values 
are given in the second column of Table 2. 

The absolute magnitudes given by Stebbins, Huffer, and Whitford" are the values ob- 
tained by Plaskett and Pearce and by Strémberg, derived from proper motions of the 
General Catalogue corrected for absorption. In addition to these, they included new 
values derived from radial velocities and from close pairs of B stars having the same color 
excess. Their results are given in the third column. Wilson,”° basing his determinations on 
the radial velocities of 987 stars and on the proper motions of 1496 stars, obtained abso- 
lute magnitudes of early-type stars which are given in the fourth column. The data of 
Table 2 were plotted, and a smooth curve was drawn. From the curve we read the 
adopted photovisual absolute magnitudes which are presented in the last column. The 
corresponding photographic magnitudes, given in the third column of Table 3, were ob- 
tained by applying the color indices given in Table 1, column 5. 


STARS OF SPECTRAL CLASSES A2-K5 


Our source for the absolute magnitudes of stars later than AO was the work of Russell 
and Moore.” Their values were reduced to photographic magnitudes by making use of 
the color indices given by Nassau and Seyfert.5 The resulting absolute magnitudes are 
collected in the third column of Table 3. 


ABSOLUTE PHOTOGRAPHIC MAGNITUDES PER UNIT VOLUME OF SPACE 
Table 4 indicates the spectral groups into which the stars were combined for statistical 
consideration. The dispersion in absolute magnitude for each individual spectral class 
18 Croningen Pub., No. 51, 1946. 


19 Medd. fran Lunds Astr. Obs., Ser. I1, No. 79, 1936. 
20 Ap. J., 94, 12, 1941. 21 4p. J., 87, 389, 1938. 
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was taken as +0.7 mag. in nearly all cases. This is in line with the values given by Rus- 
sell and Moore” but is somewhat greater than the values given by Wernberg.”* For the 
broader spectral groups proper adjustment of the dispersion was made. These values of ¢ 
are given in Table 4, third column for dwarfs and fifth column for giants. 


TABLE 2 
PHOTOVISUAL ABSOLUTE MAGNITUDES OF EARLY-TYPE STARS 


Nord- Nord- 
Stebbins Adopted “en 


—0.3 


Stebbins 


* Average for BO and BI. 
TABLE 3 
INTRINSIC RED INDICES AND PHOTOGRAPHIC ABSOLUTE MAGNITUDES 


Main SEQUENCE MAIN SEQUENCE GIANTS 


SPECTRUM SPECTRUM 


+4+++4+4+44 
SOS 
SReSLESES 


++++4++4+4+ 
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+++ 


TABLE 4 
ABSOLUTE PHOTOGRAPHIC MAGNITUDES PER UNIT VOLUME OF SPACE 


MAIN SEQUENCE GIANTS MAIN SEQUENCE GIANTS 


SPECTRUM SPecTRUM 


Me 


| 

oo 
0. +2.0} 0. 
+2.0 
| 


+2.0 


2 Ap. J., 92, 354, 1940. 23 Ann. Uppsala Obs., Vol. 1, No. 4, 1941. 
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To reduce the absolute magnitudes of Table 3 to values appropriate to a unit volume 
of space (Mo), we have applied the Malmquist* corrections. The necessary values of 
d log A(m)/dm were obtained from Seydl’s*® data. The values of Mo, rounded off to the 
nearest 0.5 mag., are shown in the second and fourth columns of Table 4. 


VI. INTERSTELLAR ABSORPTION IN THE REGION 


Guided by the general appearance of our direct photographs of the region and by the 
prints of the Ross Aflas, we divided our field into four areas within each of which the 
obscuration appeared relatively uniform. These areas are shown in Figure 2. The well- 
known parabolic obscuration is within the boundaries of area D. In this area the color 
excesses show a marked difference from those of the remaining areas. It is excluded for 
the present and will be discussed separately. The other three areas were treated together 
because no significant differences in the average color excesses were apparent. 


| | | | | | | 
E, 


Fic. 3.—Relation between observed color excesses and apparent modulus 


Figure 3 shows the relation between the observed color excess and the apparent modu- 
lus derived from 683 main-sequence stars (B2—KQ0) in the three areas. The color excesses 
for stars from B2 to A5 are represented by open circles, while those for the stars from FO 
to K5 are indicated by circles and bars. The radius of the circle in each case is propor- 
tional to the square root of the number of stars employed in forming the mean E,. The 
grouping of stars was made in half-magnitude intervals in apparent modulus. 

A curve drawn through the points in Figure 3 leads directly to the change of absorp- 
tion with distance. However, small irregularities in such a curve may produce spurious 
irregularities or clouds in the absorption. It is therefore appropriate to inquire how close- 
ly a linear relation between color excess and distance would fit the observed data. 


From equation (4) we have ; 
Ay = Kr=2.6E,, 


where K is the absorption coefficient in magnitudes per parsec and r the distance in 
parsecs. From 


M = m— 2.6E,+5—5 log r 


and 
log K = log (2.6E,) —log r, 


24 Lund Medd., Ser. 11, No. 22, 1920. *> Pub. Prague, No. 6, 1929. 
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we obtain 


where 


(7) 
v=5 log E,+ 2.6E, . (8) 


Each point in Figure 3 leads to a pair of values of uw and v and hence to a determination 
of B. The average value of B yields a final mean value of the absorption coefficient K. 
In forming the mean value of B, attention must be paid to the weights to be assigned 
to each individual value corresponding to each point in Figure 3. The weight should be 
proportional to the number of stars, m; but, in addition, it must provide for the fact that 
the same error in #, produces widely different errors in B at different parts of the curve. 


TABLE 5 
CALCULATION OF THE LINEAR ABSORPTION COEFFICIENT 


e 


| 
| 
| 


Table 5 shows the color-excess data for the B2—AS stars and the calculated values of B 
along with the corresponding weights (wz). In calculating wg we have used equations 
(6), (7), and (8), assuming that the mean errors of E, and u for a single star were 0.1 and 
1.0 mag., respectively. The weighted mean value of B is 13.50. The corresponding value 
for the FO-K5 dwarfs is 13.31, but the weight of this determination is low on account of 
the nearness of these stars to the sun. The adopted value of B based on all 683 stars is 
13.47 + 0.11, which leads to an absorption coefficient of 0.53 mag. per kiloparsec. The 
accidental error of the linear absorption coefficient is only a few hundredths of a mag- 
nitude. 

The curve in Figure 3, as well as the straight line in Figure 4, was drawn on the basis 
of this value of the linear coefficient. The fit appears to be satisfactory. 

The giant stars from FO to K3 were not included in the above discussion. Instead, 
their intrinsic colors were obtained from their observed colors and the adopted interstellar 
absorption in the region. The mean difference between these and the corresponding in- 
trinsic colors of the dwarfs is +0.05 + 0.02 mag., a value less than expected (Table 3); 
but, considering the limited material available (76 stars), we may conclude that the 
adopted absorption is not seriously in error. 


VII. SPECTRAL DATA 


The spectrum and magnitude data for the stars in the regions A, B, and C, which com- 
prise 9.9 square degrees, are shown in Tables 6a and 6b. The tabulation of photographic 
magnitudes is in intervals of half a magnitude. For stars brighter than 11.0 mag. we 
recognize three luminosity classes, namely, ITI, 1V, and V. The region had practically 


47 
B=u-?, (6) 
— 
; 
7.5.......| 0.08 | —5.1 | 12.6} 10 | 1 | 11.3......| 0.28 | 13.8} 34 | 13 
8.8.......| 14.1] 33 | 6 || 12.3... 9 13.3] 34 | 13 | 
ve 8 13.2}; 20 | 8 
9.8.......| .16| —3.6| 13.4] 65 | 20 || 13.3.. 1 12.8] 26 11 
10.3.......| .19] —3.1] 13.4] 83 | 30 |] 13.7.. 1 13.2} 23 | 10 
10.8.......| 0.21 | 13.6} 55 | 20 || 14.4.. 5 | | 13.6| 4 | 2 
r 
| 


3.6 KPS 


Fic. 4.—Total photographic absorption in regions A, B, and C. The dashed lines show the absorption 
in region D. 
TABLE 6a 
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TABLE 66 
REGIONS A, B, AND C 
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no supergiants, For stars fainter than 11.0 mag. we give two classes only; where we were 
"able to designate stars as class IV we included them with the dwarfs for stars up to G2 
‘ and with the giants for stars of class G5 and later. We were unable to assign luminosity 
classes to a number of stars, primarily because of overlapping spectra. Such stars are 
included in the column marked “Un.” 

For region D, which comprises 0.96 square degrees, the corresponding data are given 
in Table 7. FO-G2 stars classed as IV were grouped with the dwarfs, while for G5 IV and 
later this class was grouped with the giants. 


VIII. PERCENTAGES OF DWARFS 


The percentage of dwarfs for various spectral and apparent magnitude groups as ob- 
tained from the observational data of LF3b are listed in the left-hand part of Table 8. 
These data refer primarily to the faint stars and are subject to accidental fluctuations, 


TABLE 7 
REGION D 


3 
3 


24 


especially in the G8-K3 group, where the number of stars is small. Some procedure must 
be found for smoothing the figures and, more important, for extrapolating to the 
brighter stars for which little data are available. 

It is clear that the decrease in percentage of dwarfs with advancing spectral type for a 
single magnitude interval depends primarily on the mean absolute magnitudes of the 
various spectral groups and to some extent on the relative numbers of giants and dwarfs 
per unit volume of space. The second effect will be negligible if the absolute ratio of 
giants to dwarfs is constant for all spectral classes, but this cannot be assumed, and hence 
it is impossible to predict the change of percentage of dwarfs with spectral group. Our 
density analysis suggests, however, that the ratio is indeed nearly constant. 

It is also apparent that the increase in percentage of dwarfs with increasing apparent ~ 
magnitude within a spectral group is partly due to the effect of whatever absorption is 
present. The more distant giants will be obscured, while the fainter dwarfs will be un- 
affected, and their numbers will increase with apparent magnitude more rapidly than will 
the giants. In addition, variation of the space density of the giants or dwarfs with dis- 
tance will affect the percentages of dwarfs for various apparent magnitudes. 

Tt has been clearly brought out by Bok and Rendall-Arons,” however, that, if the 
numbers of dwarfs are well determined at the tenth to eleventh magnitude, it is possible 


% 1p. J., 101, 280, 1945. 
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to predict the percentages of dwarfs among the brighter stars on a very plausible assump- 
tion. The F8-G2 dwarfs brighter than m,, = 11.5 will be closer than 200 parsecs, and 
the later-type dwarfs brighter than the same limit will be still closer. It is reasonable, 
therefore, to assume in the present state of our knowledge that the density of these 
dwarfs is constant down to these magnitudes. In our case the number of dwarfs of 10-11 
mag. is well determined; hence the number of dwarfs in the brighter magnitude inter- 
vals for each spectral group can be calculated. Thus we were able to predict the percent- 
ages of giants and dwarfs and to carry out the density analysis for the giants. It was as- 
sumed that the density for the dwarf F8-G2, G5, and G8-K3 stars was uniform over the 
first 200 parsecs. The resulting percentages of dwarfs in the three spectral groups for the 
brighter stars and the smooth values for the fainter stars as given in the right-hand side 
of Table 8 were calculated on this basis. The FO-F5 dwarfs extend over too great a dis- 
tance to make the same assumption valid; the percentages of dwarfs for this group were 
assumed to remain constant at 83, as observed in the 10—12-mag. range. 


TABLE 8 


PERCENTAGE OF DWARFS 
Mog 


SPECTRA Observed Data Smoothed Values 


When the density analysis for the giants was carried out, after applying the percent- 
ages of the right-hand side of Table 8, it was found that the counts could be satisfied by 
a constant density. It is worth emphasizing that this constant density for the giants is 
not in itself an assumption. The density for the giants at 700 parsecs is quite independent 
of any assumption about the density of the dwarfs. The rapid decrease with advancing 
spectral type in the ratio of the number of dwarf to giant stars shown by the data in 
Table 8 has been noted in previous investigations. In the case of the stars near the North 
Pole the values obtained by Nassau and Seyfert® show this decrease, but it is a some- 
what more gradual one, as would be expected from the fact that this region has a galactic 
latitude of 28°. 


IX. ADOPTED VALUES OF LOG A’(m) 


In order to obtain the value of log A’(m)”" for the bright stars, we have used the data 
of Seydl** and the results of a separate analysis from the data of the Henry Draper Cata- 
logue. In the latter case the spectrum data for an area of 100 square degrees was used in 
the region of our field. The magnitudes used in the Seydl counts were reduced to photo- 
graphic. The separation into giants and dwarfs within each spectral group was made by 
means of the data of Table 8. 

The data for the fainter stars were taken directly from Tables 6a and 6b. The un- 
classified stars given in these tables were divided into giants and dwarfs in accordance 


27 4'(m) is the number of stars per 100 square degrees of the sky between m — } and m + }, where m 
is the apparent photographic magnitude. 
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with Table 8. The values of log A’(m) thus obtained were plotted against magnitude, and 
smooth curves were drawn through the points. Except for the B2 and BS groups, no 
marked discontinuities were observed where the data of the various sources overlapped. 
Table 9 presents the smoothed values of log A’(m) which were used in the succeeding 
analysis. 
X. CALCULATION OF DENSITIES 

The calculation of the space densities from the observed numbers of stars is a classical 
problem, and many procedures are available. We have used one described recently by 
Malmquist”* in which the numbers of stars in the various distance shells are read directly 
from the log A’(m)-curve. The effect of the small dispersion in absolute magnitude for 
stars of a given spectral group is allowed for in this method by entering the log A’(m)- 


TABLE 9 
ADOPTED LOG A’(m) IN LF3b 
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curve with an adjusted valve of m. The effect of the absorption was determined by the 
method of Seeliger.** 
If the variable mp is defined by the equation 


my = My—5+5 log 1, 


where M, is the mean absolute magnitude of the stars considered, per unit volume of 
space, Malmquist has shown that the number of stars at the distance ro corresponding 
to mo is given by 


o 
n (mo) A’(m). (9) 


In this formula 
_ 


m= 
\ oot om 


28 Uppsala Medd., No. 84, 1942. 
29 See Bok, The Distribution of the Stars in Space (Chicago: University of Chicago Press, 1937), p. 14. 
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The quantity ¢,, is the familiar Malmquist correction* to change Mo into Mn, and om is 
the corresponding corrected value of oo, the dispersion per unit volume of space, viz., 


én —My= 
dm 


2 [1 d?in A’ (m) ] 


2 
dm 


The assumptions under which equation (9) is valid are that the frequency function of 
the absolute magnitudes in an element of space is of the Gaussian form and independent 
of the distances of the stars and, further, that the log A’(m)-curve can be represented by 
an expression of the first or second degree in m. We have found that the second assump- 
tion does not impose a serious restriction on the analysis in this field. 

Malmquist’s formula has been modified slightly for greater convenience in the anal- 
ysis. When the second derivative of the log A’(m)-curve is small, terms in y* can be 
neglected, and equation (9) takes the form 


log n (mo) = — Mod. +log A" (mo+ 5%), (10) 
where 
log A’(m) _ 
2 Mod. dm? 


_ __o__dlog A’(m) 


If the log A’(m)-curve is a straight line, 1 = 0, and the formula is further simplified. 
The quantities u and ¢,, are to be evaluated at the point 


Em 


They are most conveniently determined by numerical differentiation of the values read 
from the smoothed log A’(m)-curve. Next a graph relating m to mp is constructed, and 
the right-hand side of equation (10) is computed for values of mo differing by five-tenths 
over the range of the counts. The effect of the increasing volume of the successive shells 
is allowed for, and the resultant densities are ready to be corrected for the effect of ab- 
sorption. 

As a test of the method, an independent analysis of the B8—A0 stars was carried out by 
means of an (m, log 7) table. No significant differences could be found in the two anal- 
yses. The Malmquist method has much to recommend it in its directness and its defini- 
tiveness. The (m, log x) table retains the advantage of showing clearly at all times how 
the stars of each magnitude group are distributed in space. 


m= Mo+ 


XI. RESULTS OF DENSITY ANALYSIS 


The results of the density analysis are given in Table 10 and Figure 5. In the latter the 
ordinates represent the logarithms of the number of stars per 1000 cubic parsecs; the 
abscissae are distances in parsecs. The dotted portions of the curves and the figures in 
parentheses in the table represent the less certain determinations. 

The densities for the early-type stars all show a decrease with distance from the sun. 
It is most marked in the case of the B8-A0 stars. The log A’(m)-curves for the B2 and B5 


, 
m) 
| 
i 
j 
and 
€ 
: 
i 
te 


TABLE 10 
SUMMARY OF DENSITIES, STARS PER 1000 CUBIC PARSECS 


MAIN SEQUENCE 
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— 68-K3 Z 
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~-- FO°-FSZ 


7.0 


600 800 1200 1400 1600 1800 2000 


F1G. 5.—Logarithm of densities per 1000 cubic parsecs plotted as a function of the distance in par- 
secs. 
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stars were less satisfactory, but there is nevertheless no doubt of the reality of a negative 
density gradient for these stars. Some of this effect can be attributed to the variation of 
density with distance from the galactic plane, since our region is at galactic lati- 
tude +3°2. 

Tt is appropriate to call attention to the uncertainties in the density analysis due to 
the lack of separation of BO, B2, and BS stars at the faint limit of the Henry Draper 
Catalogue. Stars designated as “‘B”’ in this catalogue were arbitrarily divided equally be- 
tween B2 and BS. Furthermore, the Seydl counts for B2 stars are averages for a zone 10° 
on either side of the galactic plane. Because these stars are so closely concentrated to this 
plane, the average number per square degree will not be representative except for the 
very brightest stars. This is also true for the counts obtained from the 100 square de- 
grees around our field that were taken directly from the HD. The counts of B5 stars will 
be affected in the same way but to a lesser extent. Finally, there is a certain amount of 
arbitrariness in combining our counts of B stars with the counts for the brighter stars 
averaged over all longitudes in view of the well-known clustering of these stars. 

The constant densities of the giant stars would have to be modified if different per- 
centages of giants and dwarfs are found to exist among the apparently bright stars. The 
densities of the giants at 600-800 parsecs and the densities of the dwarfs at 100-200 par- 
secs are based on observed numbers of these stars. The absorption has virtually no effect 
on the derived densities of the dwarfs. 

There is a gradual progression in the densities from one end of the spectral sequence 
to the other which will be recognized as reflecting the general character of the luminosity 
function for the stars at large. The least luminous stars outnumber the brightest by a 
factor of the order of 3000, although the latter seem to be particularly numerous in this 
field. In LF3b there are about five times as many B2 and BS stars as in LF1 or LF2.*° 

The density of the giants of classes FO-K3 is about the same as the density of the A2- 


AS stars, which have approximately the same luminosity. The latter group includes a 
number of metallic-line A stars which perhaps should have been treated as F stars in the 


analysis. 


XII. THE LUMINOSITY FUNCTION 


The luminosity function ¢(M) is defined by the statement that ¢(7)AM is the number 
of stars per cubic parsec between absolute photographic magnitude M — 34M and M 
+ 3AM. The data of Table 10 and Figure 5, along with the mean absolute magnitude 
(Mo) and dispersions of the various spectral groups (Table 4), provide us with the means 
of calculating ¢(M) at various distances. 

The numbers of stars per cubic parsec for each spectral group at the chosen distance 
were spread out over a range of absolute magnitude according to the appropriate values 
of o and M». The giants were treated as one group. The numbers for each absolute mag- 
nitude in intervals of 1 mag. were summed over all spectral groups and tabulated as 
(M) at that distance. Four such functions were obtained at distances of 100, 200, 400, 
and 600 parsecs. The formulation of the luminosity function at a distance of 100 parsecs 
is illustrated in Figure 6, in order to show the roles played by the absolute magnitudes 
and dispersions. A change in the mean absolute magnitude for a spectral group will have 
a marked effect on the luminosity function because both the density of this group and 
the absolute-magnitude range to which it contributes will be altered. This effect is dis- 
cussed in another paper of this series.*! 

The results of the calculations of ¢(M) at the four chosen distances are given in Fig- 
ure 7. The solid curves represent the well-determined portions, while the dashed por- 
tions refer to regions where less certain densities were used. These uncertainties are due to 


3° McCuskey and Seyfert, Ap. J., 106, 1, 1947. 
31 McCuskev, 4p. J., 109, 414, 1949. 
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the unknown numbers of dwarfs at distances of 400 and 600 parsecs and to the fact that 
the number of highly luminous stars at 100 and 200 parsecs is poorly determined. No 
serious extrapolation of the density-curves was made in any of the solutions for ¢(M). 


8 10 M 


Fic. 6.—Formation of the luminosity function (at 100 parsecs). The abscissa gives the absolute 
magnitude, the ordinate the values of log } @ (47). Curves: /, B8—AO stars; 2, all giants; 3, A2-A5; 4, 
F0-F5; 5, F8-G2; 6, G5; 7, G8-K3. The luminosity function is shown by the dotted curve. 


eit? T T 
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Fic. 7.—Luminosity functions at 100, 200, 400, and 600 parsecs. The ordinate log ¢ (M) for the curve 
at 100 parsecs is shown. The scales for the other distances are displaced downward by 1, 2, and 3 units. 
The abscissa gives the absolute magnitude. 


A comparison of our results with van Rhijn’s luminosity function® (dots and dashes, 
Fig. 7) shows a general agreement over nearly all the range of M. The two determina- 
tions differ in shape, but the curves crisscross one another in such a way as to leave the 
total number of stars per cubic parsec roughly the same. Between absolute magnitudes 


382 Groningen Pub., No. 47, p. 17, 1936. 
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0 and —2 (BS stars) we find an excess in our field at 200 and 400 parsecs amounting to 
three- or four-tenths in the logarithm; but this is reduced at 600 parsecs, where the de- 
creasing density of the early-type stars comes into play. 


XIII. REGION D 


COLORS 


As has been stated, the stars in region D show a marked difference in color excess 
when compared with the stars in the other regions. From the general appearance of direct 
photographs, an interstellar obscuration appears to extend over a number of degrees in 
a direction parallel to the galactic equator. Its width varies from about 1° to 4° or more. 
Our present study is confined to the northernmost portion and covers an area of 0.96 
square degrees. Table 7 presents the spectral data for this region. In all, we are dealing 
with only 138 stars, and a number of these are bound to be foreground stars with re- 
spect to the obscuration. 

In order to examine for possible field errors in the determination of the colors of these 
stars, we computed the average red index of 16 main-sequence F5 and F$8 stars brighter 
than 11.0 photographic magnitude in region D and in the adjacent region C. The aver- 
age was obtained by first reducing the colors of F8 stars to that of FS. The results were: 
RI (region D) = 0.30 mag., and RI (region C) = 0.25 mag. No serious field error in the 
colors is indicated because the intrinsic red index of an FS star is 0.28 mag. 

The open circles of Figure 4 show the results of the color-excess analysis in this region. 
The two broken lines drawn approximately through the points suggest that the cloud sets 
in at about 900 parsecs and absorbs about 0.9 mag. 


STAR COUNTS 


Another estimate of the absorption is provided by the general star counts. Four areas 
were chosen, each three-tenths of a square degree in area and one in each of the four re- 
gions. The areas in regions A and B were at galactic latitude +3°0, and the areas in re- 
gions C and D at +1°7. The blue-magnitude sequence for the counts was determined 
by Nassau and van Albada.* The limiting magnitude of the stars counted was 17.0. No 
significant differences were observed in the values of log A(m)** for regions A, B, and C; 
the average curve for the three regions is shown by full circles in Figure 8, where the 
ordinate is log A(m) per square degree and the abscissa is the apparent photographic 
magnitude. We found no discontinuity between the values obtained from our spectral 
survey and the star counts. The corresponding results for region D are shown by open 
circles. For stars fainter than 12.5 mag. we observed a deficiency of stars in region D, 
obviously due to the cloud. The amount of absorption as compared with the other regions 
slowly increases, reaching the value of 0.7 at m,, = 16.0. This result is of the right order 
of magnitude when compared with that obtained from the colors. 

The curve of log A(m) for region D (Fig. 8) shows an interesting excess of stars be- 
tween mpg = 9.5 and mpg = 11.5, which seems to be real. Table 11 gives the numbers of 
stars per square degree for different spectral groups between m,, = 9.5 and m,, = 11.5. 
The second column is for the regions A, B, and C combined; the third column is for a 
region bordering region D; and the fourth column is for region D. The percentage excess 
of stars in region D when compared with the three regions combined is given in the last 
column. It appears that the early-type stars are those which contribute to the excess. In 
our studies of LF3 and LF3a, which are adjacent areas and at lower latitude, we found 
an excess of early-type stars over other areas in the Milky Way. 

In spite of the limited material available, an analysis of the stellar density in this re- 
gion was made, first, from our spectral data and, second, from the general star counts. 


33 4(m) is the number of stars per 100 square degrees of the sky between m —} and m + }. 
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ANALYSIS OF BS~AS STARS IN REGION D 

The analysis was carried out by means of an (m, log w) table. The sudden increase 
in the absorption and the fact that the distribution of M is not Gaussian make this meth- 
od of analysis more suitable. The luminosity function of the BS—AS stars at r = 400 par- 


12 13 14 


Fic. 8.—Star counts for regions A, B, and C are shown as full circles and for region D as open circles 


TABLE 11 
NUMBER OF STARS PER SQUARE DEGREE BETWEEN mpg=9.5 AND mpg=11.5 


Excess 
Regions | Region Region (Per 


Spectrum 
A, B, C} near D Cent) 


140 
143 
96 
54 


Regions | Region Region 
A, B, C} near D 


secs was calculated from the results of our density analysis in the clear region. It has a 
maximum at Mo = +2 and is asymmetric with a larger number of stars on the side of 
smaller M. With this luminosity function an (m, log x) table was constructed, and the 
general absorption of 0.53 mag. per kiloparsec was introduced. The analysis revealed the 
gradual decrease in density of the B5—AS stars from 400 parsecs in agreement with what 
had been found before. 

A cloud of 0.9-mag. thickness was then introduced in front of the 1000-parsec shell, 
and the density was adjusted to fit the counts of the B5—AS stars in region D. It was 
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found necessary to increase the density in the 630- and 1000-parsec sl.ells to three times 
that previously obtained in order to provide a satisfactory representation. 


ANALYSIS OF GENERAL STAR COUNTS IN REGION D 


The basis of the analysis for the general star-count material was a luminosity function 
taken from the work of Bok and MacRae* at galactic latitude 2°. As before, the general 
absorption of 0.53 mag. per kiloparsec was applied to the table. For the three regions the 
observed star counts agreed with the counts predicted from the (m, log m) table, indi- 
cating a uniform density up to 2500 parsecs. For region D an additional absorption of 
0.9 mag. was introduced at a distance of 900 parsecs. An increase in the density by a 
factor of 2-3 at distances between 600 and 1000 parsecs satisfied our counts and is in 
agreement with our previous result. The stars which contribute mostly to the greater 
numbers of stars between 9 and 12 mag. are those between absolute magnitudes —1 
and +3. Miller® observed a similar excess in the number of stars in discussing star counts 
near our region. 

We may add that an examination of the individual colors of the stars in region D in- 
dicated that the more luminous stars contribute the larger color excesses, while the less 
luminous stars seem to have nearly normal colors when the general absorption is taken 
into account. 

We leave for future investigation, based on the combined data for all three LF re- 
gions, the question as to whether this density increase is closely associated with the 
parabolic obscuration which lies in region D or is a general characteristic of very low 
latitudes in this part of the sky. 


34 Ann. New York Acad. Sci., 42, 219, 1941. 
% Harvard Ann., 105, 297, 1937. 
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ABSTRACT 

Bright hydrogen lines were discovered at Mount Wilson in 1921. A few slit spectrograms at several 
observatories showed a considerable difference between the displacements of dark hydrogen lines and 
those of other elements. 

The present study with higher dispersion shows that the star is a spectroscopic binary with a period 
of 28 days and a velocity range of 209 km/sec. There are two sets of dark hydrogen lines, both of variable 
intensity, one of which usually yields the same velocity as the lines of other elements, while the other 
has a nearly constant negative displacement. 

The facts suggest that the binary orbit is enveloped by a cloud of incandescent hydrogen whose outer 
absorbing layers are rapidly expanding. The system resembles 8 Lyrae in certain respects. 


Bright hydrogen lines were discovered in the spectrum of HD 187399! from an objec- 
tive-prism photograph of the Ha region? taken at Mount Wilson on September 6, 1921. 
The first slit spectrogram was probably that taken at Mount Wilson on April 9, 1924, 
with a one-prism spectrograph attached to the 60-inch telescope. On this plate, H8 had a 
bright component, while the other hydrogen lines were wholly dark. The dark hydrogen 
lines yielded a radial velocity of —66 km/sec, lines of other elements —13 km/sec. 
Swings and Struve* found an even greater difference in the displacements on spectro- 
grams taken at the Lick and McDonald observatories in 1942. Their values were: 
H, —65 km/sec; other elements, +62 to +18 km/sec. The anomalous negative dis- 
placement of the hydrogen lines suggested an expanding shell, although other atoms 
seem not to take part.in the expansion. Moreover, the lines of Mg ti and Si m were 
stronger than in the typical shell spectrum. The star obviously deserved further ob- 
servation. 

The present account is based on grating spectrograms, dispersion 10 A/mm (20 A/mm 
in the red), taken with the 100-inch telescope on Mount Wilson during the years 1944- 
1948 (see Table 1). The emulsion of most of the plates was Eastman IIa-O. 

The measurable lines are due to H, He 1, and ionized metals (M 11}, chiefly Mg 1, 
Siu, Cau, Ti u, and Fe 11. The H lines exhibit great variations in intensity and struc- 
ture. The other stellar lines, although rather weak, are fairly narrow and are measurable 
with reasonable accuracy; with the exception of Ca 11 lines, which are fairly strong on 
some plates and scarcely visible on others, they exhibit little variation in appearance. 
Their displacements, however, differ greatly from plate to plate. The Ca 1 lines appear 
considerably stronger at the extreme velocities, especially toward the minimum velocity, 
than at the intermediate velocities when they are superposed on the interstellar lines. A 
study of the plates and of several tracings yields the impression that the apparent change 
in intensity is too great to be due wholly to contrast or Eberhard effect. There may be a 
faint Ca 11 component having about the same displacement, —95 km/sec, as that of the 
anomalous, nearly fixed hydrogen component B (see section on “‘Hydrogen Lines’’). This 


1 MWC 321; BD+29°3754; 1900 R.A. 19544"7, Dec. +29°10’; mag. 7.7; spect. AO (Harvard), B9e8 
(Mt. Wilson). 
2 Merrill, Humason, and Burwell, Mit. W. Contr., No. 294; Ap. J., 61, 389, 1925. 
3 Ap. J., 97, 194, 1943. 
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component may be responsible for some apparent enhancement of the shifting Ca 1 
lines when they have large negative velocities. The interstellar Ca 11 lines are sharp; they 
appear somewhat stronger than normal when the stellar lines are superposed on them. 


TABLE 1 


RADIAL VELOCITIES 
(Km/Sec) 


VELocity AND NuMBER oF LINES 


Ca 1 int. 


1944 Aug. 2 .75 | (—101) 1, (+25) 
Sept. 6 ; (—120) 1 
Sept. 8 : — 90 15 

1946 July 13 : — 96 16 
Aug. 15 : — 54 16 
Sept. 14 : — 85 15, +12 8° 

1947 June 30 ; (— 94) 6, +397 
Aug. 3 
Aug. 5 
Aug. 25 
Aug. 26 
Sept. 27 
Oct. 6 
Oct. 8 
Oct. 25 

1948 June 16 
June 17 
Aug. 15 
Aug. 16 
Sept. 11 | 


* See Fig. 2. 
RADIAL VELOCITIES 


All stellar lines except certain components of the hydrogen lines appear to yield the 
same velocities. For example, lines of individual elements on plate Ce 4732 yield the 
results in Table 2. 

TABLE 2 
RADIAL VELOCITIES MEASURED ON PLATE Ce 4732 
(Km/Sec) 
Velocity and 
No. of Lines 


+38.7 (7), —94 (6) ....... 
+40.4 (5) 
+38.0 (1) +39.1 (15) 
+40.1 (2), —9.4 (2)* +39.1 (29) 


Element Velocity and No. of Lines Element 


* Interstellar. 


The radial velocities of He 1 and M 11 (ionized metals) in Table 1 are well represented 
by an elliptic orbit (see Fig. 1). The orbital elements computed by Miss Louise Lowen are 
in Table 3. This orbit represents the observations of Swings and Struve’ in 1942 within the 
errors of measurement. To represent my single low-weight plate of 1924, a correction to 
the period of about +0.005 day would be required. 


Be 
CE Dare 243 | 
| H | Hex a |) 
— 8 2| — 86 23 | -10.4 SE. 
4321.......| | -132 3 | —132 24 |:-11.4 
— 39 3| — 47 19} -10.2 en 
+ 24 2/419 15| — 9.6 be 
+ 39 2} +39 —9.4 
4777.......] (— 6) 2) — 13 21] — 9.6 
4800.......| (— 44.1] -—57 5] -11.4 
4809.......| + 38 3/ +34 11| — 9.8 
4874....... 6 12] -11.5 
4919... (— 71) 1| — 66 22| -10.7 
5201.......| |- 46 1|-48 3|-81NaI 
| — 24 4] — 21 —12.6 
|} +28 3) +28 8| -11.4 
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HYDROGEN LINES 


There appear to be two independent sets of dark hydrogen lines, both of variable 
intensity. Set A yields the same velocity as the lines of M 11 and He 1; set B has a nearly 
constant displacement not far from —95 km/sec. At Ha and Hf, B is regularly domi- 
nant, but at lines of shorter wave length this is not always so. Set A,in fact, seems to have 
the slower decrement on plates where it is well separated from B. The apparently slower 
decrement may be due to superposed hydrogen emission, which tends to weaken the 
longward lines of A but not those of B. On this hypothesis the lines of B would be formed 
in an outlying absorbing layer that envelops the whole system, while the lines of A arise 
in a stellar atmosphere along with those of M m and Het. 


+80 T T T T T T 


SEC 


+40 


2432450 J.0. 2460 2470 2480 


Fic. 1.—Velocity-curve of HD 187399 from measures of metallic lines 


TABLE 3 


ORBITAL ELEMENTS OF HD 187399 


.97 days K = 104.5 km/sec 
2432465 .98 vy = —18.9 km/sec 
10°7 a sin i = 37,600,000 km 
J(m) =2.72 


Mean displacements of one or both sets are recorded for most plates in Table 1. On 
other plates the results are complicated by a relative displacement from line to line 
along the Balmer series. Several of the more striking cases are plotted in Figure 2. It 
seems probable that some of the complicated details have their explanation in the various 
degrees of blending of sets A and B and of the varying influence of hydrogen emission. 

The bright 8 line is wide, poorly defined, and mutilated by absorption components; 
hence its center is difficult to measure. The measured velocities on ten plates range from 
+33 to —47 km/sec, with a mean of +2 km/sec. The overexposed Ha line on one plate 
yielded a value of +7 km/sec. Because the absorption components lie chiefly in the 
shortward half of the line, the measurements of the emission may have a systematic 
positive error. Hence it is quite possible that the true center of the emission corresponds 
in position to the velocity of the center of mass of the binary system, y = —19 km/sec. 

The facts suggest that the binary orbit is enveloped by a cloud of incandescent hydro- 
gen whose outer absorbing layers are expanding at the rate of about 75 km/sec. The 
system seems to resemble that of 6 Lyrae in certain respects. 
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In a number of Be stars, the bright hydrogen lines do not oscillate with the dark lines 
and thus may originate in a cloud partly independent of the binary stellar system. This 
circumstance could probably best be studied by means of special observations of the 


strong /Ha line. 
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Fic, 2.—Radial velocities of two sets of dark hydrogen lines on various plates. The velocity yielded 
by the stellar lines of ionized metals is indicated by the short dashed lines at the edge of each diagram. 
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SPECTROSCOPIC OBSERVATIONS OF TW DRACONIS* 


BURKE SMITH 
McDonald Observatory 
Received March 29, 1949 


ABSTRACT 


New circular elements have been derived. The epoch of total eclipse as determined from the spectro- 
graphic material is 0.146 day later than that computed from Kordylewski’s light-minimum. The rotation 
effect is found to be small. A few lines of the secondary were observed. 


The variability of TW Draconis, a = 15"32™4, 6 = +64°14’ (1900), Sp. A6 + K2, was 
discovered in 1910 by Miss Cannon.! The normal magnitude is 7.8, but during principal 
eclipse, which is total, the light diminishes to 8.9 mag. The light-elements were given by 
Graff in 1913 as 

Primary min. = JD 2418906.453 + 2.80654E. 


A more recent determination by Kordylewski* gives 
Primary min. = JD 2423711.321+4 2.80664E. 


According to Gaposchkin,* the photometric elements are 
A,=2.55,  As=0.08, L,=0.89, 1=0.21, k=0.65, 


i=85°6, Ratio of surface brightnesses = 19.9 . 


During June and July, 1947, fifty-five spectrograms of this star were obtained with the 
82-inch McDonald reflector, equipped with the standard Cassegrain spectrograph and 
500-mm camera. The dispersion was 55 A/mm at Hy, and the emulsion was Eastman 
103a-O. Table 1 shows the list of star lines and Table 2 the measured velocities. Only the 
spectrum of the primary is visible, except on two plates, Nos. 5384 and 5385, which show 
a few lines belonging to the secondary. Since the light diminishes to 0.115 of the normal 
during totality, the spectrum of the secondary is very faint. The spectral lines generally 
are diffuse, but an average of twenty-five lines per plate was measured, except on poor 
plates. The probable error per plate, for all but the poorest plates, varied from about 
+2.2 to +4.3 km/sec. 

Table 2, fifth column, shows the measured velocities. When plotted, however, the 
scatter seems to be larger than one would expect from the probable errors. A preliminary 
orbit was computed, and the residuals, V, — V., are shown in Table 2, sixth column, 
except for eclipse points, and are plotted in Figure 1 against the time in days. It will be 
noted that the grouping of points suggests that there were systematic errors present 
when the plates were obtained, and this is confirmed by the appearance of the lines of 
the comparison spectrum, which show a slight coma on plates taken between June 14 


* Contributions from the McDonald Observatory, University of Texas, No. 169. 
1 4.N., 185, 255, 1910. 

2 A.N., 196, 395, 1913. 

3 Acta astr., Ser. c, 2, 59, 1934. 


Variable Stars (““Harvard Observatory No. 5 (Cambridge, Mass.: Harvard College 
Observatory, 1938}), p. 72. 
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and July 5. The measured velocities were therefore corrected on the basis of the smooth 
curve shown in Figure 1 and the adjusted velocities are tabulated in Table 2, seventh 
column, and plotted in Figure 2. 

The phase was computed from the epoch JD 2432340.482 with a period of 2.806 days. 
Mid-eclipse would then occur at JD 2432341.885, which is in agreement with the two 
plates showing lines of the secondary spectrum, but is 0.146 day later than the epoch of 
minimum light as computed from the light-elements of Kordylewski. 

According to Baker® the eclipse phase lasts 0.460 day. Eleven eclipse plates, taken be- 
tween phases 1.2 and 1.6 days, were accordingly omitted in computing the orbit, leaving 
a total of forty-four plates. From spectrographic orbits published by Plaskett,® and 
Pearce,’ and the plot of velocities in Figure 2 it appears that e is very small. A circular 
orbit was therefore computed, using the adjusted velocities shown in Table 2, eighth 
column. The systemic velocity, y, and semi-amplitude, K, were determined by setting 
up forty-four equations of the form 


Adjusted vel. = y+ K sin CF 


phase ) 


TABLE 1 
LIST OF STAR LINES 


Element | Element | Element 


| 3856.37 | 

| 3878.58 || 
3902.95 

3913.46 || 

| 3916.42 

| 3933.67 

3956.68 

| 3997.39 || 

4003.33 || 

..| 4005.25 |) F 


and solving for y and K by least squares, which gave 
y= + 2.75 km/sec, K = 61.23+0.85 km, 


a sin i= 2,365,000 km, sin? i = 0.0680. 


me 
(m+ m2)? 
The smooth curve in Figure 2 is based on these elements. The rotation effect is evidently 
small. For comparison the solutions by Plaskett and Pearce are shown in Table 3. 
Although eighteen lines were measured for radial velocity on the two plates which 
show the secondary spectrum, there are only a few lines which are clear enough to permit 
classification of the spectrum. The most striking feature is the strength of Fe1. The 
plates are very weak on the violet side of \ 4325, but the Cam K line is visible. Hy is 
present but is distinctly weaker than Fe 1 4325 and 4383. No other H line is seen, and the 
G band is barely visible. 


I am much indebted to Dr. O. Struve for having obtained the spectrograms and for 
his kind advice and encouragement. 


5 Laws Obs. Bull., No. 33, 1924. 
5 Pub. Dom. Ap. Obs., Victoria, 1, 145, 1922.. 7 Proc. December meeting A.A.S., 1937. 
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TABLE 2 


RADIAL VELOCITIES 


Date U.T. 
(1947) 


In Days 


PHASE 


In Period 


MEASURED | 


VELOCITIES 
(Ka /Sec) 


| 
| 


RESIDUALS 
Ve 


| 
| 


(Km /Sec) | 


ApjusTEep 
VELOCITIES 
(Ka /Sec) 


June 


— 


2252 
2:50 
2:44 
3:06 
2:40 
3:00 
3:09 
2:41 
3:01 
3:08 
2:37 
2:57 
3:22 
4: 26 


— 


0.765 
.770 


—55.4 
—72.2 


+1 


MOR 


++++ 


++++ 
ONAN 


* The letter ‘““P” in this column denotes poor spectrogram. 
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Fic. 1.—Residuals from preliminary orbit 
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Fic. 2.—Velocity-curve of TW Draconis (phase in days) 
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ANOMALOUS SPECTRA OF STARS OF CLASS A 


Miriam E. WALTHER 
Harvard College Observatory 
Received April 15, 1949 


ABSTRACT 


Objective-prism plates with a dispersion of 95 A/mm at Hy have revealed 41 peculiar A stars and 8 
metallic-line A stars. The material examined shows that, at low galactic latitudes, 13 per cent of the A 
stars between magnitudes 5.5 and 9.5 have spectra characterized by abnormal intensity of the lines of 
chromium, silicon, europium, or strontium. 


I. INTRODUCTION 
In the course of the classification of objective-prism spectra taken recently at the 
Harvard Observatory, I have found a number of A stars with anomalous spectra. The 
plates were taken as part of a program designed to adapt the Morgan-Keenan-Kellman 
system of spectral classification to objective-prism spectra. The spectra, taken with the 
12-inch Metcalf doublet (the MA) equipped with a crown-glass objective prism of 12°, 
have a dispersion of 95 A/mm at Hy. Experiment showed that Eastman emulsion 
IIa-O and development in DK20 for 15 minutes at 65° gave spectra of highest quality. 
Three luminosity classes can be assigned to most spectra obtained in good seeing, but a 
detailed luminosity classification of the A stars is possible only on plates of superior 
quality. On the other hand, certain peculiar features are very marked, and these are 
easily distinguished from luminosity effects on all plates. 


II. THE PECULIAR A STARS 
The spectra of the peculiar A stars (as distinguished from the metallic-line stars) ap- 


pearing on the MA plates fall into three groups, which are characterized by the element 
that gives rise to the most marked anomalous feature. These groups—the silicon, the 
chromium-europium, and the strontium stars—are represented by prototypes which 
have been described in detail by Morgan.! 

Table 1 contains a list of the 41 peculiar A stars detected on the MA plates; only 2 of 
these stars were previously known to be peculiar. The first five columns of the table re- 
quire no explanation. The photovisual magnitude in the sixth column is taken from the 
Henry Draper Catalogue when the star is included in that catalogue; otherwise it is taken 
from the BD catalogue. The new spectral type is based on the intensity of the K line 
relative to the near-by Balmer lines. Since this line is weak in nearly every case, all the 
stars, with one exception, have been designated either B9 or AO. This classification is not 
altogether satisfactory, since it conveys very little information about the general ap- 
pearance of the spectrum. The descriptions in the notes which follow the table are more 
informative. Following to some extent the notation adopted by Deutsch,’ I define, be- 
low, the symbols describing the peculiar features listed in the table: 


“Cr” means that 4077 and \ 4171 of Cr 1 are strong. 

“Sr” means that \ 4077 and A 4215 of Sr 1 are strong. 

‘*) 4200” means that this unidentified line is strong. 

“Ew” means that Eu 1 is implied by the fact that the general features of the spectrum match 
a typical chromium-europium star. 

“Si” indicates that the blend at AA 4128-4131 is strong and that the general features of the 
spectrum correspond to a typical silicon star. 


1 Morgan, Keenan, and Kellman, An Allas of Stellar Spectra (Chicago: University of Chicago Press, 
1943); W. W. Morgan, Ap. J., 77, 77, 1933; ibid., 74, 24, 1931. 
2 Ap. J., 105, 297, 1946. 
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TABLE 1 
LIST OF PECULIAR A STARS 


| é j Peculiar 
(1900) Features 


| 459° 163 +59° 55’ Cr, Eu, Sr 
.| +57 184 Si 

+60 194 
+61 495 
+58 538 
+59 582 
+59 597 
+60 640 
+22 518 
+44 1006 
+44 1036 
+43 1092 
+29 893 
+15 1175 
+14 1276 
— 1 1395 
— 7 4623 
+14 3854 
182381 +15 3798 
189394 +34 3817 
190275 +37 3735 
192606 +34 3934 
193344 +35 4059 
196178 +46 2977 
197374 +43 3687 
+29 4202 
+45 3303 
+48 3260 
+43 3786 
+46 3144 
+44 3701 
+52 3056 
+45 3736 
+53 2766 
+43 4119 
+44 4031 
+55 2679 
+52 3127 
219855 +57 2719 
220147 +61 2430 
221568 +57 2758 


3 
< 


ss 
SUS 
Grou 


nN 


Cr, Eu: 
d 4128, d 4131 
Cr, Eu 
Cr, Eu 


NOTES TO TABLE 1 


1. On one of the better plates both Sr 1 4077 and Sr 11 4215 are very strong, while on another plate 
of nearly equal quality these lines are weaker, especially \ 4215, which is very weak or entirely absent. 
This suggests that the star is a spectrum variable. 

8. The spectrum of this star resembles that of x Serpentis, though the feature at AA 4128-4131 is 
stronger. Except for the Balmer lines, the Sr 1 lines A 407 and d 4215 are the most prominent features in 
the spectrum. The line at \ 4161, which may possibly be dye to Sr m1, is present. The weakening of 
dA 4077, 4161, and 4215 on two of the four plates, all of relatively high quality, suggests that this star is 
a spectrum variable. When Sr 11 4215 is compared with Cr 11 4171, it is found that it appears considerably 
stronger on two plates and equal or slightly stronger on the other two. 

9. Certain changes indicate spectral variability, though neither of the two spectra available is of 
sufficient quality to establish variability. 

15. The lines \ 4026 and \ 4077 also appear, but they are not strong. 

17. The spectrum of this chromium-europium star is similar to that of BD—18°3789 at phase 0.38 day. 
Strong lines appear at \\ 4077, 4128-4131, 4171; a fainter one seems to be present at \ 4179. The char- 


No. | HD No. | (1900) 
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NOTES TO TABLE 1—Continued 


acteristic asymmetry of Hy, Hé, and He, which is produced by blends on the red side of these lines, is 
very marked. The K line is very weak, corresponding approximately to class B9. 

19. Inferior spectrum..He 1 4026 present. 

21. The peculiarities of the spectrum of this star are so unusual that it is not possible to place it in 
any of the groups described previously. The K line is weak relative to the general development of the 
spectrum; it corresponds approximately to class A5. The metallic lines are, in general, strong, correspond- 
ing to class A7 or FO. A blend at \ 4272 and a line at \ 4144 are even stronger than the corresponding 
features in the FO stars. The unusual strength of these latter two features makes one suspect that the 
spectrum is a composite, but it is hard to understand what types could combine to produce this anomaly. 
The spectrum appears on four plates. 

23. The spectrum resembles very much that of BD —18°3789 at the phase 0.38 day, where the Cr 11 
lines at \ 4077 and \ 4171 and the line at \ 4179 are at maximum strength. The blend at Ad 4128-4131 
is strong; the K line is very weak. 

33. The spectrum, with peculiarities at \\ 4063, 4077, 4128-4131, and 4171, resembles that of 
BD — 18°3789. The relative intensity of \ 4063 and A 4077 is not the same on two of the plates. While the 
spectra on these plates are not comparable, it is almost certain that they give evidence of variation. 

36, 37, 38. The inferior quality of the spectra of these stars permits only an indication of the type of 
peculiarity; hence the classification should be considered provisional. 

39. The blend at AA 4128-4131 is not very strong but probably sufficiently so to qualify this star as 


peculiar. 
41. There is some indication that the intensities of \ 4171 and the K line vary, but this observation is 
subject to confirmation. 


The classifications given are provisional and subject to modification when slit spectra 
are available. 

Peculiarities occurring in the spectra of A stars are immediately apparent, but it is 
often difficult to identify the ions responsible for the peculiarities in a given spectrum. 
Certain peculiarities are so pronounced that they may be detected even when the quality 
of the spectrum is inferior. This is especially true of the blend at AA 4128-4131, which is 
so frequently recorded in the Henry Draper Catalogue as the only peculiar feature. The 
high percentage of silicon stars doubtless arises, in part at least, from the fact that 
dA 4128-4131 appears when many of the other fainter features, such as \ 4077 and 
\ 4171, are washed out. Data from the Henry Draper Catalogue reveal that 114 of the 
188 peculiar A stars in this catalogue belong to the group of silicon stars. 

In general, the strength of the blend near \ 4130 is attributed to the silicon doublet at 
d 4128 and A 4131. Certainly, one can assume this to be true when the remainder of the 
spectrum resembles a typical silicon star, such as BD+33°1008 Aurigae, in which 4200 
is strong and Ca K corresponds to class B9. On the other hand, one of the stronger 
europium lines, Eu 11 4129, may contribute to the strength of the blend at Ad 4128-4131. 
Actually, none of the strongest lines of europium, at AX 3930, 4129, and 4205, has been 
detected in any of the spectra thus far examined. However, when both AA 4128-4131 and 
4171 appear as prominent features and the spectrum resembles that of a typical 
chromium-europium star like 78 Vir, the existence of features due to europium is implied. 
I have, therefore, designated a star as belonging to the chromium-europium group if the 
spectrum resembles that of 78 Vir. 

The intensity of the feature at AA 4128-4131 differs greatly from one silicon star to the 
next. In fact, in accordance with Morgan’s previous observation,’ it appears from the 
MA plates that the intensity of the feature at \\ 4128-4131 in the silicon stars deviates 
only in various degrees from that of a normal A star. In the case of some stars it was 
difficult to decide whether or not they should be classified as peculiar; in others the 
peculiarities are very pronounced. Most of the doubtful cases have not been included in 
Table 1. The chromium-europium and strontium stars, on the other hand, seem more 
likely to form discrete groups. ' 

It would be well to emphasize that there exists very little possibility of confusing the 


3 Pub. Verkes Obs., 7, Part II, 7, 1934. 
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highly luminous A stars with the chromium-europium stars. In the spectra of both types 
there are features at approximately \ 4130 and \ 4171, but they have distinctly different 
appearances in the two types. In the chromium-europium star, AA 4128-4131 is definite- 
ly a strong blend, and A 4171 is a fairly narrow line due to Cr 11. On the other hand, in the 
supergiant the blend at A\ 4128-4131 is relatively weaker and narrower, and the 
feature near \ 4171 consists of two strong Fe 1 lines at \ 4173 and \ 4178. The line at 
\ 4077, which is strong in the chromium-europium stars, is so faint that it is hardly 
visible in the supergiants of classes AO and A2. The hydrogen lines also are very different 
in the two types of stars under consideration. They are narrow and weak in the highly 
luminous stars; in the chromium-europium stars they are wide and strong, and some ap- 
pear asymmetric. The K line, on the other hand, is relatively stronger in the supergiants 
than in the peculiar A stars of corresponding spectral type. Only on inferior plates could 
these spectral features be so distorted that the difference would not be evident. 

Table 1 includes three stars that seem likely to be spectrum variables and two others 
that show indications of variability. Two of the first three, BD+59°163 and BD+60°640, 
belong to the chromium-europium-strontium group. For both of them it is the strong 
strontium lines at \ 4077 and d 4215 which give evidence of variability. The third one, 
BD+45°3736, appears to be a spectrum variable similar to BD — 18°3789.4 In the case of 
the other two stars—BD+22°518 and BD+57°2758, both chromium-europium stars— 
the observations are less convincing. For each of the stars the spectral variations should 
be confirmed with slit spectra which show more detail. In general, objective-prism plates 
are not satisfactory for studying spectral variations, especially when these variations are 
not very pronounced. 


III. FREQUENCY OF PECULIAR A STARS 


Of the AO and A2 stars brighter than the sixth magnitude (app. mpy), the Henry Draper 
Catalogue lists 11 per cent as peculiar. The frequency of A stars in the HD Catalogue de- 
creases to 3 per cent for stars between 6.0 and 7.0 mag., and to less than 1 per cent for stars 
between 7.0 and 8.0 mag. This rapid decrease in frequency is certainly a selection effect. 
The MA plates reveal that many more of the faint early A stars, including classes B9—-A3, 
are peculiar. All the regions listed in Table 2 are very close to the galactic circle. In the 
fourth column there is given for each of the regions the total number of stars classified on 
the MA plates in the spectral classes B9, AO, A1, A2, and A3; in the last column are the 
corresponding counts of the peculiar A stars. Table 2 does not include all the peculiar A 
stars in Table 1, since some of them are in regions for which the classifications are not 
complete. The stars range in magnitude from 5.5 to 9.0, but the material is probably only 
relatively complete down to 8.5. Although it would be difficult to be certain that the 
counts are complete to any given magnitude, it is reasonable to assume that the classi- 
fications and counts are comparable in the degree of completeness for both the total 
number of B9-A3 stars and the number of peculiar A stars. From these data it follows 
that 13 per cent of the A stars at low galactic latitudes are peculiar. Since the data are 
rather scanty, the actual percentage must be recognized as preliminary. 


IV. METALLIC-LINE A STARS 


Although the metallic-line A stars have been noted for their peculiarities since the 
earliest work on spectral classification at Harvard, Titus and Morgan’ were the first to 
recognize that these stars form a fairly homogeneous group. Usually, though not al- 
ways, Miss Cannon classified these stars as composites. Miss Maury first described the 
characteristics of the composite spectra. These puzzling spectra usually represent a 


4 Deutsch, op. cit., p. 285. 
5 Ap. J., 92, 259, 1940. § Harvard Ann., 28, Part I, 48-49, 1897. 
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combination of a late-type and an early-type star with strong hydrogen lines. When 
these lines appear in a spectrum which is otherwise a later type, the K line is invariably 
narrow and obscure. If the hydrogen lines are very strong, the K line appears as a shal- 
low band, sometimes having a fine central core. The lines characteristic of the late type 
disappear from the ultraviolet, while they are strong in the interval from HB to Hé. In 
the metallic-line A stars, on the other hand, the metallic lines appear throughout the 
whole spectral range used. The K line in the two types of spectra is entirely different. In 
the composite spectrum the K line is broad and fuzzy, while in the metallic-line A stars 
it is sharp and distinct. These characteristics usually suffice for distinguishing these two 
types of spectra. 
TABLE 2 


COUNTS OF A STARS IN THE MILKY WAY 


CO-ORDINATES OF REGION CENTER ToTaL 
No. oF 


REGION j PECULIAR 
A Stars 


a (1900) 6 (1900) 


+++ 


Wat 


+: 
+- 
+ 


Total.... 


The criteria that I have used in distinguishing the metallic-line A stars from other 
anomalous A stars follow closely those set forth in a recent report from the Yerkes 
Observatory.’ Since the present study was nearly completed at the time that this report 
came out, I did not find it possible to use the list of standards published. While there seem 
to be fundamental differences between the metallic-line and peculiar A stars, there is 
occasionally a star whose spectral features at low dispersion do not place it definitely in 
either category. Two of these stars whose classifications are uncertain are included in 
Table 3, which lists the metallic-line stars discovered on the MA plates. My two classi- 
fications depend, respectively, on the relative intensities of H and K and on the strength 
of the other metallic lines. The notes following the table indicate precisely what the ob- 
served peculiarities are. 

Most of the stars in the spectral range AO-FO are classified primarily on the basis of 
the relative strength of K and H + He; consequently, many classifications do not de- 
pend on the rest of the spectrum. Unless the spectra are of superior quality, it is unlikely 
that the anomalous strengths of the metallic lines would be apparent. It follows that 
many of the metallic-line A stars may not be detected. 

For the stars listed in Table 3, the feature at \ 4077 is strong in every spectrum except 
in that of BD+59°621, in which it does not appear. If \ 4215 were also strengthened, 


7N. G. Roman, W. W. Morgan, and O. J. Eggen, Ap. J., 107, 107, 1948. 
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one might be inclined to assign these stars to the strontium group of peculiar stars. The 
line \ 4215, however, appears normal in every star in which A 4077 is strengthened. Sev- 
eral of the stars, notably BD+30°3875, BD+43°3789, and BD+ 28°787, substantiate the 
reciprocal relation between Ca 1 K and ) 4077 recognized by Weaver.* 

References’ are occasionally made to spurious absolute-magnitude effects apparent in 
the metallic-line A stars, but in no case are the features which show these effects specified. 


TABLE 3 
LIST OF METALLIC-LINE A STARS 


3 | | MEW 
MEW No. HD No. | BD No. | My | : Type | 


1900 (Metalli 
( ) (K Line) | fetallic 


Lines) 
19342 | +58°564 | 58°22") | FO 
20095 | +59 621 3 8.7) +5937) 84 | AO | | AS? 
22538 | +18517 | 332.4| +19 1| 7.67 | AS | FO 
| 


20039 | +16 559 | 3) +16 16 | A2? 
35035 | +28 787 | 516.7} +2822] 7. A3 | A3 | FO 
43682 | +16 1080 | +15 59} 8.: A2 | AS | FO 
190537. | +30 3875 | +305 A! A3 | FO 
200407 = +43. 3789 58. +43 A2 | AS | FO 


NOTES TO TABLE 3 


1. Sr 114077 and the blend Ad 4030-4034 are relatively strong for class FO. The ratio \ 4417: \ 4481 
seems to indicate class V. 

2. Although the K line is weak, corresponding to class B9, there are many other features which indi- 
cate a much later type. In the interval between 6 and He the lines and blends match those of an FO star, 
except that \ 4077 does not appear. The lines between Hy and H6 are faint, and it is not certain whether 
they correspond to those of an F star or not. While the lines between H8 and Hy are faint, their positions 
seem to coincide with the lines of an FO spectrum. Except for the failure of \ 4077 to appear, this star 
would qualify reasonably well as an early metallic-line A star. It is possible that this enigmatic spectrum 
is a composite. 

4. The spectrum lies near the edge of the plate and is underexposed. The K line is very weak or absent; 
the metallic lines indicate spectral type A2 or later. The line \ 4077 is rather strong. 

6. The quality of the spectrum prevents definite classification as a metallic-line star. The K line 
corresponds to AS, while the remainder of the spectrum corresponds nearly to class FO, though A 4077 


is abnormally strong. 
7. Sr11 4077 is very strong. From the ratios \ 4444:) 4481 and \ 4417: 4481 the star appears to 
belong to class V. 


For the stars in Table 3 the usual ratios employed by Morgan, Keenan, and Kellman? 
in the luminosity classification of FO stars indicate luminosity class V. In every case in 
which the ratio \ 4417: 4481 could be estimated, it corresponded to luminosity class V. 


1 wish to thank Dr. Harlow Shapley, Dr. Bart J. Bok, and Dr. Armin Deutsch for 
their encouragement and advice. It is also a pleasure to acknowledge the financial as- 
sistance that several fellowships have provided for my studies. During the first year I 
held the Edward C. Pickering Fellowship and a National Fellowship of the Kappa 
Kappa Gamma Sorority. The second and third years were financed, respectively, by an 
American Association of University Women Fellowship and by the Alice Freeman 
Palmer Fellowship. 


8 Pub. A.S.P., 58, 246, 1946. 9 Op. cit. 
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ABSTRACT 


Nine new bands have been observed in the spectrum of the C2 molecule between 2370 A and 3855 A. 
They are degraded to the red and do not have pronounced heads. The wave lengths of the heads are 
included in Table 11. An analysis of the bands has shown that they are additional members of the Fox- 
Herzberg system (B*IIl,—X‘II,). Their discovery makes possible a more complete discussion of the 
molecular constants of the Bll, state. The resulting rotational constants, Bj and aj, are 1.1922 cm™ 
and 0.0242 cm, respectively. The vibrational quanta are included in Table 15, and the resulting vibra- 
tional constants in Table 14. Since the vibrational quanta show an anomalous change with increasing 
vibrational number, extrapolations to higher vibrational levels have little significance when made on the - 
basis of these constants. 

A. INTRODUCTION 

During an investigation carried out recently in this laboratory on the spectrum of 
the C. molecule, a number of faint bands were found between 2370 A and 2855 A which 
have not been previously mentioned in the literature. A low-dispersion spectrogram of 
this region is reproduced in Figure 1, spectrum a. The heads of the more prominent of 
the bands are shown. The bands degrade to the red and do not have pronounced heads. 
An enlargement of a high-dispersion spectrogram of a portion of the band at 2772 A is 
shown in Figure 1, spectrum 8. At this dispersion the fine structure of the band is well 
resolved. The lines are readily grouped into two series of triplets, forming the R and P 
branches, identified in the figure by the short and long leading lines. Furthermore, there 
is a regular alternation in the spacing of the components of successive triplets in each 
series, alternate triplets showing the same relative spacing. All the bands identified in 
spectrum a, Figure 1, have the same structure as the 2772 A band. 

The triplet structure of the bands and the alternation, or “stagger,” of the spacing of 
the component lines of the triplets are characteristics of bands of the C2 molecule arising 
from transitions between triplet energy levels. Triplet systems that have been observed 
are the Swan system (A*II, — X*If,) and the Fox-Herzberg system (B*II, — 
Four bands of the latter system, namely, the (0-3), (0-4), (0-5), and (0-6) bands at 
2855 A, 2987 A, 3129 A, and 3283 A, respectively, were found and analyzed by Fox and 
Herzberg! in 1937. They are degraded to the red and do not have pronounced heads. The 
similarity in appearance and structure of the new bands and the four previously identi- 
fied bands of the Fox-Herzberg system and the fact that they are found in the same 
spectral region suggest that the new bands are further members of this system. If this is 
actually the case, the new bands would be of considerable interest, since undoubtedly 
some of them would arise from excited vibrational levels of the B*II, state and would thus 
provide a means for making a more complete determination of the vibrational and rota- 
tional constants of this state. 

B. EXPERIMENTAL 


The apparatus and techniques used to produce the spectrum of the C2 molecule have 
been described in a previous publication.” The high-dispersion spectra were photographed 
in the fourth order of the 21-foot grating spectrograph, with a dispersion of 0.625 A per 
millimeter. The emulsions used were 103a-O and 103a-O UV. The exposure times varied 
from 5 hours at 2855 A to 10 hours at 2300 A. 


1 J. G. Fox and G. Herzberg, Phys. Rev., 52, 638, 1937. 


2 John G. Phillips, Ap. J., 107, 389, 1948. 
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Cc. ANALYSIS OF THE NEW BANDS 


For reasons given in the introduction it was suspected that the new bands were further 
members of the Fox-Herzberg system. If this was actually so, the band heads of the nine 
most prominent bands could be used to determine the vibrational transitions involved. 
Little error was made in using the band heads, since the rotational lines in the neighbor- 
hood of the heads were very weak, indicating that the band origins were very close to 
the heads. The order of magnitude of the vibrational quanta of the upper state could be 
predicted in advance, since an approximate proportionality exists between the vibration- 
al and rotational quanta of a given electronic state, according to which a state with 
smaller rotational quanta will have smaller vibrational quanta. Thus, since Fox and 
Herzberg! found that the constant Bj of the B*II, state was 1.1804 cm~, while the 
Swan bands gave B;’ = 1.6326 cm~. and w, = 1641.35 cm for the XII, state,® it 
could be predicted that w, of the B*Il, state was approximately 1000 cm~. With this 
as a guide, it was found that it was indeed possible to find a consistent set of vibrational 
numbers for all the nine bands. These tentative numbers are included in Figure 1, spec- 
trum a. 

The tentative vibrational analysis was then verified by searching for combination 
differences between the R and P branches of each new band that were in agreement with 
corresponding combination differences found for Swan bands with the same lower vi- 
brational state. This search was successful in all cases, proving conclusively that the 
identifications of the lower states that are involved are correct. The relative numbering 
of the upper vibrational levels must also be correct, since, as will be shown later, the re- 
sulting vibrational quanta of the upper state fit into a smooth series. The only uncer- 
tainty that remains is that the absolute vibrational nurabering of the upper state may 
be in error and that they should all be increased by one or two units. However, it will be 
shown later that this is unlikely. 


a) ROTATIONAL ANALYSIS 


The rotational quantum numbers of the lines of each band were found by the process 
described above. Tables 1-9 give the measured lines of each band, and their J-values. 
The assignment of the red, central, and violet components of each triplet to the *ITp—‘*IIo, 
‘JI, — *Ih, and *ITz — *ITlg subbands, respectively, was made unambiguously from a 
comparison of the corresponding combination differences of the Swan system. As an 
illustration of the degree of agreement found, combination differences of the lower states 
of the (1-2) and (2-2) bands are compared in Table 10 with those obtained previously* 
for the (0-2) Swan band. In this table, columns headed A2F}’(J) are to be compared with 
one another, etc. 

With the rotational analyses known, the usual methods were used to find the rota- 
tional constants and band origins of each of the nine bands. Three steps were taken, as 
follows: 

1. Using the method previously described in the study of the Swan bands,’ approxi- 
mate rotational constants of the upper and lower states of each band were found, using 
the following equation: 


AF n (J) 
J+} 


2. The band origins, vo, were found through the use of the equation 


R, (J -1) +P, (J) —2 +1)? 


= 4B,—8D,(J +4)? 


D’’) 


corr 


= 2v, +2 (Bi —B") J2(J+1)?. 


v ‘corr 


3 Phillips, ibid., 108, 434, 1948. 
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TABLE 1 


WAVE NUMBERS OF THE LINES IN THE (1-3) BAND 
(vo = 36,058.97 Cm) 


PJ) || 


Cm | Cm | 
36,053.60 36,027.29 
052.32 | 


| 
| 
| 
| 


| 
| Cort | Get | | 
6....| 36,056.51) 1026.33 || (064.31 | 36,034.11 
| 048. 37 013.23 || 7....| 053.60; 019.09 || 8... 060.20) 025.52 
7....| 044.92 | 36,005.53 8....4 048.97 009.36, 9 055.22; 015.98 
8....| 039.46 | 35,995.64 |) 9....| 044.15 36,000.28 | 10... 049.66 | 36,005.53 
9....| 034.58} 985.95 || 10....! 037.81 | 35,989.26 || 11... 043.12 | 35,994.61 
10....| 027.29] 974.25 || 11....| 031.43} 978.26 || 12... 035.99 982.83 
020.51! 962.76 | 12... | 023.30) 965.54 13...) 027.84 969.99 
12....} 011.70} 949.25 13....| 015.27 | 953.01 || 14...) 019.09 | 956.80 
13... .| 36,003.01} 936.12 | 14....| 36,005.53} 938.55 15...| 36,009.36 | 942.44 
14....| 35,992.48 | 921.02 || 15... .| 35,995.64 924.26 || 16...| 35,999.11 | 927.60 
15....| 982.05| 906.01 | 16....| 984.11} 908.02 || 17...| 987.67 | 911.58 
16....| 969.99} 889.24 17....| 972.77 | 892.14 || 18...| 975.76| 895.10 
17....| 957.62 | 872.43 || 18... 959.45} 874.29) 19...| 962.76| 877.42 
18....| 943.85} 854.02 | 19....| 946.51} 856.73 || 20...; 949.25 859.37 
19....| 929.86} 835.51 20....| 931.50) 837.09 | 21...) 934.35} 840.02 
| 914.40, 815.49 | 21....| 916.86) 817.97 | 22...) 919.24 | 820.33 ae 
21....; 898.65 | 795.24 || 22....; 900.08 796.60 |) 23..., 902.69 799.25 
22....| 881.57| 773.54 |] 23....! 883.83] 775.85 |} 24. 885.98 | 777.99 
23....| 864.04} 751.48 || 24....| 865.30 752.74 25...| 867.77 755.21 
24....| 845.30! 728.24 || 25.... 847.43 730.38 || 26. 849.41 732.35 
25....| 826.06! 704.44 |) 26.... 827.17 705.58 || 27. 829.49 707 .89 
26....; 805.69 679.62 || 27....| 807.71 681.64 28... 809.55) 683.40 
27....| 784.65} 654.04 |) 28 785.66 655.08 29. 787.86 657.23 
28....| 762.69 627.62 || 29.. 764.62 629.55 30. 766.27 631.21 
| 739.86 600.31 30... 740.78 601.30 31...| 742.80 603.26 
30....| 716.28 | 572.31 |} 31.. 718.09 574.11 || 32.. 719.70 575.62 
31....) 691.71. | $43.27 || 32....| 692.52 544.04 33... 694.52 546.02 
32...., 666.53) $13.63 668.22 515.37 || 34..., 669.81 | 516.87 | 
33....) 640.19 482.84 | 34....| 641.33 483.54 |) 35... 643.74 485.43 
a 34....| 613.38 | 451.60 || 35....) 614.97 453.24 36...) 616.04, 454.73 
35....| 585.28} 419.06 || 36.... 585.93 420.38 | 37...| 587.47 | 422.45 
556.84 386.20 37....| 558.40, 387.80 | 38...| 559.55 | 388.78 
$27.02 | 352.00 || 38....| 527.53 | 352.59 |) 39...| 529.06 | 354.06 
38....) 497.00, 317.55 |] 39....| 498.41 | 319.02 |, 40...) 499.50} 320.12 
39... 465.34) 281.53 || 40....| 465.80 282.05 || 41... 467.28 | 283.54 
40....| 434.14) 245.54.) 41....| 435.04 246.98  42...| 436.08 | 247.99 
41....| 399.73 | 207.77 || 42....| 400.61 | 208.77 || 43. 402.03 210.17 
42...., 366.89 170.51 43... 368.29 171.47 |) 44. 369.25 172.46 
we 43....| 331.54] 130.10 || 44.. 332.06 130.97 45. 333.30 132.38 
44....| 296.67 091.24 45.. 298.09 092.71 | 46. 299.08 093.64 
45..... 259.68 46. 260.14 47. 261.44 
| | | | | | 
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TABLE 2 


WAVE NUMBERS OF THE LINES IN THE (0-2) BAND 
(v= 36,593.78 


| 


311, | | 


a4 

= | 7....| 36,580.51 | 36,540.97 || 8...| 584.54] 544.27] "501.16 551.11 
8....| 575.32] 530.38 || 9...; $80.03} $35.43 || 10....| 585.77 | 540.97 
9....| $70.50} $21.00} 10...) 573.90 | 524.43 || 11....] $79.42) 529.86 
10...., 563.56 | 509.36 M1... 567.69} 513.55 | 12...) 572.51] 518.20 
$57.02 497.97 |) 12...| 589.84} 501.17 13....| $64.63) 505.71 
12...) $48.45} 484.83 13...) 552.11 | 488.53 14....] 556.10] 492.55 
13...) 540.09) 471.92 || 14...| 542.66} 474.26 13....| $46.70] 478.45 
15..../ $19.68 | 441.82 || 16...] 521.78 | 444.09] 525.54]... 
16...) 507.87 | 425.55 |] 17...| 510.83 | 428.49 |] 18....| $13.90 | 431.57 
496.04} 408.84 |] 18...| 497.97 | 410.83 || 19...) 501.17 | 414.16 
18....| 482.42} 390.77 |] 19...) 393.81 || 20....| 488.06 | 396.28 

21....| 438.50 332.58 || 22...) 439.90 334.06 || 23....| 442.81 336.91 
a 22...., 421.95 | 311.32 |] 23...| 424.31 | 313.66 |] 24... 426.50 | 315.96 
23... 404.82 | 289.60 || 24...| 405.89 | 290.92 |) 25...) 408.84 | 293.58 
24... 386.47} 266.71 || 25...| 388.67 | 268.94 || 26...) 390.77 | 271.08 
25....| 367.69 | 243.36 || 26...| 368.86 | 244.60 || 27....| 371.28] 247.02 
26....| 347.76 | 218.87 |] 27...| 349.81 | 220.93 || 28....| 351.75 | 222.86 

27...) 327.29 | 193.80 || 28...| 328.35 | 194.87 || 29....| 330.62 | 197.18 
305.85} 167.69 || 29...) 307.78 | 169.70 || 30....| 309.59 | 171.40 
29....| 283.70} 141.00 |} 30...) 284.57 | 141.92 31...) 286.72) 144.10, 
Es 30...) 260.57 | 113.28 |] 31...) 262.45 | 115.17 || 32....} 264.10 | 116.80 
31...:| 236.72 | 084.77 |] 32...) 237.62 | 085.67 || 33....} 239.63 | 087.72 
32... 212.13} 055.75 || 33...| 213.89} 057.50 |] 34....] 215.35 | 059.02 
33....| 186.49 | 36,025.52 || 34...| 187.26 | 36,026.33 || 35....| 189.23 | 36,028. 28 
34...) 160.34 | 35,994.61 || 35...| 161.99 | 35,996.48 || 36....| 163.39 | 35,997.93 re 
35...) 132.95 | 962.76 || 36...; 133.68 | 963.62 || 37....| 135.52 | "965.54 
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TABLE 3 


WAVE NUMBERS OF THE LINES IN THE (2-3) BAND 
(vo = 37,038.52 


3Ilo—Ilo } 31, 


37,014.03 
37,007.14 
36,998.42 
989.70 
979.00 
968.78 
955.93 
943.82 
929.58 
915.59 
899. 23 
883.79 
865.82 
848.37 


Cm" Cm || Cm Cm" || Cm" 
9....| 37,011.27 | 36,963.64 |) 10... 36,966.67 || 11....| 019.00} 971.58 
10....| 37,003.51 951.26 | 11... 955.14 || 12....) 011.27 959.34 
11....| 36,905.64 939.25 || 12...) 941.85 || 13....| 37,002.59 | 945.90 
13....| 976.54] 911.29 | || 913.62 || 15....| 982.80} 917.17 
14....| 965.39} 895.48 || 15...| 898.58 || 16....| 971.58} 901.52 
15....| 954.02 879.68 | 16...| 881.53 || 17....| 959.34] 884.90 
16....| 940.96 | 862.07 || 17...| 865.28 || 18....| 946.61 867.62 
7....) 929.03 844.43 || 18. . .| 846.25 || 19....| 932.29 | 849.18 
18....| 913.02 | 825.50 || 19... 827.50 || 20....| 918.01 830.32 
19....| 897.97 805.69 || 20... 807.25 || 21....| 902.19} 809.85 
20....| 881.53 784.50 || 21... 787.14 || 22....| 885.92 789. 36 
21....| 864.32 763.28 || 22... 764.85 || 23...| 868.28 767.01 
22....| 846.25 740.52) 23... 742.55 || 24....| 850.31 744.73 
23....| 827.50 717.35 || 24...| 828.51 718.51 || 25....| 830.90} 720.77 
24....| 807.25 692.84 | 25... 809.39 | 694.95 | 26....| 811.20) 696.73 
25....| 786.65 667.85 || 26...| 787.99 | 668.92 || 27....| 789.91 671.03 
26....| 764.85} 641.64 || 27...| 767.01 643.78 || 28....| 768.54} 645.37 
: 
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TABLE 4 


WAVE NUMBERS OF THE LINES IN THE (1-2) BAND 
630.36 


| 


i | 
R3(J) P3(J) | 


Cm | & | 
37,592.20 || 6...¢...........| 37,597.68 || 7....| 37,635. 37,605. 1. 
37,618.35 583. | 595.85 
615. 5. _..| 37,619.05 =e 5. 586.0: 
575. 
562. 
552. 


506.98 |........ 


a 
| 
| 
| 596.59 | 538 
9.... "| 7 534.40 || 13.... 387 42 525.14 
10....) 531.63 || 12... 521.34 | 377.27 510.27 
317.79 | 43... 506.35 15....| 566.38 494.93 
488.38 || 15... 550.98 474.89 | 541.99 
94 538.97 COS 33 
15....) 548. 455.54 | 17...) 5 439.90 || 19. 514.14 | 424.3: 
{ 536. 12 525. 14 20. 
16....) 12 || 18... 421.81 1 404. 
| 419.09 || 341.57 401.53 || a. | 383.91 
8.87 419. 495.81 | 482. 
508. 399.94 || 20... 48 381.64 || 22.. 465.68 362.24 
i 358.13 || 22... O1| 338.05 || 24...) 448.12 
443.68 1 | 426. 290.97 || 26... | 267.55 
| | 425.41 312.9 407.98 96 | 389.13 | 42.15 
5.80 | 288.80 386.83 | 265.26 | 28....1 368.27 | 242. 
24....| 405.8 264.22 || 26.. 240.38 || 28....| 345.661 215.0: 
385.73 || 7 366.48 J4e | 87.99 
26....; 364.43 28... 186.37 || 30... | 70 159.13 
27....| 184.43 || 29...| 157.10 | 130.37 
577 156.19 || 30...| 272.86 128.82 |) 32... 248.27 | 099.76 
127.06 || 31.. 097.78 || 33... | | 060.32 
271.09 246.3 7.91 || 34.. 36.91 
20....| 6.99 || 32...| 067.91 |) 34 5. 26 056. 
31....) 245.51 066. 04 33... .| 035.03 || 35.... | 37,004.92 
33....| 163.59 | 37,001.85 | 35...| 135.03 | 36,968.78 | 37... | 107 40 936.5 
134.42 | 36,968.20 || 36. 106.13 935.60 || 075.67 | 
33: 104.66 934.03 | 074.12 40 044.6 
73.63 864.32 || 40...) 044.67 827. 
37....) 073. 862.75 || 39... 825.99 || 41..../ 37,011.27 790.41 
| 042.18 37,009.92 42... 36,978.5 750) 93 
825.50 || 40...| 37,009.9 789.36 | 943.26 750.9 
aN: 36,977.50 ia... 3. 712 08 
39... .| 37,009.3 | 908.89] 712. : 
| 36,976.54 | 42... 941.85 8 64 
| 940.96 907.8 4s... 71. 29 
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TABLE 5 


WAVE NUMBERS OF THE LINES IN THE (2-2) BAND 
(vp = 38,609.86 


311; — 311, 


| 
| 
| 


Cm" 
38 601.36 
596.17 
590.81 
583.62 
576.33 
567.07 
557.93 
546.90 
535.88 
522.89 
510.09 
495.23 


049.14 
38,015.49 
37,981.45 _...| 37,982.40 

945.43 5. 946.80 

910.96 16 || 42....| 911.83 

873.07 19 || 43... 874.39 

836.79 |. 937.66 

5 798.10 
759.63 
718.12 
677.92 


Pa) | Ri(J) | PA(J) = 
6....| 38,596.17} 562.51 || 7...| 567.95 || 8....| 38,608.05 573.99 
7....| 592.76] 554.39]1 8...| | 557.93 || 9....| 602.31 563.71 
8....| $86.33 | 543.42 || 9...| 548.00 || 10....| 596.17 552.93 a 
9....| 580.48] 532.72 || 10...| | 536.28 |] 11....1 588.55 541.08 
10... .| 572.34 | 520.41 || 11...| 524.39 || 12....] 580.48 528.53 
11....| 564.52 507.91 || 12.. | $10.70 || 13....} 571.28 514.83 
12....| 554.39 493.47 || 13...| | 497.02 || 14....| 561.48 500.55 
Se 479.19 || 14... | | 481.46. 15....| 550.53 485.05 
14....| 532.74 462.88 || 15. | 465.96 || 16....) 538.94 469.02 
15....| 520.90 446.54 || 16...) | 448.52 || 17....| 526.18 451.73 
16....| 507.27 428.29 || 17... 431.22 || 18....| 512.79 433.89 
17....| 493.47 410.08 || 18... 411.82 19....| 498.17 414.51 
18....| 478.02 390.04 || 19...| 480.63 392.78 || 20....) 482.95 395.05 
19....) 462.29 369.92 || 20... 463.84 371.60 | 21....| 466.49 374.16 
20...) 445.02 348.21 21...| 447.41 | 350.69 || 22....| 449.50 352.80 
326.17 || 22.. 428.69 327.17 || 23....| 431.22 329.94 
22....| 408.25} 302.59 |} 23...| 410.48 | 304.85 || 24... 412.42) 306.79 
388.73 278.59 || 24...| 390.04, 279.85 |) 25...) 392.15 282.08 
24....| 367.77 253.16 || 25.. 369.92 255.38 | 26....| 371.60 257.12 
25....| 346.24} 227.39 || 26.. 347.39} 228.50 || 27...:| 349.48 | 230.58 
26....| 323.57 200.26 || 27...| 325.54 202.26 || 28....| 327.53 203.85 
27....| 300.14 172.42 || 28...} 301.11 173.43 || 29....] 302.96 175.42 
28....| 275.60 143.56 || 29...| 27.49 145.48 || 30....| 278.94 146.97 
29...., 250.27 113.80 |} 30...| 251.21 114.76 || 31....] 253.16 116.48 
30....| 223.86 083.15 || 31..| 225.74 084.93 || 32....| 227.05 086.24 
31....) 196.71 051.56 | 32.. | 197.53 052.35 |) 33... 199.22 054.05 
32....] 168.55 | 38,019.04 |) 33...| 170.29 | 38,020.84 |) 34... 171.50 | 38,022.00 
33....| 139.39 | 37,985.53 || 34.. 140.06 | 37,986.28 || 35... 141.73 | 37,987.85 
34....| 109.48 951.32 || 35.. 111.06 952.91 || 36...) 112.15 954.07 
35...., 077.49 915.88 || 36.. 078.96 | 916.53 || 37....| 080.49 918.05 
36....| 046.49 | 879.86 || | 882.46 
37... .| 38,013.44 841.51 || 38.. 844.44 
38... .| 37,979.95 804.55 | 39.. | | 807.11 
39....) 944.88 765.27 || 40._ 767.21 
40... 909.57 725.69 || 41.. | | 728.07 
41....} 872.49 684.52 | 42...) | 686.33 
42....| 835.36 643.17 || 43. 
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TABLE 6 


WAVE NUMBERS OF THE LINES IN THE (3-2) BAND 


(vo= 39,542.94 Cm=!) 


—*ITo 


P(J) | 


Pi(J) 


39,365.95 
| 347.29 
327.03 
306.22 
284.00 
261.29 
237.15 
212.41 
186.27 
159.48 
131.56 
102.61 
072.84 
042.04 
39,010.15 
38,977.43 
943.73 
908 .90 
873.26 
836.21 
799.22 
759.73 
720.38 
769.46 
637.98 
594.89 


Cm"! 
39,357.91 


39,359.78 
340.38 
319.84 
298.04 
275.61 


| Rs(J) | Pus) RJ) | 
Cm"! | Cm || Cm | Cm" 
ag 20... 271.37 || 21...| 39,368.38 273.76 || 22....| 39,370.38 = 
21. 248.28 || 22...| 348.56 249.62 || 23....| 350.89 251.82 
22. 223.59 || 23...| 329.20 225.88 || 24....| 330.98 227.56 
23. 198.62 || 24...| 307.40 199.79 || 25....| 309.48 201.77 
24. 172.05 || 25...| 286.09 174.19 || 26....| 287.55 175.66 
25... 145.02 || 26...) 262.30 146.00 || 27....| 264.29 147.98 
26... 116.69 || 27...| 239.09 118.39 || 28....| 240.61 119.94 
087.46 || 28... 213.28 088.48 || 29....) 215.14 090.15 
28.. 057.16 || 29...| 188.65 058.99 || 30....| 189.57 060. 28 
29... 39,026.06 || 30...| 160.33 | 39,027.02 || 31....| 162.07 | 39,028.94 
30. 38,993.80 || 31...) 133.42 | 38,995.70 || 32....| 134.50 | 38,996.92 
31. 960.83 || 32...| 103.60 961.55 || 33....| 105.08 963.18 
32. 926.77 || 33...| 074.69 928.32 || 34....| 075.68 929.51 
33. 891.63 || 34...| 042.74 892.30 35....| 044.25 893.75 
34.. 855.60 || 35...| 39,011.95 857.28 || 36....| 39,012.80 858.41 
35... 818.50 || 36...| 38,978.04 819.18 || 37....| 38,979.50 820.46 
36... 780.76 || 37...|...........| 782.28 || 38....1 946.45 783.17 
37. 741.50 || 38...| 908.90 742.12 || 39....| 910.68 743.45 
38... 701.86 || 39...) 874.80 703.25 || 40....| 875.78 704. 39 
39... 660.70 || 40...| 836.79 661.25 || 41....| 837.87 662.37 
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TABLE 7 


WAVE NUMBERS OF THE LINES IN THE (2-1) BAND 
(vo = 40,204.49 Cm) 


| 


J 
| 
R2(J) P2(J) 


40,170.01 ee 40,177.14 
162.11 167.69 
151.56 
141.90 
129.53 
117.27 
102.53 
088.75 
072.65 
056.62 
038.64 
020.77 

40,000. 84 

39,981.19 
959.25 


Rs(J) P3(J) | Ri(J) | Pi(J) 
Cm" Cat Car | Car Cm" 
40,172.90] 5...1 40,204.10 |. 
5....| 40,194.57 166.33 || 6...; 200.03 
6....| 190.16 156.62 |} 7...| 195.81 | 
7....| 147.84 || 8...| 190.16 | 
4... 136.57 || 9...| 184.17 | 
9....| 174.05 125.99 || 10...| 176.58 | 
10...., 164.83 113.00 || 11... 168.80 
11....| 156.62 100.20 |} 12...| 159.20 | 
12....| 146.11 | 085.17 |} 13...} 149.62 | 
13....| 135.7 070.34 || 14...| 137.96 
14....| 123.41 053.41 || 15...| 126.52 
15....| 110.94] 036.59 | 16...| 113.00 
16....| 096.79 | 40,017.92 || 17...| 099.60 
17....| 082.55 | 39,999.13 | 18... 084.17 

18....| 066.33 978.49 || 19... 068.95 - 

19....| 049.91 957.66 || 20...| 051.45 

20....| 031.92 935.28 || 21...| 034.31 937.72 || 22....| 036.59 939.83 
21....| 40,013.62 912.44 || 22...| 40,014.96 913.77 | 23....| 40,017.44 916.23 
22....| 39,993.78 888.12 || 23...| 39,995.93 890.38 | 24... 39,997.84 892.25 

23....| 973.40 863.31 || 24...| 974.59 864.57 || 25..... 976.89 866.79 

24....| 951.61 837.12 || 25...| 953.78 839.22 || 26...., 955.48 840.95 
25....| 929.30 810.43 || 26...; 930.55 811.56 || 27...) 932.57 813.54 

26..... 905.59 782.35 | 27...) 907.52 784.33 || 28....| 909.23 | 785.91 
753.62 || 28...) 882.34 754.61 || 29....) 884.14 756.54 

723.71 || 29...| 857.70 725.62 30... | 859.07 727.09 
29...., 829.52 693.02 || 30...| 830.40 693.92 31....| 832.22 695.69 
30...) 802.05 661.19 || 31...| 803.86 663.07 || 32...., 805.16 664. 32 

628.64 || 32...) 774.57 629.49 33... 776.28 631.14 im 
32....| 744.49 594.99 || 33...| 746.14 596.65 34...) 747.41 597.87 

33....{ 714.13 560.06 || 34...| 714.89 561.12 || 35... 716.58 562.62 
34....| 683.04 524.89 || 35...| 684.82 526.51 || 36.. | 685.90 527.56 
35....| 650.06 }...........]] 36...] 651.43 488.34 || 37..... 652.92 490. 38 

452.52 || 38....| 620.21 453.59 

| i | | 
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TABLE 8 


WAVE NUMBERS OF THE LINES IN THE (3-1) BAND 
(vo= 41,137.86 Cm) 


44,922.07 
117.54 


41,014.98 
40,999.50 


| 


| 976 


41,011.12 
40,994 


431.99 


| Ri(J) | PJ) | | || | 
| Cm: || Cm | Cm" || | Cm" | 

7. 080.04 |) 8... 41,121.21 083.60 || 9... 41,126.75 | 088.82 
110.49 068.50 072.99 || 10....| 119.60} 077.41 
103.71 057.39} 10...) 106.79 | 060.35) 11... 111.37 | 064.81 
10....| 094.62} 043.72 || 11...) 098.39} 047.65 12... 102-30 | 051.44 
11....| 085.71 030.58 || 12, 088.28 | 033.04 || 13... +092 13 | 036.85 
12....| 074.55 13... 077.93} 018.37 || 14...) 081.19 O2t.45 
13.... 063.52 114...) 065.70 | 41,001.70 || 15....| 069.11 | 41,004.99 
14.... 050.34} 981.99 | 15..., 053.43 | 40:985.03 || 16....! 056.26 | 40,987.83 
15..../ 037.17 | 964.41 |] 16...) 039.09 $6.31 |] 17... 042.14 | “969.30 
16... 022.05} 944.92 || 17...) 024.81} 947.09 | 15. 37 950.09 
17....) 41,006.71 925.35 || 18...| 41,008.46 | 926.97 19.. 929.57 
ip 18... 40,989.63 | 903.79 || 19. | 40,992.19 | 906.33 |) 20. 908. 44 
19....| 972.32] 882.02 || 20...) 973.78} 883.50 || 21... 26 | 885.93 Be 

| 20....| 953.17 | 858.55 21...) 955.55] 860.92 || 22....| 957.48] 862.79 
21....| 933.76 | 834.79 || 22...) 935.07 | 836.09 || 23... | 937.33 | 838.31 
22....| 912.71 809.71 || 23.../ 914.92] 811.94 916.59) 813.64 
23....| 891.18 | 783.46 || 24...| 892.35 784.68 || 25...) 804.45 | 786.74 
24.... 868.05 795.97 25...| 870.14] 758.13 |) 26... 871.69 759.64 
25....| 844.46} 728.07 || 26...| 845.58 | 729.13 || 27... 847.53 731.09 
26..... 819.50} 698.64 27...| 821.39] 700.58 || 822.79 | 701.97 
ae 27....| 793.79 | 668.72 || 28...| 794.74] 669.68 || 29. 796.64} 671.40 
28....| 766.60} 637.42 || 29...| 768.53 | 639.27 || 30... 769.79 | 640.46 
29...., 738.75 | 605.41 || 30... 739.68 | 606.31 31... | 741.37 | 607.91 
30....| 709.76! 572.09 | 573.86 || 32...) 712.71 | $74.99 
31... | 679.90 | 537.82 || 32...) 680.75 | 538.67 || 33...) 682.24 | 540.30 
32....| 648.91 502.78 | 33...) 650.63 504. 39 34...) 651.67 | 505.42 
33...) 617.10} 466.48 34... 617.76 467.23 || 35...) 619.31) 468 
34... 584.03 | 429.37 35...| $85.55] 430.99 |] 36... | 

| 
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TABLE 9 


WAVE NUMBERS OF THE LINES IN THE (4-1) BAND 
(vo = 42,032.78 


— i} 311; 


RV) | Pa) | RJ) | 


Cort} Cert 
_..| 41,983.62 || 7... 
{FPO AS. 973.96 || 8... 41,977.61 || 9....| 42,019.31 | 41,982.07 
42,002. 962.04 966.30 | 10....| 011.84 970.48 
950.39 | 10... 53. -...| 42,002.97 957.32 
41,993.41 


4 

|| || 

10... 
976.54] 922.46 12...) 979.13 925.17 13... 982.79 928. 56 
964.73 | 906.65 | 13. 968.08 | 909.97 | 14... 971.22 912.64 
13....| 952.83 890.37 || 14...) 955.17 | 892.37 || 15....| 958.34 895.43 
14....) 939.07 872.13 942.12! 875.02 || 16....| 944.70 877.64 
853.89 16...) 927.05 855.57 || 17....| 929.84 858.41 
16..... 909.09 833.65 || 17...| 911.85 836.41 | 18....| 914.09 838.69 
17....| 893.01 813.20 | 18...| 894.81 814.79 $97.01 817.33 
790.94 || 19...| 877.64 93.28 || 20.... 
856.80} 768.39 | 20..| 858.41 | 769.68 | 21....| 800.37 | 772.14 
20....| 836.41 743.85 || 21...| 838.69 746.31 | 22... 840. 56 747.93 
816.01 719.11 || 22...| 817.33 720.34 || 23.... 819.35 722.41 
22....| 793.80 692.60 | 23... 796.00 694.79 || 24... 797.58 696.35 
665.67 || 24...) 712.44 | 666.84 || 25....| 774.20 668 . 66 
637.15 || 25...| 748.84 639.18 26...., 750.30 640.60 | 
608.17 | 26... 722.94 609.07 |. 724.76 610.89 
26....| 695.55 577.55 || 27...| 697.48 579.48 || 28....| 698.77 580.57 
27....| 668.66 546.30 | 28...) 669.44 547.20 | 670.95 548.92 
an 28....| 639.98 513.68 |729...| 641.99 514.72 || 30....| 643.13 516.49 
29....; 610.89 480.26 || 30... 611.60 481.05 |, 613.21 482.79 
30....| 580.57 445.72 || 31...) 582.23 447.32 32....| 583.25 448.30 a 
31....| $48.92 409.99 | 32...| 549.80 410.76 || 33....| 551.17 412.27 
32.... 516.49 373.29 | 33...) 518.20 374.99 || 34....] 519.22 375.81 
335.43 || 34...) 483.55 336.38 | 35...., 484.85 337.63 
34....| 448.30 296.90 35...) 449.76 299.05 || 36....| 450.62 299. 32 
35;.;.) 256.98 | 36.. 413.08 257.66 || 37....| 414.22 258.85 q 
216.24 || 37...| 376.95 217.82 | 218.50 
174.10 || 38...} 337.63 174.68 || 39....| 339.36 175.88 
40.... 217.82 |...........]] 209.33. 420...) 220.06 
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In this equation the approximate (B; — B,’) and (D, — D,’) were derived from the 
approximate constants B, and D, found in step 1. The quantities R»(J — 1) and 
P»(J) were found by averaging the wave numbers of lines with the same J. Thus every- 
thing on the left side of the equation was known. The differences (B; — B; ‘eorr and 
(Di — Di ‘eorr are corrections to be added to the approximate differences (B; — Bj’) 
and (D; — Di’) to give more accurate differences. Since the last term of the equation 
was small for small J, it was neglected. A plot of the left side of the equation versus /” 
produced a straight line for small J , from which 2v9 was found by extrapolating to 

J = 0. The origins found in this way are included in Table 11. 


TABLE 10 


COMPARISON OF COMBINATION DIFFERENCES OF THE LOWER STATE 


(1-2) Bann | (2-2) BAND (0-2) Swan BAND 


| | | 
| AsFi'(J) | | | AeF2'(J) AaFs'(J) | AaFi’(J) | A2Fs'(J) 
(Cm=) | (Cm=!) | (C (Cm-*) | (C (Cm-!) (Cm!) | (Cm=!) | 


73.65 | 79.23 : 73.72 | 
79.99 | 85.85 79.93 
86.32 | 91.98 86.23 
92.49 | 98.35 92.46 
98.87 | 104.64 ; 98.80 
104.91 111.08 105.05 
111.39 | 117.16 111.67 
117.66 | 123.61 | 129. 117.72 
123.91 | 129.93 | 124.01 
130.23 | 136.22 130.15 
136.47 | 142.43 56 | 136.55 
142.70 | 148.92 142.71 
149.09 | 155.04 149.14 
155.23 | 161.40 155.30 
| 161.59 | 167.60 161.57 
167.75 | 173.92 167.75 

| 174.10 | 180.11 174.06 

| 180.28 | 186.40 | 180.56 

186.53 192.55 186.48 
| 192.69 | 198.86 | 205. 192.70 

| 198.94 | 204.95 | ; 199.11 | 
4205.15 | 214-301 205.05 
211.36 | 217.35 | 59 | 211.37 


217.47 | 223.72 | 229.56 | 217.43 | 223.52 | 229.62 | 217.60 | 223.60 | 229.61 


3. With the band origins known, the previous equation could be rearranged as follows: 


Ry (J —1) +Pm (J) — 2 (BL —Bi’) 2 (Di = De) J? (J +1)? = 2n 
J? 


= 2 (By — Br’) corr — 2 (Dg — De’ corr (J +1)? 


A plot of the left side of the equation versus (J + 1)? gave a straight line, from which 
— Be’ )corr and — Dh’) corr were found from the intercept and slope. The cor- 
rected differences (B — Bi’) and (D; — Dz’) were then computed by adding these 
corrections to the approximate differences B — By’) and (D7 — D7’) found in step 1 
above. 


Be) 
85.61] 91.51]........] 85.29 | 91.87 
91.97 | 97.98 | 85.74 | 91.98) 98.06 
15.. 98.38 | 104.45 | 92.16 | 98.39 | 104. 35 
104.66 | 110.82 | 98.84 104.67 110.78 
111.07 | 117.23 | 105.09 | 111.05 117.10 
117.31 | 123.55 | 111.40 | 117.32 | 123.46 
19... 123.63 | 129.80 | 117.65 | 123.68 | 129.73 
129.94 | 136.12 | 123.97 | 129.94 136.12 
136.67 | 142.43 | 130.21 | 136.31 | 142.35 
142.56 | 148.74 | 136.54 | 142.56 | 148.72 
148.84 | 155.09 | 142.76 | 148.73 | 154.90 
155.10 | 161.34 | 149.09 | 155.08 | 161.29 
161.54 | 167.51 | 155.28 | 161.43 | 167.43 
167.66 | 173.82 | 161.64 | 167.54 173.80 
173.96 | 180.01 | 167.81 | 173.90 179.99 
= ee 180.06 | 186.34 | 174.06 | 180.04 | 186.06 
oy 186.35 | 192.45 | 180.45 | 186.33 | 192.47 ae 
192.56 | 198.71 | 186.53 | 192.27 | 198.73 
a. 3... 198.86 | 204.82 | 192.50 | 198.87 | 204.90 B34 
32... 204.90 | 211.18 | 199.12 | 204.92 | 211.17 
211.25 |.217.23 | 205.14 | 211.24 217.26 
Ge 34... P1738 | 293 54 | 223s | 
| 
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From the differences (B; — B;’) and (D/ — D;’), which were derived in this way for 
each band, it was possible to find improved constants, B, and D,, for the upper and 
lower vibrational levels that were consistent with these differences and which also gave 
the best agreement with the constants B;’ and D,’ of the lower state, as found previously* 
from a study of the Swan bands. Tables 12 and 13 give the resulting constants. For com- 
parison purposes, the constants of the lower state resulting from the analysis of the 
Swan bands are included in Table 13. 

The analysis of the (0-2) band gave 1.1801 cm~ for Bj. This was in good agreement 
with the value of 1.1804 cm~ found by Fox and Herzberg.! However, since their de- 


TABLE 12 
ROTATIONAL CONSTANTS OF THE B'Il, STATE 


Computed 
Dy+ 106 
(Cm~') 


6.4 


aB(Cm-1) | Boomp(Cm-! 
Bo(Cm>!) AB(Cm"!) ‘eomp(Cm~') 108) 


| 


1. 1803 —0.3 


0.0241 
6.7 


7.0 
733 
7.6 


1.1562 
1.1311 
.0239 


0.0234 


1.1072 


1.0838 +0.0005 


| 
| 
| 


TABLE 13 
ROTATIONAL CONSTANTS OF THE XII, STATE 


D,+ 106(Cm=~!) 


B,(Cm~!) 
Swan System 


Dy* 10®(Cm~!) 


Swan System 


1.6068 


6.6 


1.5909 6.6 
1.5740 6.6 


termination was based upon the analysis of two well-observed bands and two incom- 
pletely observed bands, their value has undoubtedly greater weight. Accordingly, 
1.1803 cm was adopted as the most probable value of the constant Bj. 

The third column of Table 12 gives the differences between successive rotational con- 
stants of the upper state. These differences are not constant but have a tendency to de- 
crease with increasing v. However, there is a pronounced scatter of the differences, which 
is greater than the probable error of the determinations. This scatter makes it impossible 
to fit an exact curve to the data. The best that can be done is to represent the observed 
constants by the usual linear equation, 

B’ = —al(v'+3). 

By graphical methods it was found that the constant B; is 1.1922 cm and a; is 
(0.0242 cm. The computed B;’s given by this equation are included in the fourth column 
of Table 12 and the differences (O — C) in the fifth column. 

Similarly, a linear approximation to the constants D} gave the equation 


Di = [6.3+0.29 (0'+4)] X10-%cm™. 


on 
as 
f 
is 
| 
Pit | 
2 
1.1317 | — .0006 | 7.0 0 
3. — | 1.1075 — 7.5 +2 
| 
| | —0.5 
} 
2 
| 
| 
j 
| 
i 
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The computed Dj’s are given in the seventh column of Table 12 and the differences 
(O — C) in the eighth column. For the sake of convenience the molecular constants in 
these equations have been collected in Table 14. 

The method of Budd‘ was used to investigate the triplet splitting of the vibrational 
levels of the upper state. For *II — *II transitions, Budé has shown that the coupling 
constants, A, could be found from the combination differences. The scatter of the final 
results was so great that numerical figures cannot be given, but it can be said, in general, 
that the triplet splitting of the BII, state is much smaller than that of the ground state. 


b) VIBRATIONAL ANALYSIS 


The determination of the band origins, vo, of the nine new bands is described in the 
previous section. They are included in Table 11, which also includes the origins and heads 
of the four bands, (0-3), (0-4), (0-5), and (0-6), found by Fox and Herzberg. The table 
also gives the vibrational quanta as horizontal and vertical differences. The vibrational 
quanta presented in Table 11 show some large discrepancies. For instance, there are two 
measured values of AG;/., namely, 1594.63 and 1594.92 cm. A study of the observation- 
al data showed that this discrepancy arose through an uncertainty in the origin of the 
(3-2) band. The head of this band is blended with the strong Hg line at 2536 A and its 


TABLE 14 
MOLECULAR CONSTANTS OF THE Bll, STATE 


. 1.1803 39806 .46 
0.0242 cm™ 1106.56 cm“! 
. +0.29X10~* cm™ 


acco mpanying system of ghosts. Therefore, 1594.63 cm~ is to be preferred as the value 


of AG;/.. This compares favorably with the value 1594.57 cm= derived from the Swan 
bands, For the same reason, 933.37 cm™ is to be preferred over 933.08 cm~ as the value 
of AG; 5/2° 

A similar discrepancy is observed in the case of the quantum AG;/,. The value given 
by the Swan bands, 1571.38 cm™, is in close agreement with two of the three determina- 
tions, while the value 1571.64 cm, found as the difference between the origins of the 
(0-2) and (0-3) bands, is much too large. A study of the observational data shows that 
the internal uncertainties in the determinations of the origins of these two bands are of 
the order of 0.02 cm. The only explanation of the discrepancy seems to be that there is 
a systematic difference between the measurements made in the present study and in that 
of Fox and Herzberg. Since this possibility exists, no reliance can be placed upon the 
value of 1036.83 cm~ for AG{/. found as the difference between the origins of the (1-3) 
and the (0-3) bands. Accordingly, the value of 1036.58 cm~ is to be preferred. The vibra- 
tional quanta finally chosen as being most probably correct are collected in Table 15 
and plotted against vibrational number in Figure 2. 

The relative numbers assigned to the vibrational levels of the upper state must be 
correct, otherwise the AG-versus-v-curve of Figure 2 would not be smooth. However, the 
absolute vibrational numbering may be open to question. The locations of band origins 
in the Deslandres table presented in Table 11 show that the corresponding Condon parab- 
ola is very open, as, indeed, might be anticipated from the great difference between the 
internuclear distances of the upper and lower states. Therefore, if the vibrational num- 
bering of the upper state were too low by gne unit, the additional level would be observed 
only through transitions from it to lower levels with o> 3, say. By extrapolating from 


4A. Budd, Zs. f. Phys., 98, 437, 1936. 


i 
at 
Ke 
3 
B 
ve 
— 
4 
i 
id 
He 


88 JOHN G. PHILLIPS 


the AG-versus-v-curve, it is possible to predict quite accurately the wave lengths of the 
heads that these new bands would have. They were looked for very carefully, and, while 
the spectrum in that region is very complex, no features were found that would corre- 
spond to this prediction. Thus it is improbable that there is an error in the vibrational 
numbering of the upper state. 

TABLE 15 


VIBRATIONAL QUANTA IN THE UPPER AND LOWER STATES 


| 
Upper State Bll, Lower State 


| 
AGt+1/2= AGt'+1/2= AGo'+1/2 
G'(v+ 1) —G’(v) G’’(v+1) —G’’(v) Swan System 


1036.58 cm™ 
979.52 1594.63 cm™ | 1594.57 cm™ 
933.37 1571.36 1571.38 


ann 1 

2 3 4 Vv 
_ Fic. 2.—B, and AG curves of the upper (B*II,) state. Observed points on the B,-curve are indicated 
by closed circles, and points on the AG-curve are given as open circles. 


The differences between successive vibrational quanta of the upper state, given in the 
third column of Table 15, decrease systematically with increasing vibrational number. 
This is also shown by the upward curvature of the AG-versus-v-curve in Figure 2. An 
equation that would represent the observed quanta accurately would have to have as 
many terms as there are quanta. The resulting equation for the origins of the nine bands 
is as follows: 


vo = 40080.41 + 1106.56 +3) —39.260 ( 0’ + 2.8050 + 
— 0.1271 +4) 4— [1641.35 (0 +3) —11.67 +3) 
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The constants for the lower state were found previously from the study of the Swan 
bands. Any extrapolation made to higher vibrational levels on the basis of this equation 
would be very uncertain, since it is questionable whether the upward curvature of the 
AG-versus-v-curve is maintained for higher v. It is preferable to regard the equation mere- 
ly as a representation of the observed levels only. The constants of the upper state in the 
equation are included in Table 14. 

The upward curvature of the AG-versus-v-curve of the B*II, state is quite anomalous. 
Usually the curve is either linear or shows a downward curvature. Similar upward curva- 
tures have been observed in the case of excited states of a few molecules. An example is 
AgH. There have, in general, been two explanations given for these upward curvatures. 

1. The excited state in question is being perturbed by another electronic state. If this 
is the correct explanation in the present case, the perturbing state must have a lower 
equilibrium energy than the B*Il, state in order to produce the observed perturbation. 
Furthermore, such a vibrational perturbation can occur only if the Kronig selection rules 
are fulfilled and also if AA of the two perturbing states is zero. Thus the perturbing state 
must be a II, state, and its presence would be known through the observation of the 
band system resulting from transitions to the X*II, state. The only observed state ful- 
filling these conditions is the A*II, state, which is the upper state of the Swan system. 
Therefore, if the anomaly observed in the case of the B*II, state is due to a vibrational 
perturbation, then there exists a strong possibility that the A*II, state is the perturbing 
state. This possibility is strengthened somewhat when the equilibrium internuclear dis- 
tances of the two states are considered. The equilibrium distance of the A*II, state is 
1.2660 A, while that of the B*Il, state is 1.5350 A. Thus the point of closest approach of 
the corresponding potential curves would occur for low vibrational levels of the B*IT, 
state, as is actually observed. 

2. A second explanation that has been given for the anomalous behavior of the excited 
states of some molecules is that ionic forces between the two atoms play a role.® It is 
difficult to decide whether or not this effect would be important in the case of the C; 
molecule. 

At the present time no decision can be made as to which of the two explanations 
given above is the correct one in the case of the B*II, state. Such a decision must await 
the discovery and analysis of bands from much higher vibrational levels of this state 
and an accurate knowledge of the dissociation energies. 


I wish to acknowledge many helpful suggestions made by Dr. G. Herzberg during the 
course of this study and also his generosity in reading the manuscript. 


5 See, e.g., A. G. Gaydon, Dissociation Energies and Spectra of Diatomic Molecules (London: Chap- 
man & Hall, Ltd., 1947). 


6A. G. Gaydon, Proc. Phys. Soc., 58, 525, 1946. 
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THE PROBABLE BEHAVIOR OF THE NEXT SUNSPOT CYCLE 


W. GLEISSBERG 
Istanbul University Observatory 
Received A pril 18, 1949 


ABSTRACT 


The probability laws of sunspot variations, which have yielded successful predictions for the present 
sunspot cycle, lead to the conclusion that, with a high degree of probability, a steep ascent to a very high 
maximum is to be expected also in the next cycle and that the period of low activity preceding the next 
cycle probably will be extremely short. 


In the September, 1942, issue of this Journal' I published four probability laws of sun- 
spot variations, which permit us to predict the probable features of a sunspot cycle from 
the observed features of the preceding ones. From these probability laws I deduced the 
following three predictions for the then next sunspot cycle, the eighteenth in the Ziirich 
statistics: (1) The probability that the highest smoothed relative number of the coming 
cycle will exceed 145 is 86 per cent. (2) The probability that in the coming cycle the re- 
duced period of rising will be shorter than 32 months is 91 per cent. (3) It may be ex- 
pected with a probability of 93 per cent that the period of low activity which will precede 
the next spot cycle will be shorter than 40 months. 

Now the maximum of cycle 18 has passed. It occurred in May, 1947, and the highest 
smoothed relative number was 151.8. The reduced period of rising, i.e., the time during 
which the smoothed relative sunspot numbers increase from a quarter of their maximum 
value to the maximum, began in August, 1945, and ended in May, 1947; it therefore 
lasted 21 months. And as, in cycle 17, the reduced period of falling had ended? in July, 
1942, so the period of low activity preceding cycle 18 lasted from July, 1942, to August, 
1945, i.e., 37 months. 

From these observational data we see that all the above predictions have come true. 
As far as I know, this is the first time that essential features of a sunspot cycle have been 
successfully predicted before its beginning. Moreover, the success of these predictions 
proves that the hypothesis of the occurrence of alternately high and low maxima cannot 
be maintained; for, according to this hypothesis, one would have expected a low maxi- 
mum in cycle 18. Actually, a prediction which was based upon this hypothesis® has failed. 

As the length of the period of low activity, the length of the reduced period of rising, 
and the height of the maximum now are known for the present sunspot cycle (i.e., cycle 
18), we are able to establish predictions concerning the corresponding characteristics of 
the next cycle (i.e., cycle 19). 

When we wish to predict the probable features of a solar cycle, we have first to in- 
vestigate whether ¢{*”, the four-cycle average of the reduced period of rising, is increasing 
or decreasing. For this purpose let us complete the table of the characteristics of the sun- 
spot cycles published in my earlier paper cited above.! By adding, to Table 1 of that 
paper, the data yielded by the observations in the meantime, we obtain the accompany- 
ing Table 1. In this table Rsz denotes the highest smoothed relative number, ¢; the period 
of low activity, ¢, the reduced period of rising, ¢; the reduced period of falling, and the 
quantities marked with the upper index (4) are the averages for every four consecutive 


1 Ap. J., 96, 234, 1942. 
2 W. Gleissberg, Beob.-Zirk., 25, 89, 1943. 
3 A. H. Shapley, Terr. Mag., 49, 43, 1944. 


ine 
> 
\ 
Bas 
: 
| 
| j 
H 
3 
q 
{ 
90 


91 


BEHAVIOR OF NEXT SUNSPOT CYCLE 


cycles. The unit of ¢,, ¢,, and ¢; is 1 month. As the arrangement of this table is the same as 
that of Table 1 in my previous paper,' no further explanation of its contents is needed. 

When I was making the predictions for cycle 18, the last known value of ¢{*) was 33, 
and there was no doubt that ¢{*) would continue decreasing. The present situation is 


TABLE 1 


THE CHARACTERISTICS OF THE SUNSPOT CYCLES 
AND THEIR FOUR-CYCLE AVERAGES 


Cycle No. Ry t, RiP 
1 86.5 41 52 48 
49 
2 115.8 22 29 58 
125.5 33 31 59 
3 | 158.5 29 17 48 
| 116.2 40 33 54 
ee 141.2 41 28 81 
99.4 52 34 53 
5 Pores 66 57 27 
77.7 59 43 49 
6 48.7 73 35 54 
79.1 56 42 42 
99.7 51 36 52 
112.0 44 35 56 
131.6 45 34 68 
129.2 37 29 60 
97.9 45 31 70 
111.2 45 34 56 
140.5 30 26 51 
100.2 47 35 53 
44 74.6 62 47 34 
91.8 49 36 49 
87.9 51 36 57 
83.0 55 38 48 
14 64.2 53 35 54 
83.9 49 37 49 
15 105.4 52 32 46 
91.7 47 33 51 
16 78.1 39 43 40 
113.6 43 30 
119.2 45 23 63 
18 | 151.8 37 21 


more difficult. Now the last known value of ¢*) is 30; and, as the two minima of ¢{*) which 
have been observed hitherto have been 31 and 29, respectively, it seems uncertain 
whether the next value of ¢{) will be below 30 or above 30. We can, however, base our 
judgment upon the course of the values in columns Rj? and ¢{*), which are strongly 
correlated to the values of ¢{*. From these columns we see that the last known values of 
Ri and ¢{* still are far from their previous maximum and minimum values, respectively. 
Thus we may su pose that the increase of R{}) and the decrease of ¢{*) will continue and 


that, therefore, ¢{*) is still decreasing. 
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On the supposition that ¢{“) is still decreasing, the following equations can be deduced‘ 
from the probability laws of sunspot variations: 


prob —1] =0.95, (1) 
prob[R,, > 1080 — —S,(R,)] = 0.95, 2) 
prob[t, < — S, (t,) — 150] = 0.95, (3) 
where the symbol “prob” denotes the probability that the inequality which follows this 
symbol in brackets will be fulfilled in the next sunspot cycle; S; denotes the sum of the 


quantity in parentheses for the last three cycles; and /‘*)* denotes the last known value 
of t{*. From Table 1 we have 


*=30, =214+23443 =87, 
= 151.84119.24 78.1 = 349.1, 
S3(t) =37445439=121. 


Substituting these values in equations (1), (2), and (3), we obtain 
prob[t, <32] =0.95, prob [Ry > 130.9] =0.95 , prob[t,; < 29] =0.95 . 


These equations contain the following predictions: 

It can be expected with a probability of 0.95 that: (1) in cycle 19 the smoothed rela- 
tive sunspot numbers will rise from a quarter of their maximum value to the maximum 
within less than 32 months; (2) the highest smoothed relative number of cycle 19 will 
be higher than 130; and (3) the period of low activity which will precede cycle 19 will 
be shorter than 29 months. 

Predictions 1 and 2 indicate that in the next sunspot cycle a steep ascent of solar ac- 
tivity to a very high maximum will again occur. Prediction 3 is of special interest, since 
a period of low activity which was shorter than 29 months has not been observed since 
cycle 2, i.e., for more than 180 years (see Table 1, column ¢;). From this short duration of 
the period of low activity we can conclude that the coming minimum will not be a deep 


one. 


4 A detailed proof of these equations will be published in Pub. Istanbul U. Obs., No. 36. 
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NOTES 


THE DIAMETER OF NEPTUNE* 


The American Ephemeris and Nautical Almanac for 1950 uses 2744 for the diameter 
of Neptune at its mean distance, 30.07 a.u. In the revised edition of Russell, Dugan, and 
Stewart’s Astronomy, Vol. 1, the value 2"28, due to W. Rabe,! has been adopted. Rabe’s 
value is based on seven series of observations made over a century after personal equa- 
tions estimated by Rabe were applied to the original measures. All but one series (by 
Lassell) were based on micrometer measures. 

Experiments carried out with the 82-inch telescope have convinced the writer that 
micrometer measures on small disks are far less trustworthy than measures obtained by 
either of the following two types of instruments: (1) a disk meter, showing an artificial 
disk matching the appearance of the planetary disk in diameter, brightness, and color; 
and (2) a double-image micrometer in which the two images produced are brought into 
contact. The writer used an instrument of the first type in measuring satellite diameters,” 
following earlier similar work by Camichel; while an instrument of the second type was 
acquired from Dr. B. Lyot in 1946. 

The diameter of Neptune was measured with the 82-inch telescope on four occasions, 
as shown in the accompanying tabulation. The magnification in all cases was 720 diame- 


Date (U.T.) Instrument Seei d” (Mean distance) 


March 1, 1948 Disk meter 2.079 
March 24, 1948 Disk meter 2.041 
Feb. Double-image 2.013 
Feb. 24, 1949...... Disk meter 2.054 


Weighted average 2.044+0.006 (m.e.) 


ters. The measures with the double-image micrometer were somewhat difficult because 
of the light-losses inherent in this irstrument and the comparative faintness of Neptune. 
The disk-meter measures were very easy and were made in series of several settings each, 
separated by other activities. The number of series is given in the table in parentheses. 
No pronounced ellipticity of the disk was noted, and the diameter measures refer to the 
average for the visible disk. 

The scale value of the attachments was determined as follows. A preliminary calibra- 
tion was obtained from the dimensions of the disk meter, from its optics, and from micro- 
scopic measurement of the holes producing the artificial disks. This calibration is called 
the “provisional” one and is attached in tabular form to the instrument. The double- 
image micrometer was used mounted on the disk meter and tied in with it system 
through measurement of the largest of the artificial disks by the double-image method. 


* Contributions from the McDonald Observatory, University of Texas, No. 179. 


1A.N., 234, 153, 1928. 
2 The Atmospheres of the Earth and Planets (Chicago: University of Chicago Press, 1949), pp. 305-306; 
Contr. McDonald Obs., No. 161. 
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This could be done very accurately. Thereupon, three visual binaries, for which good 
photographic positions were available, were measured with the double-image microme- 
ter (a direct comparison of artificial disks with visual binaries was much less satisfactory, 
although it had led to a rough verification of the provisional scale). This calibration led 


| Measured 
37408 +0"040 
3.792 — .051 
4.113 +0.009 


m.e. of the average of the measured positions +07027, or +0.7 per cent. The correction 
to the provisional scale is therefore +0.5 + 0.7 per cent (m.e.), which, while small and 
uncertain, was applied to all measures quoted. 

Combining measuring errors and scale errors we have: 


d (Neptune) = 27044+07016 (m.e.) . 


It is believed that systematic errors are negligible. Accordingly, we adopt d = 2704 + 
0702; diameter 44,600 + 400 km = 3.50 + 0.03 Earth; and density equal to 2.22 cgs. 
The correction to the density is large and appreciably increases the diversity among the 
Jovian planets. 

The magnitude of Neptune was measured photoelectrically by A. Shatzel on May 3, 
1949, who found 7.95 IPv reduced to mean opposition (29.07 a.u.). This value is 0.30 
fainter than that previously assumed, but the albedo is changed only slightly, to 0.49, 
owing to a similar correction to the disk size. However, it has been suggested by W. 


Becker that Neptune is somewhat variable. 
GERARD P. KuIPER 


YERKES AND MCDONALD OBSERVATORIES 
July 12, 1949 


NOTE ON THE PHYSICAL SIGNIFICANCE OF STEWART’S 
SUNSPOT FORMULA 


Stewart and his collaborators! have given the following mathematical expression for 

the course of solar activity within a spot cycle: 

where R denotes Wolf’s relative sunspot number and 6 the time as reckoned from the 
beginning of the cycle; F, a, and 6 are constants whose values vary from cycle to cycle. 
The aim of this note is to show that the above formula, which has been established by 
Stewart and his collaborators without proof, can be deduced from very simple and rea- 
sonable suppositions concerning the rates of formation and disintegration of sunspots. 

Since the variation of solar activity arises from the formation of new spots and from 
the disintegration of existing spots, we can write 


dR 
fi (R, 0) — fo(R, 6) . 


1J. Q. Stewart and H. A. A. Panofsky, Ap. J., 88, 385, 1938; J. Q. Stewart and F. C. Eggleston, 
Ap. J., 91, 72, 1940. 
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where /fi(R, 6) and f2(R, @) represent the formation and the disintegration of spots, 
respectively. The functions /i(R, 6) and f2(R, @) must be positive for all possible values of 
R and 0. The higher the solar activity, the greater the number of new formed and of 
disintegrating spots. Thus we may suppose that fi(R, @) and f2(R, @) are proportional 


to R: 
fi(R, 0) =Rf; (6), fo(I, 0) =Rf,(6), 


where /3(@) and f4(@) again are positive functions. On the other hand, it is obvious that, 
while the rate of disintegration need not change with time, the intensity of formation of 
new spots is gradually decreasing during the spot cycle. Thus the most simple supposi- 
tion will be as follows: 

a 


= 0, 


f3(0) = 


where a and 6 are positive constants. Hence 


dR a 
OR. 


Integration of this differential equation yields 
R=F@e*, 


where F is the constant of integration. This is Stewart’s formula. Its physical sig- 
nificance, therefore, can be expressed as follows: The rate of formation of new spots and 
the rate of disintegration of existing spots both are proportional to the intensity of solar 
activity; and the rate of formation of new spots, moreover, is inversely proportional to 
the time from the beginning of the spot cycle. 

W. GLEISSBERG 


ISTANBUL UNIVERSITY OBSERVATORY 
April 9, 1949 


AN ULTRAVIOLET LIGHT-CURVE OF V 444 CYGNI* 


In the course of photometric observations of Wolf-Rayet spectroscopic binaries with 
the 82-inch reflector of the McDonald Observatory, a few observations in the ultraviolet 
were obtained of V 444 Cygni. This system is a Wolf-Rayet eclipsing binary in which the 
primary minimum (Wolf-Rayet star in front) is appreciably broader than the second- 
ary.! The difference in width cannot be explained by eccentricity of the orbit but must 
have its origin in a tenuous scattering envelope about the Wolf-Rayet star. If this 
envelope has a high electron density and is responsible for scattering of light from the 
secondary star at primary minimum,’ the width of the primary minimum should be 
independent of wave length. This has been investigated by Kron and Gordon’ at effective 
wave lengths \ 4500 and \ 7200. Their observations showed that, after correcting for a 
slight difference in depth, there was no significant difference between the two wave 
lengths. The observations presented here give further evidence that the character of the 
primary minimum is independent of wave length. 

The photometer employed was that previously described. A UG2 filter and 1P21 


* Contributions from the McDonald Observatory, University of Texas, No. 170. 
1G. E. Kron and Katherine C. Gordon, Ap. J., 97, 311, 1943. 


2 Zdenék Kopal and Martha B. Shapley, Ap. /., 104, 160, 1946. 
3 A.J., 154, 1947. 4 Ap. J., 108, 58, 1948. 
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photocell were employed for isolating the ultraviolet. The effective wave length was not 
carefully determined but is estimated to be \ 3550 A from the known transmission of 
the filter and published spectral sensitivity of the 1P21 photocell. 

The individual observations are recorded in Table 1. They are presented in detail in 
order to facilitate future studies of the time of primary minima. The phases were com- 
puted from the formula 


JD = 2428771.379+ 4.21238. 


The individual observations at primary minima have been plotted in Figure 1, along 
with the light-curves in the blue and infrared as determined by Kron and Gordon, so 


TABLE 1 
PHOTOMETRIC OBSERVATIONS OF HD 193576 AT A 3550 


| 
| | 


Heliocentric Am | | Heliocentric Am 
Phase | HD 193032— Phase HD 193032- 
(In Days) V 444 Cyg \| ‘ | (In Days) | V 444 Cyg 


2.968 +0.006 
2.971 | .009 
2.973 
3.519 
3.500 
0.024 
0.027 
0.036 
0.038 
0.042 
0.111 
0.115 
0.117 


ooo 
S38 


4 

JD | | on 
2432+ | 
} 
430.693........| || 448.664........) 4.0900 | 
430.696........| 448.673........| 4.099 | 195 
430.698........| || 448.676........; 4.102 | 201 
439.669........| || 448.678........) 4.104 | 
439.672........| || 448.736........) 4.162 | 269 
444,598........| 448.742........| 4.168 | 265 
444.610........| | 449.588........| 0.801 | 016 
444.612........| 449.590... 0.803 013 
444.616........| 449.594........| 0.807. | 014 
444.685........| || 449.598........) O.811 | 013 
ya 444.689....... | 451.578........| 2.791 | 012 
444.691........| 451.581........| 2.704 | 015 
| 222 451.583........| 2.796 | 016 

444.777........| 125 ‘|| 451.587........| 2.800 | 004 
444.780........| 451.589........| 2.802 O11 
444.783........1 0.209 | 114 ‘||: 452.580... | 021 
445.590........| 1.016 | 000 452.585........] 3.798 | 021 
445.594........} 1.020 | ) 5) 
445.600........) 1.026 | 

aS 1.149 | 

445.734......,.] 1.160 

447.586........| 3.012 

447.589........| 3.015 
447.591........| 3.017 
447.594........| 3.020 
4.006 

448.587........| 4.013 fas 
448.500........| 4.016 
448.658........| 4.084 
448.661........} 4.087 | | | | 


Phose in days 


Fic. 1.—Light-curves of V 444 Cygni. The observations at \ 3550 were obtained by Hiltner and 
those at \ 4500 and \ 7200 by Kron and Gordon. 
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that comparison of the minimum at the three different wave lengths is facilitated.® It is 
obvious that, after a slight correction for difference in depth, the character of the ultra- 
violet-curve is essentially identical with the others. Consequently, this further supports 
the conclusion that the scattering in the envelope is a consequence of high electron 
density. 

It had been observed by Kron and Gordon’ that V 444 Cygni was subject to night- 
to-night fluctuations by as much as 0.1 mag. in the infrared and 0.06 mag. in the blue. 
Similar fluctuations, but smaller, have been noted by the author in CQ Cephei,® another 
eclipsing binary in which one component is a Wolf-Rayet star. This point was given care- 
ful attention in the present investigation. The comparison star HD 193032 was checked 
against a second comparison star, HD 193681. The results are given in Table 2. On the 
assumption that there is no variation of the two comparison stars, the probable error of 
a single observation is + 0.0037 mag., which is not abnormal. It is apparent from Figure 1 
and from the unplotted points of Table 1 that the nightly fluctuations were never greater 
than 0.01 mag. That no large fluctuations were observed here may have one of two 
explanations. Either the fluctuations are a function of wave length, or, which is more 


TABLE 2 


COMPARISON OF HD 193032 AND HD 193681 
Am Am 
JD 2432+ HD 193681—HD 193032 JD 2432+ HD 193681—HD 193032 
+0.029 449.591... +0.029 
439.676...... 449 596. . 
444.617... , .035 028 
444.695........ .033 451.588... 
+ 019 452.589... 
.033 455.588 035 
448.591... + .035 456.582 + .036 
+0.021 459.591... +0.033 


probable, the present observations, which were continued for only one month, were made 


during a period of inactivity. 
W. A. HILTNER 


YERKES AND MCDONALD OBSERVATORIES 
March 23, 1949 
® The author is grateful to Dr. Kron and Miss Gordon for permission to reproduce their observations 
in the blue and infrared. 


A.J., 53, 198, 1948. 
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REVIEWS 


Einfiihrung in die Dynamik von Sternsystemen. By E. VAN DER PAHLEN. Basel: Verlag Birk- 
hiuser, 1947. Pp. 240. Swiss Fr. 36. 


Investigations on the mathematical theory of stellar dynamics have centered largely on two 
problems. The first relates to stellar encounter, its importance for stellar dynamics and the 
manner in which it is incorporated into the framework of a general theory. The second relates 
to the consequences for the dynamical theory which the known kinematics of stellar motions 
implies. The present book is devoted to summarizing and in part extending the results of known 
investigations on these two problems. 

In the treatment of stellar encounters the author follows an early investigation of Charlier. 
This method consists in evaluating the energy, AE, exchanged between two stars describing 
hyperbolic orbits about each other and finding the average rate of increase of AE?. But in the 
treatments which have followed Charlier’s, AE is first averaged over some of the parameters of 
the encounter, squared, and then averaged over the remaining parameters. The reviewer cannot 
see that this procedure gives any meaningful result. It would seem far more rational to analyze 
the effect of encounters on the motion of a star directly by finding the mean rate of change of the 
velocity in the direction of motion and at right angles to the direction of motion. When this 
is done, it is found that the effect of encounters is systematically to reduce the mean velocity 
to zero while at the same time making the mean dispersion tend to a finite limit. Treatments of 
stellar encounters along these lines appeared during the war, but they were apparently inac- 
cessible to the author. 

‘The treatment of the second major problem in the book is exceptionally complete and satis- 
factory. As is well known, the central problem here is equivalent to one of defining the potential 
in a dynamical system which admits an integral (other than the energy integral) which is quad- 
ratic in the velocities. Detailed investigations of this problem have appeared in the literature, 
but Dr. van der Pahlen’s treatment is more concise and in certain parts more elegant than the 
published investigations. 

The book is attractively written and is a very useful addition to astronomical literature. 


S. CHANDRASEKHAR 
Yerkes Observatory 


Kinematic Relativity. By E. A. Mine. Oxford: Clarendon Press, 1949. Pp. vi+-238. $6.50. 


One of the most fully documented alternative theories to the Einstein theory of gravitation 
which has appeared in recent years is the “kinematic relativity” theory of Professor E. A. Milne, 
of Oxford University. The book under review is the second written by him on this subject; the 
first, entitled Relativity, Gravitation, and World-Structure, was reviewed in an earlier issue of this 
Journal.’ The basic formulation of the theory remains unchanged, but a considerable extension 
has been made in the discussion of the gravitational problem and in the extension to electrody- 
namics. The exposition of the preliminary kinematical part of the theory has been materially 
altered from that of the earlier book by the introduction of a group theoretical discussion along 
the lines of the paper by Milne and Whitrow.? The reader will find it profitable to refer also to 
the general treatment by Robertson® in this Journal, to which, strangely enough, no reference is 
made in this book. In consideration of the extensive discussion of the kinematical aspects‘ of 


1 Ap. J., 83, 61, 1936 (reviewed by H. P. Robertson). 
2 Zs. f. Ap., 15, 263, 1938. 
3 4 p. J., 82, 284, 1935; 83, 187 and 257, 1936. 


4 The reader may be interested to refer to the older discussion by G. D. Birkhoff in his book, Rela- 
tivity and Modern Physics (Cambridge, Mass.: Harvard University Press, 1927), which proceeds from 
arguments very similar to those developed by Milne. 
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Milne’s theory which has already appeared in the literature, we shall not undertake to summa- 
rize it again here but will aim rather at giving the reader unfamiliar with Milne’s ideas some 
notion of their scope and intent in relation to the more general corpus of present-day theoretical 
physics. 

It is perhaps not unfair to say that in Milne’s view the most fundamental fact about the uni- 
verse is that it is “expanding” and that this fact is of such transcendent importance that we must 
build all our theories of physics to be conformable with it rather than consider the expansion to 
be derivable from other more prosaic laws of physics. We are challenged to develop all our theo- 
ries of dynamics, gravitation, electrodynamics, cosmic rays, electron spin, etc., on the fixed pat- 
tern provided by the theory. It is this contention which makes of Milne’s theory much more 
than just a model for the expanding universe and which requires that its fundamental postulates 
be subjected to full exploration from all angles before assent is accorded to them. In a word, 
Milne proposes to give us a theory comparable in scope to that originally proposed for Einstein’s 
general theory of relativity as a comprehensive theory of the physical universe, with the as yet 
missing correlation with quantum mechanics presumably to be filled in later. 

The basic element of Milne’s universe, which is called by him the substratum, consists of a 
distribution of moving points, or ‘fundamental observers,” concerning which three basic as- 
sumptions are made: (a) each fundamental observer moves with constant velocity in a radial di- 
rection directly away from any other fundamenta! observer, as observed in the reference system 
of the latter; (6) the Hubble speed-distance law of nebular recession, giving a linear homogeneous 
relation between the recession speed and the distance between any two observers, is exactly 
correct out to the confines of the observable universe, which is then bounded by observers hav- 
ing a recession speed c (velocity of light); and (c) the distribution of the fundamental observers 
in space is preserved under the only continuous group of symmetry transformations between 
fundamental observers which is allowed in the theory, which is the 6-parameter homogeneous 
Lorentz group of transformations (i.e., not including changes of origin of the space and time 
co-ordinates). The theory necessarily exhibits a singular epoch of ‘‘creation,”’ which on the time 
scale of this expanding universe is at ¢ = 0, the same for every fundamental observer; at this 
instant the particles of this model universe began to pour out from the common point of creation 
x = y =z =0. Milne shows that these conditions can be satisfied in just one way. 

The manner in which the concept of gravitation comes into the theory can be indicated 
briefly by a consideration of the argument used by Milne for the establishment of the equation 
of motion for a free test particle in the presence of the substratum. The equation itself is 


0 (1) 


v2 
’ 


c 


dv 
ar (r — vt) 


where r = xi + yj + sk is the space co-ordinate vector of the particle and v = dr/dt is 
its velocity vector. The singularity associated with the epoch of creation requires that the solu- 
tions of equation (1) satisfy the conditions 


Ir|<ct, |vl<e, as 


Equation (1) is actually established by Milne with no thought of gravitational action be- 
tween the test particle and the substratum in mind, but simply from the requirement that it 
shall be as similar as possible to Newton’s equation of motion of a free particle and that it shall 
be form-invariant under the allowed group of symmetry transformations between equivalent 
fundamental observers. In fact, Milne regards his arguments as a derivation of the dynamical 
equation of motion (1) from purely kinematical principles. But if we now regard equation (1) as 
a Newtonian equation of motion for the particle, then the second member would be considered 
as being of the nature of a force exerted by the substratum on the particle; and we therefore 
attribute it to the gravitational action of the substratum on the particle! It must be said in fair- 
ness that Milne makes a much more extended analysis of the analytical character of this assump- 
tion by the consideration of more general statistical distributions of fundamental observers and 
of the corresponding equations of motion for a particle; but the reviewer believes that the state- 
ment given above is a fair summary of the essential physical argument by which Milne converts 
his kinematic assumptions into a theory of gravitation. It constitutes Milne’s substitute for 
Einstein’s “principle of equivalence.” 

One can but be struck at once by the fact that, since in Milne’s argument the test particle 
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and the substratum (or more general statistical system) form an indissoluble entity and since 
there would presumably never arise the concept of a particle moving freely in the presence of 
nothing else, we should be at a loss to know what the Newtonian equation dv/dt = 0 could pos- 
sibly mean, despite the fact that it is evidently being used as the basis for the interpretation of 
the second member of equation (1) as a gravitational force. It would seem to be much more 
consonant with the logic of Milne’s argument to regard equation (1) as the natural equation of 
motion of the particle under no forces; but then the whole connection between kinematics and 
gravitation would appear to be dissolved, or at least pushed onto some other plane. In any event 
the reviewer will readily admit that Milne’s arguments on this point seem to him to be very pre- 
carious as a physical basis for such an important element of the theory. 

One of the further points in Milne’s theory which has been much discussed is the introduction 
of a second time scale (r-scale). If one makes the following transformation of variables on the 
equations of motion of a particle as given by Milne 


r=Rich «sinh —, 
0 


where R= x'j + y'i + 2’k and R, = R/R, then to a first approximation for low velocities these 

_ equations can be reduced to Newtonian form (sec. 79). From this Milne contends that we must 
consider two time scales— (a) the ¢-scale appropriate to the description of ‘expanding universe” 
phenomena and (6) the r-scale appropriate to comparisons with Newtonian dynamics. In this 
view the time variables appearing in Newton’s equations and in the Maxwell-Lorentz equations 
of the electromagnetic field are of quite different natures. This is perhaps not a question to which 
there is a quick and easy answer, but the following remark seems to be in order. Since the struc- 
ture of the underlying symmetry group of the theory cannot be altered merely by making a 
transformation of co-ordinates, it hardly appears feasible to argue that the 6-parameter homo- 
geneous Lorentz group of the ¢-scale theory is transformable in this manner into the 10-param- 
eter symmetry group of the Newtonian theory. The proposed reduction to Newtonian mechan- 
ics would therefore be expected to be incomplete in the sense that, to maintain the correlation 
with the ¢-scale theory, we could not expect to make use of the full symmetry group of Newton’s 
equations in the ¢-scale theory. 

Turning to the introduction of a theory of the electromagnetic field within the framework of 
Milne’s ideas, one might anticipate that, since the (homogeneous) Lorentz group is a part of the 
theory from the beginning, there would be no particular difficulty until he tried to establish con- 
tact with dynamics through the expression for the ponderomotive forces of the fields on the 
charges and currents, with the resultant determination of the physical meaning of the field vec- 
tors. This seems to be the case, in fact; but, unfortunately, the argument is too detailed to sum- 
marize here in an intelligible manner, and the reader must refer to Milne’s treatment. 

One’s reaction to Milne’s theory as a whole is almost certain to be determined by his reaction 
to Milne’s fundamental philosophy. It is the reviewer’s opinion that we may willingly allow the 
logical possibility of constructing a self-consistent model universe along the lines assumed and 
may even grant that features will arise in it which may be made to correspond with more or less 
exactitude to such familiar ideas as gravitation, electricity, etc. But it is clear from the mathe- 
matical structure of the theory that these correspondences must be forced in nature, owing to 
the restricted character of the underlying symmetry group allowed in the theory. It does not 
seem at all surprising, for instance, that the attempt to formulate the dynamical equations in 
Lagrangian and Hamiltonian form remains incomplete. Quite apart from the comparison of 
Milne’s and Einstein’s theories, there remains the far more important fact that neither of these 
theories seems to offer a very hopeful basis as a framework for a generalized formulation of 
quantum mechanics. The entirely different mathematical structure of modern quantum theory 
makes it more likely that they will prove to be incompatible, except in some limiting sense. If 
this view is correct, then it is probable that Milne’s theory has been introduced too late into the 
scheme of modern physics to be of much influence on the development of quantum-mechanical 
theory. Its influence in other directions remains to be determined. — 

. L. Hit 


Department of Physics 
University of Minnesota 
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The Atmospheres of the Earth and Planets: Papers Presented at the Fiftieth Anniversary Symposi- 
um of the Yerkes Observatory, September, 1947. Edited by GERARD P. Kuiper. Chicago: Uni- 
versity of Chicago Press, 1949. Pp. viit+366. $7.50. 


This volume is an impressive illustration of the extent to which modern physical methods of 
study have expanded our knowledge of what, not so many years ago, seemed to be small regions 
of study. It is confined to those branches of the field which are definitely of astronomical interest 
—and therefore deals only incidentally, as required for its purpose, with the main mass of studies 
of the ionosphere. Nevertheless, it contains so much important material that a full study, espe- 
cially of the admirable research surveys included in it, would occupy a large part of a year’s work 
for an energetic graduate student. 

The strongest impression produced by reading it is of the extreme complexity of most of the 
problems. This is astrophysics at room temperature—or at least at laboratory temperature— 
where the properties of matter are complicated by chemical combination, condensation of liq- 
uid, and solid particles, etc., and the physical conditions undergo large diurnal and seasonal 
variations. Eddington wrote wisely of the hope of understanding “‘so simple a thing as a star.” 

No one reviewer could speak with authority upon all the problems discussed in this volume. 
The writer’s qualification is that he was unable to attend the symposium, and so can at least 
view the work from the outside—so far outside that Rossby’s important meteorologica! paper 
can only be mentioned. 

The thorough discussion of atmospheric scattering by van de Hulst is a full-dress research 
survey. He regards photometric observations during twilight as the most promising method for 
— study of scattering in the earth’s atmosphere and calls for more observations upon the 
planets. 

The observations of the ultraviolet solar spectrum from rockets— doubtless the most striking 
feature of the original symposium—are now familiar, but Greenstein’s estimate of the intensity 
to be anticipated in the ‘‘continuous” background, though published in 1943,! is less so than it 
should be. The rapid drop in this occurs at a longer wave length than he then anticipated, prob- 
ably owing to absorption by QO, in the solar atmosphere, as he later suggested.? Beyond A 900, 
‘ta direct photograph would show the corona and not the sun.” In chapter v, Whipple, Jacchia, 
and Kopal present evidence for seasonal variations as the density of the upper atmosphere, de- 
rived from Harvard observations of meteors; but they wisely refrain from theoretical explana- 
tion. 
Chapter vi contains a second research survey by Swings on the spectra of the night sky and 
the aurora. Previously proposed identifications are critically examined. Those accepted for the 
night sky are the first positive and Vegard-Kaplan bands of Nz and the Herzberg system of O2 
and the familiar lines of [0 1] and Na. The conditions of their appearance and the evidence for 
other possible bands are discussed in detail, and also the twilight phenomena, which include the 
appearance of bands of N2*. All these appear in the auroral spectrum as well as lines of [NV 1] and 
at times the Balmer lines and lines of He I or even He 11. A thorough discussion of the conditions 
for excitation of both spectra ends with the sentence: ‘“‘We are on the opening stages of the ex- 
ploration of one of the most fascinating fields in geophysics.” In the following chapter Spitzer 
treats of the extreme upper part of the earth’s atmosphere. The ‘“‘exosphere,” where the free 
path is so long that upward-moving molecules are brought back by gravity rather than by colli- 
sion, is reached at a density not far from 10-? that at the surface. For the corresponding altitude 
‘a value between 500 and 1000 km seems likely.” The uncertain factor is the temperature at 
the high levels, which depends upon the intensity of solar radiation beyond 900 A. A new calcu- 
lation of the rate of escape of gases from an atmosphere much hotter at the top than lower down 
indicates that a steady temperature of about 1500° C is required to dissipate helium into space 
as fast as it is released by the weathering of igneous rocks, but one as low as 500° is possible if the 
escape takes place at intervals when the temperature of the critical layers is raised to 2000° or 
more by the strong radiation of solar flares. 

The late R. T. Chamberlin, in a short and lucid presentation of the geological evidence, leaves 
open the question whether the earth was originally a single molten mass or was formed by ac- 
cretion, and he presents the evidence that the changes in climate and atmospheric conditions 
since early geologic time have been small and fluctuating. 

Next, Harrison Brown presents strong evidence that the inert ‘‘rare” gases are very much 


1 Pub. A.A.S., 10, 307, 1943. 2 At the meeting of the A.A.S., December, 1947. 
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rarer on earth than cosmically, which may be explained if the earth’s atmosphere is almost en- 
tirely of secondary origin—the result of chemical processes that took place since the formation 
of the planet. It may well be, as he suggests, that much of the terrestrial argon A*° may have 
been produced during geologic time by the radioactive decay of K*°. But this decay is a cosmic 
process which would increase the absolute abundance of argon, but not its relative abundance in 
different systems unless the initial abundance of K*° in these were different. Indeed, since free 
potassium atoms exist in interstellar space, while most of those on earth are locked up deep in 
the crust, the change might well increase the ratio of observable argon in the nebulae as com- 
pared with the earth. 

Three following papers deal with spectroscopic observations. Adel presents bolometric re- 
sults, including the first evidence for telluric VO and HDO; Migeotte the discovery of telluric 
CH,; and Dunham gives an admirable summary of the well-known photographic work at 
Mount Wilson. 

Kuiper’s ‘‘Survey of Planetary Atmospheres” is the most striking chapter in a notable book. 
Once more it exhibits its author’s remarkable flair for studying a problem thoroughly, predicting 
what new phenomena may be expected to be observable, and then observing and finding them! 

The study begins with insolation and escape-velocity and with estimates of the “deficiency 
factor” representing the ratio of the relative abundance on the given body to the general “‘cos- 
mical”’ values. (The known values for the earth range from unity for the metals to 10- for Ne 
and 10-4 for He.) With the resulting composition for the ‘“‘atmophile” portion of the body, the 
probable abundance of all plausible compounds is studied, with allowance for the chemical and 
photochemical reactions to which they may be exposed. 

The observations extend from A 3000 to A 27,000 (with a PbS cell). The resulting discoveries 
of CH, on Titan and CO; on Mars are well known, as is also the evidence that the polar caps of 
Mars are “composed of H,0 frost at temperatures much below 0° C,”’ while the rings of Saturn 
are also probably frost-covered, and that the green areas of Mars differ greatly in their reflecting 
power in the infrared from the higher plants but not from lichens. 

The reflecting power of the ‘‘desert’’ areas of Mars was found to differ widely from that of 
terrestrial red soil and rock but was closely matched by a “brownish, fine-grained felsite.” 
The conclusions that the higher regions of Mars are composed of similar rocks is somewhat star- 
tling. But the traditional ‘‘red”’ color of Mars is measured on the stellar scale in which sunlight is 
yellow. It might be interesting to set up some device whereby rock samples illuminated by a 
“daylight” lamp giving a distribution similar to sunlight in the visible region could be viewed at 
— = such a distance that they appeared as points about as bright as Mars and compared 

or color. 

New and better estimates of the amounts of gases in planetary atmospheres are given. The 
greatest values (in kilometers at S.P.T.) are 1.0 for CO, on Venus, and 1.5 and 2.5 for CH‘ on 
Uranus and Neptune—the last two much smaller than the early estimates. 

On the assumption that the surface temperature is that of a “gray” body exposed to solar ra- 
diation, it is impossible to explain how so much CH, can be observable in Uranus and Neptune, 
or NH; in Jupiter and Saturn. Kuiper, discussing Neptune, rejects, for good reasons, the hy- 
pothesis of a methane ocean and accepts that of an atmosphere of almost pure hydrogen, in 
which the “‘dewpoint” for formation of clouds of methane is reached at a depth exceeding 120 km 
and at a temperature which may be 40° K higher than that at the outer limit. Opacity of the 
atmosphere for long waves may diminish the observed radiometric temperature almost to the 
boundary value. 

It may be added that, on all four planets, the clouds are composed of frozen gas and that the 
possibility of influence of internal planetary heat still deserves investigation.’ 

The relation between gravity and pressure leads to a small point which even Kuiper has 
missed. Surface-gravity on Titan is one-seventh that on Saturn. Hence, for an atmosphere of 
pure methane, seven times the “thickness” could exist without condensation on Titan—other 
things being equal. The actual amount is less than on Saturn and could not condense on the sun- 
lit side even if no hydrogen were present—though there may be condensation on the dark side. 

Some problems, such as the composition of the cloud-layer on Venus, remain unsolved. The 
laboratory absorption spectra, next described by Herzberg, were obtained with paths up to 5 km 
in length, using multiple reflections in long tubes and low gas pressures, and provided funda- 


3 Compare the discussion of the older data (Russell, Dugan, and Stewart, Astronomy, 1 (rev. ed.; 
Boston: Ginn & Co., 1945), 367, 383, 401. 
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mental data for Kuiper’s work. Forbidden (quadrupole) absorption lines of H; should appear 
in high-dispersion spectra of Uranus and Neptune if these atmospheres are deep enough. 
A short note by Franck on the possibility of photosynthesis on Mars—leaving the question 
open for the present—concludes a volume of remarkably high and varied interest. 
HENRY Norris RUSSELL 
Princeton University Observatory 


ERRATA 


In Volume 107, page 287, 1948 (Babcock, Moore, and Coffeen, ‘The Ultraviolet Solar Spec- 
trum, AA 2935-306,” Mt. W. Contr., No. 745), the undersigned introduced an error by applying 
+0.001 A for reducing to standard atmospheric conditions. The correction should have been 
—0.001 A. All wave lengths in Tables 1 and 3 should be diminished by 0.002 A. 

In Table 2 the first line reads: 

it should read: 

A new second line should be added, reading: 

3133-3370 A (—0.0010)f, 


and the footnote concerning AX 3133-3370 should be deleted. 
Haro tp D. BABcock 


The legend for Figure 1 of Nassau and van Albada, “A Study of M-Type Stars in Cygnus” 
(Ap. J., 109, 393, 1949) is defective; it should read as follows: 
Fic. 1.—Principal features in infrared spectra of M stars. a, @ Ser (K4); 6, 83 UMa (M2); 


c,a Her A (M5); d, BD+49°2999 (M7); e, U Her, April 7, 1947 (M7) shows band near A 7900; 
f, S UMa, April 8, 1947 (S2e) shows LaO band; g, Y CVn, April 8, 1947 (N3 or C54). 


The following paragraph should be added to the note by Miller and van Dien, “A Large New 
Planetary Nebula” (Ap. J., 109, 537, 1949): 

“We are indebted to Dr. Minkowski for informing us that he learned of this object by letter 
from the Harvard College Observatory in 1941, that it was found on the Metcalf plate by Miss 
Rebecca Jones, and that no description of it has been published. The data presented above may 
therefore be of some interest.” 
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BY RALPH B. BALDWIN 


An answer, in terms of the meteor- 
ite theory, to the question of how 
the moon came to exist in its pres- 
ent form. 

The study of the moon—a mirror 
for the study of the earth. 


256 pages. lUlustrated. $5.00. 
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ASTROPHYSICAL MONOGRAPHS 


Sponsored by the Editors of the Astrophysical Journal — 


PRINCIPLES OF STELLAR DYNAMICS 
By S. CHANDRASEKHAR 

spiral nebulae, and star clusters are developed from a coherent point of view. “This 
book should exert a profound influence on the future developments in the field of 
galactic dynamics. I can recommend its study unreservedly to newcomers in the field 
and to those who already have a passing acquaintance with its problems. The experts 
can profit from reading it.”—Barr J. Box, in Science. 1942. $5.00 


SPECTRA OF LONG-PERIOD VARIABLE : STARS 
By PAUL W. MERRILL 


discussion of the physical significance of various spectral features is included. A brief 
introduction and a historical chapter are inserted for the benefit of the reader who is 
not a specialist in stellar spectra, “To one wishing an intimate view of our present 
knowledge of the spectra of long-period variable stars, whether he be a student, as- 
tronomer, physicist, or chemist, the present monograph is recommended.” J. H. 
Moore in Astronomical Society of the Pacific. 1940. Paper, $2.50 


THE MASSES OF THE STARS, WITH A GENERAL 
CATALOGUE OF DYNAMICAL PARALLAXES 
By HENRY NORRIS RUSSELL and CHARLOTTE E. MOORE 

“Astrophysicists will be grateful to the authors for so thorough an examination of ~ 
available knowledge. . .. We may be permitted to congratulate the senior author on 
this magnificent outgrowth from an idea which germinated in his mind. thirty years 
ago.”-—Prorrssor A. S. Eppincton, The 1940. Paper, $4.50 


THE DISTRIBUTION OF THE STARS IN 5 cai 
By BART J. BOK 

“. . . an extremely important contribution to the field of galactic structure, invaluable 
alike to both students and research workers in the field.” —Science. The large number 
and variety of researches carried out in this field during the past fifteen years are dis- 
cussed critically in this study with clarity and dispatch. 1937. Paper, $2.50 
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